Ising n-fold integrals as diagonals of rational 
functions and integrality of series expansions: 
integrality versus modularity 

November 26th, 2012 

A. Bostan^, S. Boukraa||, G. Christol*, S. Hassani§, J-M. 
Maillard^ 

^ INRIA, Batimcnt Alan Turing, 1 rue Honore d'Estienne d'Orves, Campus de 
I'Ecolc Polytechnique, 91120 Palaiseau, France 

j I LPTHIRM and Departement d'Aeronautique, Universite de Blida, Algeria 
1 Institut de Mathematiques de Jussieu, UPMC, Tour 25, 4eme etage, 4 Place 
Jussieu, 75252 Paris Cedex 05, France 

§ Centre de Recherche Nucleaire d'Alger, 2 Bd. Frantz Fanon, B.P. 399, 16000 
Alger, Algeria 

^ LPTMC, UMR 7600 CNRS, Universite de Paris 6, Tour 23, 5cmc ctage, case 
121, 4 Place Jussieu, 75252 Paris Cedex 05, France 

Abstract. 

We show that the n-fold integrals x^"' of the magnetic susceptibility of 
the Ising model, as well as various other n-fold integrals of the "Ising class", 
or n-fold integrals from cnumcrative combinatorics, like lattice Green functions, 
are actually diagonals of rational functions. As a consequence, the power 
series expansions of these solutions of linear differential equations "Derived From 
Geometry" are globally bounded, which means that, after just one rcscaling of 
the expansion variable, they can be cast into series expansions with integer 
coefficients. Besides, in a more enumerative combinatorics context, we show that 
generating functions whose coefficients are expressed in terms of nested sums 
of products of binomial terms can also be shown to be diagonals of rational 
functions. We give a large set of results illustrating the fact that the unique 
analytical solution of Calabi-Yau DDEs, and more generally of MUM DDEs, is, 
almost always, diagonal of rational functions. We revisit Christol's conjecture 
that globally bounded series of G-operators are necessarily diagonals of rational 
functions. We provide a large set of examples of globally bounded series, or series 
with integer coefficients, associated with modular forms, or Hadamard product of 
modular forms, or associated with Calabi-Yau ODEs, underlying the concept of 
modularity. We finally address the question of the relations between the notion of 
integrality (series with integer coefficients, or, more generally, globally bounded 
series) and the modularity (in particular integrality of the Taylor coefficients of 
mirror map), introducing new representations of Yukawa couplings. 
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1. Introduction 

The series expansions of many magnetic susceptibilities (or many other quantities, 
hke the spontaneous magnetisation) of the Ising model on various lattices in arbitrary 
dimensions are actually series with integer coefficients [1, 2, 3]. This is a consequence 
of the fact that, in a van der Waerden type expansion of the susceptibility, all 
the contributing graphs are the ones with exactly two odd-degree vertices and the 
number of such graphs is an integer. When series expansions in theoretical physics, 
or mathematical physics, do not have such an obvious counting interpretation, the 
puzzling emergence of series with integer coefficients is a strong indication that some 
fundamental structure, symmetry, concept have been overlooked, and that a deeper 
understanding of the problem remains to be discovered! . Algebraic functions are 
known to produce series with integer coefficients. Eisenstein's theorem [6, 7] states 
that the Taylor series of a (branch of an) algebraic function can be recast into a 
series with integer coefficients, up to a rescaling by a constant (Eiscnstein constant). 
An intriguing result due to Fatou [8] (see pp. 368-373) states that a power series 
with integer coefficients and radius of convergence (at least) 1, is either rational, 
or transcendental. This result also appears in Polya and Szego's famous Aufgaben 
book [9] (see Problem Vin-167). Polya [10] conjectured a stronger result, namely 
that a power series with integer coefficients which converges in the open unit disk is 
either rational, or admits the unit circle as a natural boundary (i.e. it has no analytic 
continuation beyond the unit disk). This was eventually proved^ by Carlson [12, 13]. 
Along this natural boundary line, it is worth recalling [14, 15, 16, 17, 18] that the 
series expansions of the full magnetic susceptibility of the 2D Ising model corresponds 
to a power series with integer coefficientsf . For them, a unit circle natural boundary 
certainly arises [19] (with respect to the modulus variable fc), but, unfortunately, this 
cannot be justified by Carlson's theoremff. 

A series with natural boundaries cannot be D-finite, i.e. solution of a linear 
differential equation with polynomial coefficients [20, 21]^. For simplicity, let us 
restrict to series with integer coefficients (or series that have integer coefficients up 
to a variable rescaling), that are series expansions of D-finite functions. Wu, McCoy, 
Tracy and Barouch [25] have shown that the previous full magnetic susceptibility 
of the 2D Ising model can be expressed (up to a normalisation factor (1 — sY^'^/s, 
see [16, 26]) as an infinite sum of n-fold integrals, denoted by x'"'': which are actually 

X The emergence of positive integer coefficients corresponds to the existence of some underlying 
measure [4] (see also the concept of Mahler measures [5]). 

^ The Polya-Carlson result can be used to prove that some integer sequences, such as the sequence 
of prime numbers [pn] [11], do not satisfy any linear recurrence relation with polynomial coefficients, 
f In some variable w [15, 14, 16, 17]. In the modulus variable k, one needs to perform a simple 
rescaling by a factor 2 or 4 according to the type of (high, or low temperature) expansions, 
ft The radius of convergence is 1 with respect to the modulus variable k, in which the series does not 
have integer coefficients, being globally bounded only (this means that it can be recast into a scries 
with integer coefficients by one rescaling of the variable fc). If one considers the series expansion with 
respect to another variable (such as ui) in which the series does have integer coefficients, then the 
radius of convergence is not 1. 

tt D-finite series are sometimes called holonomic. A priori, these notions differ: a function 
f{xi , . . . , Xr) is called D-finite if all its partial derivatives D"^ ■ ■ ■ Z?""" ■/ generate a finite dimensional 
space over Q(a;i, . . . , x^), and holonomic if the functions x^i ■ ■ ■ i^^'-Dfi ■ • ■ D^-- • / obtained by 
multiplying monomials in the variables and higher-order derivatives of / subject to the constraint 
oi + ■ ■ ■ +oir -|-/3i -|-- ■ ■ -|-/3r < N span a vector space whose dimension over Q grows like 0{N^). The 
equivalence of these notions is proved by profound results of Bernstein [22] and Kashiwara [23, 24]. 
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D- finite^. We found out that the corresponding (minimal order) differential operators 
are Fuchsian [14, 16], and, in fact, "special" Fuchsian operators: the critical exponents 
for all their singularities are rational numbers, and their Wronskians are iV-th roots 
of rational functions [27]. Furthermore, it has been shown later that these x*^")'s are, 
in fact, solutions of globally nilpotent operators [28], or G-operators [29, 30]. It is 
worth noting that the scries expansions, at the origin, of the x''"''s, in a well-suited 
variable [16, 26] w, actually have integer coefficients, even if this result does not have 
an immediate prooff for all integers n (in contrast with the full susceptibility). 

From the first truncated series expansions of x^"\ the coefficients for generic n 
can be inferred [28] 

X^")(w) = 2"-w"' • (l + 4n2.w2 + 2-(4n'* +13n2 +1)-m;'' 

+ ^ ■ {n^ +4){4:n^ +23n^ + 3) -w^ (1) 

+ - ■ {32n^ + 624 + 4006 + 8643 + 1404) • 
3 

+ ^ - + (32n^ +784n^ + 6238n^ + 16271 + 3180) • w^^ + 
15 / 

Note that the coefficients of the expansion of x(")(i(j)/2" depend on n^. Note 
that these coefficients are integer coefficients when n is any integer, this integrality 
property [31] of the coefficients becoming straightforward to see when one remarks 
that {4n^ +23n^ + 3) and (32 n« + 624 + 4006 + 8643 + 1404) are of the 
form n ■ in? — 1) f{n) + 3 g(n) (respectively f{n) = An and g(n) = +1, and 
f{n) = 2n- {16n^ +328n'^ + 2331) and g{n) = 4435ri2 +468), and, hence, are 
always divisible by 3, that (32 + 784 + 6238 n'^ + 16271 + 318O) is of the 
form n- (n^ - 1) (n^ - 4) • f{n) + 3 • 5 • g{n), with f{n) = 16 n • (2 + 59) and 
g{n) = 722 + 833 + 212, hence, always divisible by 3 and 5. 

These coefficients are valid up to for n > 3, w'^ for n > 5, for n > 7, 
uu^ for n > 9, and w^^ for n > 11 (in particular it should be noted that x^"-* is an even 
function of w only for even n). Further studies on these x^^^'s showed the fundamental 
role played by the theory of elliptic functions^ (elliptic integrals, modular forms) and, 
much more unexpectedly, Calabi-Yau ODEs [32, 33]. These recent structure results 
thus suggest to see the occurrence of series with integer coefficients as a consequence of 
modularity [34] (modular forms, mirror maps [32, 33, 34, 35], etc) in the Ising model. 

Along this line, many other examples of series with integer coefficients emerged 
in mathematical physics (differential geometry, lattice statistical physics, enumerative 
combinatorics, replicable functions]), ...). One must, of course, also recall Apery's 
results [44]. We give, in Appendix A, a list of modular forms, and their associated 
series with integer coefficients, corresponding to various lattice Green functions [45, 

fj For Ising models on higher dimensional lattices [1, 2, 3] no such decomposition of susceptibilities, 
as an infinite sum of D-finite functions, should be expected at first sight. 

t We are interested in this paper in the emergence of integers as coefficients of D-finite series. In 
general, this emergence is not obvious: it cannot be simply explained at the level of the linear 
recurrence satisfied by the coefficients, as illustrated by the case of Apery's calculations (see also 
Appendix G below). 

^ Which is not a surprise for Yang-Baxter integrability specialists. 

(t The concept of replicable functions [36] is closely related to modular functions [37, 38], (see 
the replicability of Hauptmoduls [38]), Calabi-Yau threefolds, and more generally the concept of 
modularity [34, 39, 40, 41, 42, 43] (the third etale cohomology of a rigid Calabi-Yau threefold comes 
from a modular form of weight 4, ...). 
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46, 47, 48]. This integrality is also seen in the nome and in other quantities like the 
Yukawa coupling [32]. 

We restrict to series with integer coefficients, or, more generally, globally 
bounded [49] series of one complex variable, but it is clear that this integrality property 
does also occur in physics with several complex variables: they can, for instance, 
be seen for the previous (D-finitc§) n-fold integrals x*-"-* for the anisotropic Ising 
model [50] (or for the Ising model on the checkerboard lattice), or on the example of 
the lattice Ising models with a magnetic field| (see for instance, Bessis et al. [4]). 

One purpose of this paper is to "disentangle" the notion of series with integer 
coefficients {integrality) and the notion of modularity [34, 39, 40, 41, 42, 43]. In 
this down-to-earth paper we will use the wording of "modularity", not to refer to 
the modularity conjecture and other Serre's results that certain Geometric Galois 
representations are modular, but as a quick proxy word to say that a series solution of 
a linear differential operator as well as the nome, and hopefully other series (Yukawa 
coupling, ...) are all series with integer coefficients. 

We will show that the ^'"•''s are globally bounded series, as a consequence of the 
fact that they actually are diagonals of rational functions for any value of the integer 
n. We will generalise these ideas, and show that an extremely large class of problems 
of mathematical physics can be interpreted in terms of diagonal of rational functions: 
n-fold integrals with algebraic integrand of a certain type that we will characterise, 
Calabi-Yau DDEs, MUM linear DDEs [52], series whose coefficients are nested sums 
of binomials, etc. We take, here, a learn-by-example approach: on such questions one 
gets a much deeper understanding from highly non-trivial examples than from general 
mathematical demonstrations [53, 54]. 

2. Series integrality 

2.1. Globally bounded series 

Let us recall the definition of being globally bounded [49] for a series. Consider a series 
expansion with rational coefficients, with non-zero radius of convergencef . The series 
is said to be globally bounded if there exists an integer TV such that the scries can be 
recast into a series with integer coefficients with just one rescaling x ^ N x. 

A necessary condition for being globally bounded is that only a finite number 
of primes occurs for the factors of the denominators of the rational number series 
coefficients. There is also a condition on the growth of these denominators, that must 
be bounded exponentially [49], in such a way that the series has a non-zero p-adic 
radius of convergence for all primes p. 

When this is the case, it is easy to see that these series can be recast, with just 
one rescaling, into series with integer coefficients^. 

It will be seen, in a forthcoming section (see (3) below), that the series expansion 
of diagonals of rational functions [55, 56, 57] are necessarily globally bounded. 

§ For several complex variables the ODEs of the paper are replaced by Picard-Fuchs systems. 

I Along this line original alternative representations of the partition function of the Ising model in a 

magnetic field are also worth recalling [51]. 

f A series like the Euler-series X]^o ' ^" which has integer coefficients is excluded. 

^ For a first set of series with integer coefficients, see Appendix A, where a set of such series with 

integer coefficients corresponding to modular forms is displayed. 
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2.2. Globally logarithmically bounded series 

There is another notion, weaker than being globally bounded, namely the notion of 
being globally logarithmically^ bounded. 

As an example consider the series expansion of 2^'i([l/4, 1/2], [5/4],4a;). This 
series is not globally bounded 
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When looking at the denominators of the series coefficients, one finds that almost all 
primes of the form A£ +1 occur. There is no way to recast this scries into a series 
with integer coefficients with one rescaling. 

Let us denote 6 = x-d/dx. The hypcrgcomctric function 2^1 ([1/4, 1/2], [5/4], x) 
is solution of the operator 

i) (3) 

which clearly factors f 40 + 1 at the right: 

(29- X ■ {29 + 1)) • {49 + 1). (4) 

Consequently, the action of 40 + 1 on 2-Fi([l/4, 1/2], [5/4], x) becomes the solution of 
the order-one globally nilpotent operator 29 —x ■ {2 9 + 1), and is, thus, an algebraic 
function (rational or iV-th root of rational), namely (1 — x)^^/^. The hypergeometric 
function 2^i([l/4, 1/2], [5/4], x) is not globally bounded. We will see below that, 
consequently, it cannot be the diagonal of a rational function, however it is not a 
"wild" series, the denominators do not grow "too fast": it is actually such that a 
simple order-one operator, namely 4 0-1-1, acting on this series, changes it into a 
diagonal of rational function. Other examples of globally logarithmically bounded 
hypergeometric series are given in Appendix B. 



3. Minimal recalls on diagonals of rational functions 

Let us recall here the concept of diagonal of a "function"^, and some of its most 
important properties. 

3.1. Definition of the diagonal of a rational function 

Assume that J^(zi, . . . , z„) = P{zi, . . . , z„)/Q{zi, . . . , Zn) is a rational function, 
where P and Q are polynomials with rational coefficients such that (5(0, . . . , 0) 7^ 0. 
This assumption implies that J-' can be expanded as a Taylor series at the origin 

t(^Zi, Z2, • ■ ■ , Zn^ = (5) 

00 00 00 

E E ■■■ E ^r^2™^ ■■■ e q[[^i,...,^„]]. 

mi — m2 — m„ — 

§ For a series ^a„x" being "bounded" means bounded by 1, p-adically: \an\p < 1, i.e. a„ has 
no p factor at the denominator.) Being logarithmically bounded means "with logarithmic grows", 
i.e. \an\p < n. 

t This result can, of course, straightforwardly be generalised to 2^1 ([a, ''l, [1 -I- <i]i a^). 

^ This is an abuse of language: the "functions" are in fact defined by series of several complex 

variables: they have to be analytical (no Puiseux series). 
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The diagonal of J- is defined as the series of one variable 

oo 

Diag(^(zi, Z2, z„)) = ^ ™- z" e QM. (6) 

m = 

More generally, one can define, in a similar way, the diagonal of any multivariate 
power series T S ^^[[zi, • ■ • , Zn]], with coefHcients in an arbitrary field K (possibly of 
positive characteristic)!. 

3.2. Main properties of diagonals 

The concept of diagonal of a function has a lot of interesting properties (see for 
instance [59]). Let us recall, through examples, some of the most important ones. 

The study of diagonals goes back, at least, to Polya [60], in a combinatorial 
context, and to Cameron and Martin [61] in an analytical context related to Hadamard 
products [62]. Polya showed that the diagonal of a rational function in two variables 
is always an algebraic function. The most basic example is T = 1/(1 — zi — Z2), for 
which 

Diag(^) = 




E 

m=0 



''"^ z™ = (7) 



The proof of Polya's result is based on the simple observation that the diagonal 
Diag(J^) is equal to the coefficient of z^ in the expansion of T{zi, z / zi). Therefore, 
by Cauchy's integral theorem, Diag(J^) is given by the contour integral 

Diag(^) = [z-^]T{z,,z/z,)/z, = -I^ /j-(zi,zM)^, (8) 

2m zi 

where [z"] meansf extracting the n-th coefficient of a power series, and where the 
contour 7 is a small circle around the origin. Therefore, by Cauchy's residue theorem, 
Diag(J^) is the sum of the residues of the rational function Q = J-"(zi, zj z\)l z\ at all 
its singularities s(z) with zero limit at z = 0. Since the residues of a rational function 
of two variables are algebraic functions, Diag(J^) is itself an algebraic function. 

For instance, when J" = 1/(1 — zi — Z2), then Q = T {z\, z j z\) j z\ has two 
poles at s = 5(1 ± -s/l — 4z). The only one approaching zero when z — > is 
So = ?(1 ^ \/l — 4z). If p{s)/q{s) has a simple pole at sq, then its residue at sq is 
p{so) / q' {so) . Therefore 

1 / dzi 



Zl — z 



dzi 1 1 

= Res^n o = = , (9) 

Zl - z2 - z 1 - 2so VT^^ 

When passing from two to more variables, diagonalisation may still be interpreted 
using contour integration of a multiple complex integral over a so-called vanishing 

X The definition even extends to multivariate Laurent power series, see e.g. [58]. 
f This is a convenient notation, very often used in combinatorics [63]. 
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cycle [64]. However, the result is not an algebraic function anymore. A simple example 
is = 1/(1 — Z2 — Z3 — z\Zi — 2123), for which 

Diag( J") = 1 + 4z + 362^ + 400^^ + 4900z'' + 63504z'5 + 8537762^ 

+ 11778624z^ + • • • (10) 

is equal to the complete elliptic integral of the first kind 

2 

• z"- 

\ m / 

)ri>0 

^""l"^ „ /A 1, 



1 - 16zsin^(i^) 



-], [l],16z), (11) 



which is a transcendental function. 

Less obvious examples (see [65]) are 

Diag (- ) (12) 

\ 1 - 21 - Z2 - 23 - 21 22 - 22 Z3 - 2321 - 212223 



1 „ /rl 2, 542 
• 2-Fl 



3' 3^' ' " (1 -2)3 



and 



Digg/" (l-2i)(l-22)(l-23) \ 

\1 - 2 (21 + 22 + 23) +3 (2122 + 2223 + 2321) - 4 212223 / 

Jo V3' 3^' (l-4u;)3 J ' {I ^ iw) {1 ^ 6iw) ' 

It was shown by Christol [66, 67, 68] that the diagonal Diag(J^) of any rational 
function T is D-finite, in the sense that it satisfies a linear differential equation with 
polynomial coefficients^. Moreover, the diagonal of any algebraic power series in 
Q[[2i, . . . , 2„]] is a G-function coming from geometry, i.e. it satisfies the Picard-Fuchs 
type differential equation associated with some one-parameter family of algebraic 
varieties. Diagonals of algebraic power series thus appear to be a distinguished class of 
G-functionsft. It will be seen below (see (3.5)) that algebraic functions with n variables 
can be seen as diagonals of rational functions with 2n variables. Thus diagonals of 
rational functions also appear to be a distinguished class of G-functions. It is worth 
noting that this distinguished class is stable by the Hadamard product: the Hadamard 
product of two diagonals of rational functions is the diagonal of rational function. 

An immediate, but important property of diagonals of rational functions in 
Q[[2i, . . . , 2„]] is that they are globally hounded, which means that they have integer 
coefficients up to a simple change of variable z ^ N z, where N ^ "L. 

Furstenberg [69] showed that if K has positive characteristic, then the diagonal 
of any rational power series in ii'[[2i, . . . , 2„]] is algebraic. Deligne [64, 58] extended 
this result to diagonals of algebraic functions. For instance, when F = 1/(1 — 22 — 
23 — 2i22 — 2123), one gets modulo 3, modulo 5 and modulo 7 respectively 

^ A more general result was proved by Lipshitz [57] : the diagonal of any D-finite series is D-finite, 
see also [56]. 

(t Such diagonals are solutions of G-operators. They are functions that arc always algebraic mod. 
a prime p. They fill the gap between algebraic functions and G-scrics: they can be seen as 
generalisations of algebraic functions. 
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Diag(J') mod 3 = 1 + z + + + + + z^^ + 

1 



mod 3, 



Diag(J") mod 5 = 1 + 4 x + + + + + z^" + 4z" + + 

1 

mod 5, 



^1 - z + 

Diag(J") mod? = 1 + 4 z + z^ + z^ + 4z'^ + 2z^ + 4z^ +••■ 

= ^, = mod 7. 

v^l + 4z + z2 + z3 

More generally, for any prime p, one has 

Diag(J") modp = P(z)^/(^-p) mod p (14) 
where the polynomial P{z) is nothing but [70, 71, 72] 

^'(^) = 2Fi([^, [1], 16z)' " mod p 

(p-l)/2 . _ 1 /0\ 2 

^ [ I • (16^)" mod p. (15) 

For instance, modulo 11, the polynomial (15) reads: 



-], [l],16z 



mod 11 = l+4z+3z^+4z'^ + 5z''+z'' mod 11. 



Interestingly enough, the polynomial modulo p 

(p-l)/2 / -, /o\ 2 

Pi ' ' 



n=0 ^ ^ 

is [73, 74], up to a sign (—1)'^"^^/^, the 77asse mwarianif of j/^ = x- (1 — a;) - (A — cc). 

Note, however, that the Furstensberg-Deligne result [69, 64], that we illustrate, 
here, with J- ~ 1/(1 — Z2 — z^ — z\Z2 ~ z\z-i)^ goes far beyond the case of 
hypergeometric functions for which simple closed formulae can be displayed. 

3. 3. Hadamard, and other products 

Let us also recall the notion of Hadamard product [62, 75] of two series, that we will 
denote by a star. 

oo oo 

If fix) ^ ^a„-a;", gix) ^ ^ 6„ • x", then: 

n=0 n=0 

oo 

f{x)*g{x) = ^ x". (17) 

n=0 

f Note Igusa's sentence "Hence the elliptic differential of the first kind has only one period, and that 
is A{X), up to an arbitrary differential constant. This version of Hassc invariant has not yet been 
explicitly remarked." 
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The notion of diagonal of a function and the notion of Hadamard product are 
obviously related: 

Diag(^/i(a;i) • /2(x2) • • • /„(2:„)) = fi{x) ★ /2(x) ★ • • • ★ /„(x). (18) 

In other words, the diagonal of a product of functions with separate variables is equal 
to the Hadamard product of these functions in a common variable. In particular, 
the Hadamard product of n rational (or algebraic, or even D-finite) power series is 
D-finite. 

The Hadamard product of two series with integer coefficients is straightforwardly 
a series with integer coefficients. Furthermore, the Hadamard product of two operators, 
annihilating two series, defined as the (minimal order, monic) linear differential 
operator annihilating the Hadamard product of these two series, is a product compatible 
with a large number of structures and concepts that naturally occur in lattice statistical 
mechanics. For instance, the Hadamard product of two globally nilpotent [28] 
operators is also globally nilpotent. 

Let us introduce another product, namely the Hurwitz (shuffle) product] of two 
series which is defined as [77, 78, 76]: 



Hurwitz Prod ( a„ • x", ^ /3„ • x" j = X] X] 



n + m 
n 



A very simple example is, for instance, 

HurwitzProd(- , \ = . (19) 



Vl— a-x' 1 — b ■ x' 1 — {a +b) ■ X 

Again, we have a remarkable compatibility property between the diagonal and the 
Hurwitz product. The Hurwitz product of two series that are diagonals of two power 
series A{xi, X2, ■ ■ ■ , x„) and B{yi, y2, . . . , ym) can itself be seen as the diagonal of 
a power series [68], that is very close to the product of these two power series^: 

HurwitzProd(D\ag{A{xi, X2, x„)), Diag(B(yi, j/2, . ■ . , J/m))) (20) 

/A{xi, X2, . . . , Xn) ■ B{yi, y2, . . . , ?/,„) 
V I - t ■ xiX2 ■■■ x„yiy2 ■■■ ym 
where t is an additional variable. In fact there exists an infinite number of products 
that enjoy a similar property of compatibility with the diagonal. The most general 
products of series compatible with the diagonal are, beyond the Hadamard and 
Hurwitz products: 



GeneralProd 



n n n m 

where the p{n^ m)'s are coefficients of any rational function of two variables R{x^ y): 
Rix,y) = ^^p(n,TO). x^y". (21) 

n ra 

Special cases are Lamperti's product [79, 76] and Trjitzinsky's product [80]. 

f The Hurwitz product of two algebraic functions is not algebraic in general, e.g. the Hurwitz square 
of (1 - 42)-l/2 jg equal to 2-Pi{[l/2, 1/2], [1], 162 (1 - 42;)). However, the Hurwitz product of two 
algebraic functions is actually an algebraic function modulo a prime p, cf. Prop. 8 in [76]. The 
Hurwitz product of a rational and of an algebraic power series with coefficients in Q is algebraic, cf. 
Prop. 3 & 7 of [76]. 

I Note a misprint in the second equation of page 68 of [68] . 



Diagonals of rational functions 
3.4- Chiral Potts examples 

Let us consider, for instance, the Hadamard cube of a simple algebraic function 
3F2 
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1 1 1 



[1, 1] ; ^) = Diag((l - x)-i/3 ^ (1 _ y)-l/3 ^ (1 _ ^)-i/3^ 

iFo([i], []; z) * iFo([i], []; z) * iFo([i], []; z) 
(1 - * (1 - z)"i/3 ^ (1 _ ^)-i/3^ 



(22) 



It is globahy bounded 
"1 1 1 

3' 3' 3 



3F2 



1, 1] ; 35 



= 1 +9a; +648a:2 +74088x3 + 10418625 x'' 

6 



+ 1648059777 x^ + 281268868608 x^ + 50621016116736 x^ + ••• 

Other examples, related to the chiral Potts model and its associated Format 
curves [59] are 2Fi([l/2, 7V/3], [1], 36t) and 2Fi([l/2, 7V/5], [1], lOOt) which have 
series expansions with integer coefficients, or, more generally: 

t q 



TV' N' N 



= (1 - x)-*/^ * (1 - x)-"'" * (1 - x) 



q/N 



-s/N 



3.5. Furstenberg's result on algebraic functions 

It was shown by Furstcnberg [69] that any algebraic series in one variable can 
be written as the diagonal of a rational function of two variables (however, this 
representation is, by no means unique). For instance, 

/ = = X +-x^ +-x^ +—X* +—x' +— x6 + ... (23) 

2 8 16 128 256 ^ ' 

is the diagonal of (2 x y — cx + cy) /(x + y +2) for any rational number c. 

The basis of Furstenberg's result is the fact that if /(x) is a power series without 
constant term, and is a root of a polynomial P{x, y) such that Py{0, 0) 7^ 0, then 

When Pj,(0,0) — 0, this formula is not true anymore. For instance, it does not 
apply to the algebraic function / = .x/ \/l — x, annihilated by P = (x — l)y'^ + x^, 
since the diagonal of y'^ Py{xy,y)/ P{xy,y) = 2y {xy — l)/(x^ + xy — 1) is zero. 
However, Furstenberg's result still holds. A way of seeing this on our example is to 
observe that <? = / — x — ^x^ is an algebraic series annihilated by a polynomial 
Q such that (5(0,0) = and Qy{0,0) 7^ 0. This reasoning, which extends to the 
general situation, yields the following rational function whose diagonal equals /: 

xy ■ where: (25) 

Q[x, y) 

V{x,y) = 16x^2/^ +4. (3x-4) . x^y'' +4 . (3 + x) . x2y3 

+ (12x-24 + x2) X2/^ +5y.x2 +6x -16, 

Q(x, y) = 8x^y3 + 8 . (x - 1) . xt/^ + 2 . (x + 4) . xy +6x -16. 

When compared to (2 x ?/ — cx + cy)/{x +y +2), whose diagonal is also /, this shows 
that Furstenberg's proof does not necessarily produce the easiest rational function. 
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Furstcnbcrg's result has been generalised to algebraic power series in an arbitrary 
number of variables: any power series % in Q[[a;i, ...,a;„]] algebraic over Q(a;i, ...,a;„) 
is the diagonal of a rational function with 2n variables (see Denef and Lipshitz [81]). 

4. Selected n-fold integrals are diagonals of rational functions 

Among many multiple integrals that are important in various domains of mathematical 
physics, let us consider the 7i-particle contribution to the magnetic susceptibility of 
the Ising model which we denote They are given by (n — l)-dimensional 

integrals [14, 82]: 

^ Mi(nf^)-.-i±| 

where, defining $o by J27=o ~ 

n — 1 2 n— 1 



G= n ^'ff^'^'p^'^ ^^"^"^ A, = 1 -2z.cos($,). (27) 

The integrality property (1) had been checked [15] for the first x'""* and inferred [28] 
for generic n. We are going to prove itX for any integer n, showing a much fundamental 
result, namely that all the {n^l)-fold integrals x*""^ 's are very special: they are actually 
diagonals of rational functions. 

4-1- X^^^ '^'5 ^ ioy example 

At first sight the x*-"^'s are involved transcendental holonomic functions. Could it be 
possible that they correspond to the distinguished class [58] of G-functions, generalising 
algebraic functions, which have an interpretation as diagonals of multivariate algebraic 
functions (and consequently diagonals of rational functions with twice more variables)? 
If this is the case, then the series of the x^"'''s must necessarily reduce modulo any 
prime to an algebraic function (see (??))■ The x''^'' and x''^'' contributions being too 
degenerate (a rational function and a too simple elliptic function), let us consider the 
first non-trivial case, namely x^^''- Its series expansion has already been displayed 
in [14]. It reads x^^V^ = w"^ ■ F{w) with: 

F{w) = 1 + 36^2 + Aw^ + 884 w^'"^ + 196 + 18532^6 + 6084 + ••• 

Since we have obtained the exact ODE satisfied by x^'^'' '^e can produce as many 
coefficients as we want in its series expansion. Let us consider this series modulo the 
prime p = 2. It now reads the lacunary series 

F{w) mod 2 = 1 + + + w''^ + w^^" + w^^^ + w^"^ + w^^^^ +■■■ , 

^ In the one-variable case, Puiscux scries could be considered but only after ramifying the variable. 
J Actually we only prove global boundedness for the Taylor expansion x^"\w) = a^io''. 
However, looking et the process more carefully, and, in particular, adding corresponding properties 
on the set Tn below, one can find out that only powers of 2 appear in the denominators of the aj, : 
the rescaling factor ( "Eisenstein constant" ) is 2 or 4 according to the fact that one considers high or 
low temperature series [15, 27], and that J^afeUi'^ do converge for |ui| < 1/4. 
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solution of the functional equations on F{w) or, with z = w^, on G{z) = 
1 + ■ F{w) 

F{w) = 1 + ■ F{w'^), G{z) = z + G(z2), (28) 

where one recognises, with equation G{z) = z + G{z^), Furstenberg's example [69] 
of the simplest algebraic function in characteristic 2§. In fact H{w) = F{w) is 
solution of the quadratic equation: 

H{wy + w ■ H{w) + u)i" = mod 2. (29) 

The calculations are more involved modulo p = 3. Indeed, Hiw) ~ x'^^^(u')/8 
satisfies, modulo 3, the polynomial equation of degree nine 

P9 • H{wf + ■ P3- H{wf + w^" • pi ■ H{w) 

+ w'' ■ pi'^ ■ p'^^ = 0, (30) 

where: 

P3 = (u;2 + 1)18 (1 „ u;)15 (y;4 _ ^2 _ ^ (w2 + l)20(i_^)13^ 

p\^'^ = +W'^ ~V? ~W + l, (31) 

= u;37 _ ^36 ^ ^35 _ ^33 ^ ^31 _ ^30 ^ ^28 ^ ^27 ^ ^24 _ ^23 ^ ^22 
- - _ ^16 ^ ^14 _ ^12 _ ^11 _ ^10 + _ ^5 _ ^3 _ 

The calculations are even more involved modulo larger primes. The scries for 
T^(3) niod. 5 reads: 

^(3) = ur^ + w^^ +Aw^'^ +Aw^^ + w^^ +2w^^ +Aw^^ + ■■■ (32) 

The (minimal order) linear differential operator annihilating the x^'^^ series mod. 5, 
reads f: 

{x + l){x^ + x + l){x + 2) - x^ ■ Dl +2x^- (x^ + 2x^ + 4x + A)ix + A) ■ 

+ • (x* + 3 a;^ + 4) • + 4 • (x^ + 3)- x - +3 (33) 

If one can easily get this linear differential operator, finding the minimal polynomial 
of x^^^ modulo 5, generalising (29) or (30), such that P{x^^\w), w) = mod. 5, 
requires a very large number of coefficients. Since the series (32) starts with , it is 
more convenient to consider the polynomial P{k, w), relating k = x'^'^^ (m;)/w^ and 
w. This (minimal) polynomial! is a polynomial of degree 50 in k, and degree 832 in 
w, sum of 4058 monomials. This (minimal) polynomial of the form: 

+ P!f'"\w) ■ + P^'''\w) ■ K + P^'''\w), (34) 

where the Pn"^\wys are polynomials of degree m in w, and where the head 
polynomial reads: 

d(8 
50 



P5(832) ^ ^382 . ^ 2)20 {w^ + 2 W + 4)^^ 1)70 ^ 4)20 



X (u;^ +3w + 4)'^(w'*+4w^ + w + l)'". (35) 

§ Modulo the prime p = 2, the previous functional equation becomes G(z) = z + G{z)^. 
f This operator is of zero 5-curvature [28] . 

I This polynomial has been checked with a series (32) of 380000 coefficients. 
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This (minimal) polynomial is a factor of a much larger polynomial (in 1, k, k^, 
, and ) of a more "p-adic nature" , which is of degree 125 in k, sum of 3559 
monomials, and of degree 1941 in w. 

One can imagine, in a first step that the x^^'^ series mod. any prime p are also 
algebraic functions, and, in a second step, that x^^"^ niay be the diagonal of a rational 
function. In fact we are going to show, in the next section, a stronger result: the 
^(ra) 'g g^j.g actually diagonals of rational functions, for any integer n. 

4-2. The X^"'' 's ire diagonals of rational functions 

Let us, now. consider the general case where n is an arbitrary integer. 

With the change of variable Zi = exp(t(f>i) (where = —1), one clearly gets 

Y[z^ = 1 . — = id^j , (36) 



1=0 



1 z, 
2cos($i) = Zi -\ , 2cos($i - (f>j) 



Zi Zj Zi 



n-1 



and (26) becomes 

^H(^) = (fl ±1^). Fiw,z„...,z^.,) , (37) 

n! V ZLTT Jq Zj / 

where C is the path "turning once counterclockwise around the unit circle" and where 
F is algebraic over Q(w, zi, . . . , Zn-i) and reads: 

F(«;, Zl,..., z„_i) = r.X"-i.i^.G. (38) 

Now, let us suppose that F is analytic] at the origin, namely that it has a Taylor 
expansion (5). Then applying {n — 1) times the residue formula, one finds 

X^''\w) = Diag( ^^^^'', ■ f(zi...z„_iz„, Zl, z„_i)). (39) 

To check that this is actually true, we introduce an auxiliary set, namely %i the subset 
of Q[zi, . . . , z„_i, zf ^, . . . , z~-^][[u;]], consisting of series 

oo 

/(W, Zl, . . . , Z„_i) = ^ P,„ • W™, 

m=0 

where Pm belongs to Q[zi, . . . , z„_i, z^^, . . . , z~^] and is such that /(zi, • • • , z„) 
belongsto Q[zi, z„_i][[z„]] C Q[[zi, z„]], where /is defined by: 

/(zi, •••,z„) = /(zi---z„, Zl, ••• ,z„_i) (40) 

In other words, we ask the degree of Pm, in each of the z~^, to be at most to. 
Then to prove that F has a Taylor expansion, we only have to verify that F belongs 

f One could consider Laurent, instead of Taylor, expansions, but this is a slight generalisation. The 
rational function l/{x + xy) would become allowed but rational functions like l/{x + y) would remain 
forbidden. For similar purposes, B. Adamczewski and J. P. Bell [58] recently used a generalised notion 
of Laurent expansions in the several variables case due to Sathaye [83] , see also [84] for various other 
generalisations . 
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to %i- Checking this is a straightforward step- by-step computation on auxihary 
functions: 

Ai = 1 — w ■ {zi -\ ), for: 1 1, 

A, ^ 1 -w-( \- zi ■ ■ ■ z„-i) , for: z= 0, (41) 

Ai = 1 -zi ■■■ • Zi+iZn ■ {zf -1-1), for: 1 < i < n-1, 

Ai = 1 - z„ • (1 + zi • • • 4-0 > for: i ^ 0. 

Hence Ai e Tn- The set Tn being clearly a Q[w]-algebra, A^ — e Tn- But 
%i is complete for the w-adic valuation. In particular, it is stable by the operations^ 

f{w) ^ = I -w f{w) + • • • , and: (42) 

\ +w ■ j(w) 

1 



fiw) + w f{w) = 1 +-wf{w) + ••• (43) 

So to be sure that the inverse or the square root of some function in Tn is also in Tn 
we have only to check that its first Taylor coefficient is actually 1: 

A^-Aw^ = 1 -2w {z, +-) ■ {z, - -)\ (44) 



hence yAf-4u;2 = 1 + ••■ ^ Tn, 

= , . ^ = I + ■■■ e Tn, ¥=!+■■■ e Tn, 



^JA\ - 4u;2 
2w 



Tji 1 *^ i j — ^ . . . ^ "J^ 



A, + ^Aj - Aw 



2 



n , 



x= + ••• G r„, = 1 +■■• e r„, 

G = n n J'' "5^' = n ^^"^ + •■■ e r„. 

0<z<j<n-l ^'^3) ■ 0<»<j<n-l ^'^^ 

From the definition of $o which implies JlILo^ = 1 wc also have 

n — 1 T 
/ ■ \ n — 1 

n = (n^o = ^^^^ 

0<i<j<n-l 1=0 

enabling to rewrite G in other ways. 

Thus, F belongs to Tn and it makes sense to take its diagonal. As F is algebraic, 
X^"^ is the diagonal of an algebraic function of n variables and, consequently, the 
diagonal of a rational function of 2 n variables. 

We thus see that we can actually find explicitly the algebraic function such that 
its diagonal is the n-fold integrals x^"'' • it is nothing but the integrand of the n-fold 
integral, up to trivial transformations, namely (38). 

Remark : x'"^ is a solution of a linear differential equation, and has a radius 
of convergence equal to 1/4 in w. Among the other solutions of this equation, there 

J More generally it is stable by the operations / — >■ (1 -|- ui f)^ for 5 G Q. 
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is the function obtained by changing the radical appearing in Xi into its opposite. A 
priori there are 2" ways to do this, hence 2" new solutions but, not all distinct. At 
first sight, for these new solutions, the x,;'s are no longer in Tn- 

In fact, we find some quite interesting structure. Let us consider, for instance, the 
case of x^'^''- If one considers other choices of sign in front of the nested square roots 
in the integrand, the series expansions of the corresponding n-fold integrals read: 

w + 6w^ + 28w^ + 124 w* + 536 + 2280 + 9604 + 40164 

+167066 + 692060 + 2857148 + • • • 

+6w^ + 30w^ + 140 + 628 + 2754 + 11890 + 50765 

+214958 + 904286 u)" +••• (46) 

These two series expansions (46) are solutions of the same order-seven operator [14] 
1/7 as x^^'. 

We know that other forms (equivalent for integration purposes) of G exist 
(see [85, 86, 87]). For these forms, other choices of sign in front of the nested square 
roots give 

S{+,-,+) = w + 6^2 + 28w^ + 126w'^ + 552w^ + 2388w^ + 10192 u-^ 
+ 43238 + 181936 + 762836 + 3180800 + • • • 

S{+,+,-) = w-2w2-20w^-110w^- 552w^- 2536w^- 11428w^ 

- 49898 - 216016 - 920776 - 3905764 + • • • 

S{-,-,+) = w +6w^ + 8 + Uw-^ ~8iw^ - 596w^ - AOOAw^ -19610 

- 99148 - 447332 - 2068492 + • • • 
S{-,+,-) = w +10^2 + 44^3 + 202^4 + 848^^ + 3672 wS + 15200 

+ 64310 + 264424 + 1104872 + 4523656 + • • • 

S{+,-,-) = u;^ +3u>3 + 13u;4 + 47ii;^ + 189w6 + 707u;'' + 2800w® 
+ 10637 + 41865 + 160535 + • • • 

S'(-,+,+) = +4 w^ + 25 w^ + 103 w^ + 496^^ + 2042 10^ + 9013 
+ 36931 +••• 

to be compared with x^^^ ~ S{+,+,+) — S{— ,—,—). These alternative series 
are not solutions of the order-seven operator [14] annihilating x^^\ however, 
some linear combination of these series are solutions of Lj. For instance the 
linear combination S{+, +, — ) — 5(— , +, +) — S{+, — , +), which is actually equal to 
S{-, +, -) + 4 S{+, -, -) - 5(-, -, +), is solution of Lj. 

Another form of G (again equivalent for integration purposes) gives: 

S{+,-,+) = u; +8^2 + 36 u)3 + 164 10"* + 704 w^ + 3041^^ + 127861/;^ 
+ 54067 + 224864 + 939709 + 3881708 ui" + • • • 

5(+,+,-) = w +2^2 -6^3 -48^^-314^5 - 1555w^- 7626w^ 

- 34461 - 155898 - 678199 - 2957648 ui" + • • • 
S'(-, +, +) = +3 + Uw^ + 51 + 214 + 810 + 3296 

+ 12679 + 50878 + 197466 + • • • 
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Again, these other ahernative series are not sohitions of Lj, but the hnear combination 
2 S{-,+, +) + +) - S{+, +, -) is solution of Lj. 

Ah these ahernative series are, in faet, solution of a higher order linear differential 
operator that iy rightdivides. 

One does remark that all these alternative series are, as x'^^K series with integer 
coefficients. 



4-3. More n-fold integrals of the Ising class and a simple integral of the Ising class 

It is clear that the demonstration we have performed on the x'^"''s can also be 
performed straightforwardly, mutatis mutandis, with other n-fold integrals of the 
"Ising class" like the ri-fold integrals in [82] , which amounts to getting rid of the 
fermionic term G (see (41)), the corresponding to n-io\d integrals associated 

with the diagonalf susceptibility [18, 33] (the magnetic field is located on a diagonal 
of the square lattice), the '&^'''s in [17] which are simple integrals, and also for all the 
lattice Green functions displayed in [48, 52], and the list is far from being exhaustive. 

in) 

For instance, the simple integral q is the diagonal of the algebraic function: 



4 • (1 - t^r^'^ - , where: (47) 

F„ = 1 -2w + (1 -4w +4^2 -4^2^2)1/2^ (48) 
G„ = 1 -2wt-T„J\) +((1 -2wt- T„_j\)f -Aw'' ■ t^^'^^ 



.t. 

and where Tn-i(t) is the (n — l)-th Chcbyshev polynomial of the first kind. The way 
we have obtained these Chebyshev results (47) is displayed in Appendix C. 

The integral ^^^\w) is the diagonal of an algebraic function of two variables, 
and also the diagonal of a rational function of four variables, and this independently 
of the actual value of n. 



4.4- More general n-fold integrals as diagonals 

More generally the demonstration we have performed on the x'^"^'s can be performed 
for any n-fold integral that can be recast in the following form: 

— — ••• — • ^(x, zi, Z2, • • • , z„), (49) 



where the subscript C denotes the unit circle, and where A denotes an algebraic 
function of the n variables, which (this is the crucial ingredient), as a function of 
several variables x and the z^'s, has an analytical expansion at {x, zi, Z2, ■ ■ ■ Zn) ~ 
(0, 0, 0, ••• , 0): 



A(2 



X, Zi, Z2, ■■■ , Zn] = (50) 



C« CXD CX3 OO 

V V V • • • V /I . ma 

m — mi — m2 — rUn — 



X 



f Of course this "diagonal wording" should not be eonfused with the notion of diagonal of a function. 
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Consequently, an extremely large set of n-fold integrals occurring in theoretical 
physics (lattice statistical mechanics, enumerative combinatorics, number theory, 
differential geometry, ...) can actually be seen to be diagonal of rational functions. 
Consequently these n-fold integrals correspond to series expansions (in the variation 
parameter x) that are globally bounded (can be written after one rescaling into series 
with integer coefficients), and are solutions of globally nilpotent [28] linear differential 
operators. 

Such a general n-fold integral is, thus, the diagonal of an algebraic function (or 
of a rational function with twice more variables [81]) which is essentially the integrand 
of such n-fold integral. Furthermore, such a general n-fold integral is solution of a 
(globally nilpotent) linear differential operator, that can be obtained exactly from the 
integrand, using the creative telescoping method (see Appendix D). 

Finally, in the case of Calabi-Yau ODEs (see below), these functions can 
be interpreted as periods of Calabi-Yau varieties, these algebraic varieties being 
essentially the integrand of such n-fold integrals. The integrand is thus the key 
ingredient to wrap, in the same bag, algebraic geometry viewpoint, differential 
geometry viewpoint and analytic and arithmetic approaches (series with integer 
coefRcients). 



5. Calabi-Yau ODEs 



Calabi-Yau ODEs have been defined in [88] as order-four linear differential ODEs 
that satisfy the following conditions: they are maximal unipotent monodromy [89, 90] 
(MUM), they satisfy a "Calabi-Yau condition" which amounts to imposing that the 
exterior squares of these order-four operators are of order five (instead of the order 
six one expects in the generic case), the series solution, analytic at x = 0, is globally 
bounded (can be reduced to integer coefficients), the series of their nome and Yukawa 
coupling are globally bounded]:. In the literature, one finds also a cyclotomic condition 
on the monodromy at the point at cxj, a; = oo, and/or the conifoldf character of one 
of the singularities [92]. 

Let us recall that a linear ODE has MUM (maximal unipotent monodromy [32, 
91]) if all the exponents at (for instance) x = are zero. 

In a hypergeometric framework the MUM condition amounts to restricting to 
hypergeometric functions of the type „+iF„([ai, 02, • • • a„], [1, 1, • • • l];a;), since the 
indicial exponents at a; are the solutions of p (p + 61 — 1) • • • (p -f- &„ — 1) = 
= 0, where the bj arc the lower parameters which arc here all equal to 1. 

Let us consider a MUM order-four linear differential operator. The four solutions 
yo, 2/1: y2, ys of this order-four linear differential operator read: 

ln(x)^ 

yo, yi = yo- ln(a;) -I- yi, 2/2 = yo ^ h yi ■ ln{x) + 2/2, 

ln(a;)^ ln(x)^ , , , 

ys = yo — ^ — + yi — ^ — + y2 ■ M^) + vs, 

where yo, TJi, 2/2, 2/3 are analytical at x — (with also yi(0) = ^2(0) = 2/3(0) = 0). 

X The instantons numbers are integers. 

f The local exponents are 0,1,1,2. For the cyclotomic condition on the monodromy at 00 see 
Proposition 3 in [91]. 
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The nomc of this hncar differential operator reads: 

q{x) = expf^) = s-expf^V (51) 

Calabi-Yau ODEs have been defined as being MUM, thus having one solution 
analytical at a; = 0. As far as Calabi-Yau ODEs are concerned, the fact that this 
solution analytical at a; = has an integral representation, and, furthermore, an 
integral representation of the form (49) together with (50), is far from clear, even if 
one may have a "Geometry-prejudice" that this solution, analytical at a; = 0, can be 
interpreted as a "Period" and "Derived From Geometry" [29, 30, 93]. 

Large tables of Calabi-Yau ODEs have been obtained by Almkvist et al. [91, 
94, 95]. It is worth noting that the coefficients An of the series corresponding to 
the solution analytical at a; = 0, are, most of the time, nested sums of product 
of binomials, less frequently nested sums of product of binomials and of harmonic 
numbers^ ff„, and, in rare cases, no "closed formula" is known for these coefficients. 

Let us show, in the case of An coefficients being nested sums of product of 
binomials, that the solution of the Calabi-Yau ODE, analytical at a; = 0, which is 
by construction a series with integer coefficients, is actually a diagonal of rational 
function, and furthermore, that this rational function can actually be easily built. 



5.1. Calculating the rational function for nested product of binomials 

For pedagogical reasons we will just consider, here, a very simple example§ of a series 
S{x), with integer coefficients, given by a sum of product of binomials 

3 



Si-) - EE 



ri=0 fe=0 

= ifeunG(-8,-2,l,l,l,l, 8x) = i7eiinG(-l/8, 1/4, 1, 1, 1, 1, -a;) 
= 1 +2a: +10a;^ +56x^ +346x^ 4- 2252 x-'^ -I- 15184 a;^ -f- 104960 

+ 739162 -h 5280932 -I- 38165260 a;^" (52) 

This is the generating function of sequence A in Zagier's tables of binomial coefficients 
sums (see p. 354 in [96]). 

The reader can easily get convinced that the calculations of this section can 
straightforwardly (sometimes tediously) be generalised to more complicated [97] nested 
sums of product of binomialsf. 

The diagonal of a rational function P/Q is written using Deligne's trick 

DiagQ = {^y-f (53) 

\Q/ \2tn/ Jc Q Zi 22 

where C avanishingf cycle [98], which is, with everyday words, the n- variables residue 
formula. Finding that a series is a diagonal of a rational function amounts to framing 
it into a residue form like (49). In order to achieve this, we write the binomial (JJ) as 
the residue 

1 f (1 + z)" dz 

^ The generating function of Harmonic numbers is H{x) = Hn ■ = — ln(l — x)/(l — x). 
§ See Proposition 7.3.2 in [90]. 

f Not necessarily corresponding to modular forms as can be seen on (81), (82). 
f Cycle evanescent in french. 
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and, thus, wc can rewrite S{x) as 
(2iTrf ■ S{x) = 

-til /Er-^^•((^ + ^^)(^ + ^^)(l + ^^)••^)"•^^^'^ 

t^aJ J J t^a (^1^223)'= V ; Z1Z2Z3 

/ \ (n+1) 

- 1 - (1/(Z1Z2^3)) . x« dzidZ2dZ3 

t'o 1 - (1/(Z1Z2Z3)) ^ ^ ^1^2^3 



/■ ^1 Z2 Z3 (, , , \« dzidz2dz: 



3 



^1 ^2 23 

1 / (I + zi) (1 + Z2) (1 + 2:3) ■ dzi dz3 

1 - Zl Z2 Z3 ^ Zi Z2 Z3 J Zi Z2 Z3 

ui \ '^^i dz2 dz3 

R(x; zi, Z2, Z3) ■ , (55) 

Zl Z2 Z3 



where R{x; zi, Z2, Z3) reads: 

Zl Z2 Z3 



(1 - a; • (1 + zi)(l + Z2)(l + Z3)) (zi Z2 Z3 - a; • (1 + zi)(l + Z2)(l + Z3)) ' 
From this last rcsuh one deduces immediately that (52) is actually the diagonal of: 

1 

(1 - Zo • (1 + Zl)(l + Z2)(l + Z3)) • (1 - Zo Zl Z2 Z3 (1 + Zl)(l + Z2)(l + Z3)) ' 

Note that, as a consequence of a combinatorial identity due to Strehl and 
Schmidt [99, 100, 101], S{x) can also be written as 

n=0 k=0 V / \ / „=0 fe=[,i/2]-l V / \ / 

Calculations similar to (55) on this other binomial representation (56), enable to 
express (52) as the diagonal of an alternative rational function: 

1 



(1 - Zo • (1 + Zl)(l + Z2)(l + Z3)2) • (1 - Zo Zl Z2 • (1 + Zl)(l + Z2)) 



(57) 



We thus see that we can actually get explicitly, from straightforward calculations, 
the rational function (56) for the Calabi-Yau-like ODEs (occurring from differential 
geometry or enumerative combinatorics) when series with nested sums of binomials 
take place, and, more generally, for enumerative combinatorics problems (related or 
not to Calabi-Yau manifolds) where series with nested sums of binomials take place. 

Remark: These straightforward effective calculations guarantee to obtain an 
explicit expression for the rational function (56), however the rational function is far 
from being unique, and worse, the number of variables the rational function depends 
on is far from being the smallest possible number. Finding the "minimal" rational 
function (whatever the meaning of "minimal" may be) is a very difficult problem. 



X 
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Recalling the well-known Apery series A{x), and its rewriting due to Strehl and 
Schmidt [99, 100, 101], 

n=0 k=0 \ / \ / n=0 k=0 j=0 \ / \ / 

1 + 5 a: +73x2 + 1445a;^ + 33001 x"* +•••, (58) 

A{x) is known to be the diagonal of the rational function in five variables I/R1/R2 
where i?i, R2 read [66]: 

i?l = 1 - Zq, i?2 = (1 - Zi){l - Z2){1 - Z3)(l - Z4) - ZQZ1Z2, 

as well as the diagonal of the rational function in five variables I/Q1/Q2 where Qi, Q2 
read [67, 49]: 

Ql = 1 - Z1Z2Z3Z4, Q2 = (1 - Z3)(l - 24) - zo • (1 + Zl)(l + Z2), 

and also the diagonal of the rational function in six variables I/P1/P2/P3 where 
Pi, P2, P3 read [66]: 

Pi = 1 - Zq Zl, P2 = 1 - Z2 - Z3 - ^0 Z2 ^3, P3 = 1 - Z4 " ^5 " ^1 2^4 ^5- 

A yet different diagonal representation for the Apery series, due to Delayguef, is 
provided by the diagonal of the rational function in eight variables: 

1 

(1 - Z4Z5Z6ZJ) • (1 - Zo • (1 + Z4)) • (1 - Zl • (1 + Z5)) • (1 - Z2 - Zg) • (1 - Z3 - Zj) ' 

Calculations similar to (55) on these new binomial expressions provides two new 
rational functions such that (58) can be written as the diagonal of one of these two 
rational functions. One is a rational function of five variables, of the form l/Q^i^ /Q'^^ 

Q[^^ = 1 -Z0Z1Z2Z3Z4 • (1 + Zi)(l + Z2)(l + Z3)(l + Z4), 

of = 1 -Zo • (1 + Zi)(l + Z2)(l + Z3)2(l + Z4)2, (59) 

and the other one, is a rational function of six variables, of the form l/Qi^V02^V'33^'' 

Qf'' = 1 - Zo Z3 Z4 Z5 • (1 + Zl) (1 + Z2f (1 + Z3) (1 + Z4) (1 + Z5), 
= 1 -Z0Z1Z2Z3Z4Z5 • (1 + Zi)(l + Z2), 
^ = 1 - zo ■ (1 + Zl) (1 + Z2f (1 + Z3) (1 + Z4) (1 + Z5). (60) 

We thus see that, when a given function is a diagonal of a rational function, the 
rational function is far from being unique, the "simplest" representation (minimal 
number of variables, lowest degree polynomials, ...) being hard to find. 

Similar computations show that the generating function of sequence E in Zagicr's 
list [96]: 

= 1 +4x +20x^ + 112a;3 +676x'' +4304.t^ + 28496 .t^ + 194240 x^ 
+ 1353508 X* +9593104x9 +68906320x1° + ••• (61) 

f Private communication. 
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is the diagonal of the rational function in four variables 

1 

(1 - 4Z0Z1Z2Z3 ■ (1 + Zi)) • (1 - Z^Z2Z3 ■ (1 + Z2)2(l + ^3)^(1 + ^l)^) ' 

while the generating function of Zagier's sequence B 

^ |e(-i)-3"-(;)c:)c: 

= 1 +3x +9x^ +21x^ +9x'^ - 297a;'5 - 2421x^ - 12933 

- 52407 X* - 145293 a;^ - 35091 a;i° + ••• (62) 
is the diagonal of the rational function in four variables 

^ (63) 



(1 - 3 Z0ZIZ2Z3 ■ {I + Zl)) ■ {1 +Z0^Z|Z3^ -(1 + 21)3(1+22)3(1 + 23)2) 



Such calculations can systematically be performed on any series defined by nested 
sums of product of binomials. We have performed such calculations on a large number 
of the series corresponding to the list of Almkvist et al [91], that are given by such 
nested sums of product of binomials. 

In fact, any s-nested sum of products of binomials raised to powers ii, . . . , it can 
be written as the diagonal of a rational function in £i + ■ ■ ■ + it + 1 variables, of 
the form ((1 — Qo){l — Qi) ■ ■ ■ (1 — Qs))~^, where the Qi's are products of powers 
of the variables Zi and of the linear forms 1 + 2^. 

For instance, when s ~ 1, it is easy to prove using the same technique that the 
power series 



oo n 



n=0 k=0 4=1 

is the diagonal of the rational function in p + 1 variables 

1 



(1 — • 2o 2l • • • Zp) • (1 — 0"+^ • UV Zq Zl ■ ■ ■ Zp) ' 

where u and b read: 



(65) 



—vr^ ' ^ = 11 — ■ (66) 

i=i 1=1 ^% 

The same machinery provides in some cases diagonal representations for algebraic 
power series (as diagonals of bivariate rational functions) that are much simpler than 
those produced by Furstenberg's result sketched in Section 3.5. For instance, using 
the fact that 

^2n - 2\ J_ I (1 +2)2"-^ dz 

n-l ) 211:' Jc z""^ ' -2 ' 

the algebraic function 



2 v-^/2n-2\/2\" 

' = 7T=. - "-i: „-i)(i) ■ 

^ n— 

is readily seen to be the diagonal of the rational function 



1 - 2o-(l + 2l)V4 
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which is much simpler than (25). 
Similarly, the power series 



E : - 

n=0 ^ ' 



(70) 



is seen to be the diagonal of the rational function 

1 



(71) 



1 -zo- (1 + zi)"' 



6. Comments and speculations 

6.1. Christol's theorem 

In [49] (page 61 Theorem 12, see also Proposition 7 in page 50 of [68] ) it is proved 
that any power series with an integral representation (as defined in (72) see below) 
and of maximal weight for the corresponding Picard-Fuchs linear differential equation 
(denoted by Ly below) is the diagonal of a rational function and, in particular, is 
globally bounded. 

The technical nature of the original papers is such that the result itself is 
difficult to find. This paragraph is devoted to explain, in down-to-earth terms, the 
somewhat esoteric expressions used in its wording, and to explain what it means on 
explicit examples. As the original proof is very obfuscated its principle is sketched in 
Appendix E. 

A function /, analytic near 0, is said to have an "integral representation" if it can 
be written in the following formf : 



where F is an algebraic function, hence living on some (projective) complex {n + 1)- 
fold V, and C is a "cycle" , namely an (oriented) compact (i.e. without boundary) real 
n-fold contained in V. In (72) x must be seen as a parameter. Then one integrates 
the n-differential F{x; . . .) dxi ■ ■ ■ dxn, that depends on x, on a "constant" cycle| C . 
If it exists, the integral representation is far from unique. Formula (49) shows that 
any diagonal of a rational function has an integral representation for which C is the 
so-called "vanishing [104] cycle". This is straightforwardly extended to diagonals of 
algebraic functions. 

In practical examples, to obtain a cycle, one often needs to complete the 
integration domain by means of symmetries in the usual way when dealing with the 
method of residues. That happens for instance for the hypergeometric function f{x) = 
B{a2, bi —02) • 2^i([ai,a2], [bi]; x) [B is the beta function, 3?(5i) > 3^(02) > 0), 
which has the following Euler integral representation: 



f Following [102] or, better, its version with parameter [103], one could define integral representations 
by integrating in (72) on domains C defined by polynomial (in xi, . . . , Xn) equalities or inequalities. 
Actually the case C = {xi S [0, 1] ; 1 < « < n} would be enough. To connect this definition with 
ours, one needs to use the Stokes theorem. However, the dimension n of the underlying complex 
manifold is basic for our next definitions and it is unfortunately not preserved by Stokes theorem. 
X One could, more generally, integrate an m-differential for m < n. But Lefschetz theorems assert 
that, up to taking hyperplane sections, one can reduce to the case n = m. 




(72) 




(73) 
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More generally the only hypergeometric functions pFg having an (Euler) integral 
representation are the n+iFn with rationalj: parameters [105, 106] a,;, bj: 

n+lF„i[ai, 02, ■■ ■ , On, On+l], [bl, 62, • • • , x) 

t Xl' ■■■ • (1 - Xi)"'-'''-^ ■■■ (1 - Xnf"-^"-' 



^-a ^, dxi dXn 

X (1 - a;2 • • • a;„ • x) °"+^ • • • . (74) 

Xi Xji 

Moreover, the integration cycle C of (72) can involve points at infinity (we are 
dealing with projective geometry). 

As C is a cycle, adding an exact differential to F dxi ■ ■ ■ dxn clearly do not change 
/. But a famous Grothendieck's theorem [107] asserts that n-differentials on V up to 
exact ones built up a finite C(a;)-space. Moreover^ differentiation under the integral 
endows this space with a connection, namely the Gauss-Manin one. In other words, / 
is solution of an ordinary linear differential equation (ODE) Ly, namely the Picard- 
Fuchs differential equation. This ODE depends only on V and do not involve any 
cycle but choosing a particular cycle C amount to choosing a particular solution / of 

Ly. 

The Picard-Fuchs ODE Ly is known to have only regular singularities with 
rational exponents. Now L being a linear ODE for which is a regular singularity, 
solutions near of L are endowed with a "monodromyl weight" : / is of weight W if 
L has + 1 solutions that are built in the following way: 

fix) , fix) - \ogix) +/i(x) , ••• , (75) 

log^(x) hg'^-'ix) 

where the fi are analytic near 0, and no solution involving fix) ■ log^^^(a;). For 
instance, a MUM (maximal unipotent monodromy) ODE L of order fi has a unique (up 
to a multiplicative constant) analytic solution near and this solution is of logarithm 
weight /i — 1. 

Geometric considerations imply that, for solutions of Ly, the maximum 
monodromy weight is n — 1. So we will say that / is of maximal weight for Ly if 
it is of weight n — 1 . 

The Picard-Fuchs linear ODE is difficult to determine and, moreover, depends on 
the particular integral representation. What is well defined is the minimal linear ODE 
L/ of which /, with an integral representation, is solution. Then its Picard Fuchs 
ODE Ly = ML/ is a left multiple of Lf and the monodromy weight of / for Ly 
is at least the monodromy weight of / for Lf (it is likely that the two monodromy 
weights are actually the same). The order of L/ is smaller, and often much smaller, 
than the order of Ly. 

(t The parameters have to be rational for the n-form in the integral representation (72) be 
algebraic. More deeply, when globally bounded, the hypergeometric function is a G-function and the 
corresponding linear differential operator is globally nilpotent hence it has only regular singularities 
with rational exponents. The parameters and bj are directly linked to exponents at and oo. 
^ See [108] specially § 1 and § 4 for instance but notice that the weight filtration used there and the 
monodromy weight we will use are quite distinct even if both constructions could look similar. 
J The monodromy operator T, "turning once counterclockwise around 0", acts on the space S of 
solutions of L by T{x) = xe-^"^ = x and T(log(x)) = log(a;) +2i-K. The solution / is of weight W 
if and only if it is in the image (T - 1)^(5) but not in (T - 1)^+^(5'). 
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There is absolutely no reason for the order oi Lf to be the number n of variables 
of the integral representation, but this effectively do happen in examples, notably for 
hypergeometric functions n+iFn and for certain Calabi-Yau linear ODE. Under these 
circumstances, if is MUM then / is of maximum weight for Ly and the theorem 
asserts that / is the diagonal of a rational function. Let us remark that in that case 
/ is the unique (up to a multiplicative constant) analytic solution near of Lf. 

Disappointingly, when this result can be applied to it becomes somewhat 

trivial. More precisely, the hypergeometric function is of maximal weight if and only 
if bj = 1 for all j (there is only n!'s in the denominator of coefficients). In that case 
it is obviously the Hadamard product of algebraic functions: 

„K_i([ai, a2, ■ ■ ■ , an], [1, 1, • • ■ 1], x) (76) 
= (1 -x)-"i ★ (1 -xy"^ ■■■ ★ (1 -x)-"". 

Therefore, we now have (at least) three sets of problems yielding diagonal of 
rational functions: the n-fold integrals of the form (49) with (50), the Picard-Fuchs 
linear ODEs with solution of maximal monodromy weight and, finally, the problems 
of enumerative combinatorics where nested sums of products of binomials take place. 
Diagonal of rational functions, thus, occur in a quite large set of problems of theoretical 
physics. 

6.2. Christol's conjecture 

The diagonal of a rational function is globally bounded (i.e. it has non zero radius of 
convergence and integer coefficients up to one rescaling) and D-finite (i.e. solution of 
a linear differential equation with polynomial coefficients)!. 

The reciprocal statement is the "Christol's conjecture" [49] saying that any D- 
finite, globally bounded series is necessarily the diagonal of a rational function. 

A fantastic Chudnovski theorem ([110] page 267) asserts that the minimal 
linear differential operator of a G function (and in particular of a D-finite globally 
bounded series) is a G-operator (i.e. at least conjecturally, a globally nilpotent 
operator) [28, 29, 30]. "Christol's conjecture" amounts to saying something more: 
if the solution of this globally nilpotent linear differential operator is, not only a G- 
series, but a globally bounded series, then it is the diagonal of a rational function. 

Conversely the solution, analytical at 0, of a globally nilpotent linear differential 
operator is necessarily a G- function [29, 30]. Moreover, a "classical" conjecture, with 
numerous avatars, claims that any G-function comes from geometry i.e. roughly 
speaking, it has an integral representation§. 

To test the validity of Christol's conjecture we look for counter-examples not 
contradicting classical conjectures. Then we search D-finite power series with integer 
coefficients which are not algebraic but have an integral representation and are not of 
maximal weight for the corresponding Picard-Fuchs linear ODE. 

As a first step let us limit ourself to hypergeometric functions „+ii^„. The 
monodromy weight W is exactly the number of 1 among the hi . 

When n+iFn is globally bounded and has no integer parameters hi [W — 0), its 
minimal ODE has ap-curvature zero for almost all primes p. However, a Grothendieck 

X The series expansion of the susceptibility of the isotropic 2-D Ising model can be recast into a series 
with integer coefficients (see [15, 27, 109, 19]), but it cannot be the diagonal of rational functions 
since the full susceptibility is not a D-finite function [109] . 
§ Bombieri-Dwork conjecture see for instance [30]. 
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3F2 



1 

-, 1 
3 



36 , 



x] = 1 +60 a; + 20475 + 9373650 



conjecture, proved for 3F2 in [108], and generalised to „+iF„ in [111], asserts that, 
under these circumstances, the hypergeometric function is algebraic. 

So we are looking for globally bounded hypergeometric functions satisfying 
1 < W < n — 1. In general such hypergeometric functions are G-series but are very 
far from being globally bounded. The hypergeometric world extends largely outside 
the world of diagonal of rational functions. 

Such an example in the first case n = 2, W = 1 was given in [49]: 

"1 4 5" 

9' 9' 9 

+ 4881796920 a;'' + 2734407111744 + 1605040007778900 x*' + ••• (77) 
The integer coefficients read with the rising factorial (or Pochhammer) symbol 
(l/9)„ ■ (4/9)„ ■ (5/9)„ 6„ ^ 

(l/3)„-(l)„-n! p(0)' ^ ' 

where: 

, , _ r(l/9 + n)r(4/9 + 7i)r(5/9 + n) 

^^"^ " r(i/3 + n)r(i + n)r(i + n) ■ ^^^> 

Note that, at first sight, it is far from cZear§ on (79), or on the simple recursion on 
the p{n) coefficients (with the initial value p{0) = 1) 

Pjn + l) ^ (1 + 9n) (4 + 9n) (5 + 9n) 
p{n) (i + 3n)(i + n)2 ' y'''') 

to see that the p{n)'s are actually integers. A sketch of the (quite arithmetic) proof 
that the p{nys are actually integers, is given in Appendix G. 

Because of the 1/3 in the right (lower) parameters of (77), the hypergeometric 
function (77) is not an obvious Hadamard product of algebraic functions (and thus 
a diagonal of a rational function), and one can sec that it is not an algebraic 
hypergeometric function either by calculating its p-curvature and finding that it is 
not zero [93], or using [112]. Proving that an algebraic function is the diagonal of a 
rational function and proving that a solution of maximal weight for a Picard-Fuchs 
equation is the diagonal of a rational function use two entirely distinct ways. The 
hope is to combine both techniques to conclude in the intermediate situation. 

This example remained for twenty years, the only "blind spot" on Christol's 
conjecture. We have recently found many other 3F2 examples^, such that their series 
expansion have integer coefficients but are not obviously diagonal of rational functions. 
These new hypergeometric examples are displayed in Appendix F. Unfortunately these 
hypergeometric examples are on the same "frustrating footing" as Christol's example 
(77): we are not able to show that one of them is actually a diagonal of a rational 
function, or, conversely, to show that one of them cannot be the diagonal of a rational 
function. 



7. Integrality versus modularity: learning by examples 

A large number of examples of integrality of series-solutions comes from modular 
forms. Let us just display two such modular forms associated with HeunG functions 
of the form HeunG{a, q, 1, 1, 1, 1; x). 

§ In contrast with cases where binomial (and thus integers) expressions take place. 

X 2F1 cases are straightforward, and cannot provide counterexamples to Christol's conjecture. 
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7.1. First modular form example 

One can, for instance, rewrite the example (52) of subsection (5.1), namely 
FeMnG(-l/8, 1/4,1,1,1,1 hypergeometric functionf with two rational 

pullbacks: 

oo n / \ 3 

iJe7.nG(-l/8,l/4,l,l,l,l,-x) = ^ E ( J (^1) 

71=0 k=0 ^ ^ 



(( 



1 +4x) • (1 + 228a: + 48 a:^ + 6Ax^)^ 

, /, 1 5,.^. 1728- (1 - 8a;)6 • (1 +x)3 -a; 

^l[l2' 12]'[^' 



X 2F1 



(^(1 -2a;) • (1 -6a; + 228a'2 -8a;^)) 



(1 + 228a; +48a;2 + 64a;3)3 • (1 +4x)3 

-1/4 



,1 5 1728.(1 - 8a;)3.(l +a;)«.a;2 

X 2^1 [-7:, -r7:\, [IJ 



'12' 12" ' " (1 - 2a;)3 • (1 - 6x + 228a;2 -8a;3)3. 

The relation between the two pullbacks, that are related by the "Atkin^" involution§ 
X O —1/8/x, gives the modular curve: 

1953125 - 187500 ■ {y + z) + 375 y z • (16 z2 - 4027 z y + 16 j/^) 

-U- {z + y) ■ {y^ + z^ zy) + 110592- zy = 0. (82) 

Series (52) is solution of the (exactly) self adjoint linear differential operator VL 
where {9 = x ■ D^): 

x-n = 02 - a;- (76*2 + 761 + 2) - 8 a;^ • (6* + 1)^. (83) 

The relevance of the "Atkin" involution x o — l/8/a;is also clear on the operator: 
note that operator is invariant by changing 6* f^- — 1 — 6*, and x O —1/8/x. Also 
note that changing £7 by a pullback x <-> — 1/8/a;, amounts to changing into 

X ■ n ■ x~^. 

7.1.1. Modular invariance 

Do note that these pullbacks are respectively of the form (o denotes the compo- 
sition of functions): 

1728a; (l~8xf , ~ 1728a; (l-8a;)3 

M2 = 7 -TTT o At and: M2 - 7 ^7777^ o ^— 84 

(a; + 16)3 a;(l + .T)3 (.t + 256)^ x {1 + x)^ 

f The relation between modular forms and hypergeometric functions can be simply seen in the 
identity 2-Fi([l/12, 5/12], [1], 1728/ jir))* = E4,{t), where £4 is an Eiscnstein series. 
^ In previous papers [113, 32], with some abuse of language, we called such an involution an Atkin- 
Lehner involution. In fact this terminology is commonly used in the mathematical community for 
an involution t — > —N/t, on r, the ratio of periods, and not for our x-involution. However, lu/ien 
the modular curve is of genus zero one has a parametrisation in term of the variable xi^(t) = 
(r;(T)/77{ArT))24/(^-i) (see eq. (27) in [114]), which actually transforms as an involution (xjv — > 
A/xff) under the Atkin-Lehner involution. This is why we switch to the wording "Atkin" involution. 
§ The relevance of the "Atkin" involution x •<-*■ —1/8/x is also clear on the operator: note that 
operator Q is invariant by changing 9 •o- —1 —6, and x -o- —l/8/x. Also note that changing H by 
a pullback x <-> —1/8/x, amounts to changing Q, into x ■ ■ x^^ . 
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where one recognises the two Hauptmoduls of the modular curve corresponding to 
T 2 T. Introducing the Dcdekind-hke parametrisation: 



3t- (1 -t^) 

1-^6 ■ 



a;(g) = q. where: ^ a(n) • 

n—1 n—1 

one can rewrite the Hauptmoduls as: 

~ / (l-8x(g))3 X ^ / (l-8x(g^))^ X 

Ha;(<7) • (1 + x(g))3y ' Vx((;2) . (1 + a;((j2))3y ' ^"^^ 

7.1.2. Other representations 

In fact, using Kummer's relation and other relations on 2-^1 's: 

.1 1 „ /4 2 



1/4 / I 5 64x3. (l^x) 



.9 - 82;/ V^12' 12^' ' (9 - 8x)3 

HeunG{— 1/8, 1/4, 1, 1, 1, 1, —x) can be written in many different ways as hypergeo- 
metric 2F1 with two puUbacks: 

1 / 1 2 27 a; 

H™nG(-i/8.i/4,i,i,i.i.-i) = rn^-^'^'Os'si-w-oTi^F 

= (1 + -8.,-"- 1], 111 

1 ^ /.I 1. 108a;(l + a;)(l -8a;)2 



1 + 4a; ' 'V^6' 3^' ' " (1 + 4a;)6 

The relation between the two pullbacks in (87), namely u = 27x/{l + Ax)^, 
and w = 27xV(l- 2x)3),isthe genus-zero modular curve: 

8u^v^ -Uu'^v^ ■ (u + v) +3uv {2u^ + 2v'^ + 13uv) (88) 
- {u + v) ■ {v'^ + 29 uv + u'^) +27uv = 0. 

Let us consider the modular curve (88) for u ^ x, seeing u as a algebraic function of 
X. One of the three root-solutions expands as: 

~ / X 1 2 10 . 256 4 18928 . 154000 . 

vo(x) = —x^ H x^ -\ x^ ~\ a;^ H +■■■, (89) 

' 27 243 6561 531441 4782969 ' v ^ 

the two other ones being Puiseux series (here t denotes a;^/2); 

r- o 119 ^ , 1904 4 50701 ^ r 

6 27 648 

the third root, V2ix) corresponding to change t into —t. These series expansions 
compose nicely: vo{vi{x)) = vo{v2{x)) = x. 
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Remark: Ramanujan's cubic transformation. For the hypcrgcometric 
function 2^i([l/3, 2/3], [1], x), the existence of two pullbacks is also reminiscent of 
relation (28) in [33], or to the known Ramanujan's cubic transformation formula [115]: 

(l + M-,F,([i|l.|ll,27.») ^ ||, 11], 1 - 

1 +6x +6a;3 +36x^ +90x^ + 540 + 1680 + ••• (91) 

Noting that 2i^i([l/3, 2/3], [1], x) and (^[1/3, 2/3], [1], 1 -a;) are solutions of the 

same second order linear ODE, one can change one of the two pullbacks. Pi ~ 27 x^ 
and P2 = 1 - (1 - 3a;)3/(l + 6x)^ in (91), into 1 - Pi, i ^ I, 2. The relation 
between the two pullbacks u = 1 — Pi and v ~ P2, is a simple (genus zero, (w, w)- 
symmetric) modular curve: 

512 ■ ■ ~ 1728 ■ ■ ■ (u + v) +216- uv {17- uv + 9- {u^ + v^)) 

- 243- (w + u)- (3- (u^ + w^) +S-UV) + 729 ■ (u^ + uv + v^) = 0. (92) 



7.1.3. Schwarzian condition 



A necessary condition for two different rational (resp. algebraic) pullbacks to exist 
for a hypergeometric function like 2Pi([l/3, 2/3], [1], x), i.e. a necessary condition 
for a relation 



([i, ^], [l],pi(x)) = ri,2(x). 2Pi([^, [l],P2(x)), (93) 

for some algebraic functions ^1^2(2;), is the (symmetric) condition 
1 8pi{x)^ -Spijx) +9 /dpi(x)\2 

iPl[X),Xi + p^(^)2.(p^(2,)_l)2 ■{ dx J 

where {pi{x), x} denotes the Schwarzian derivative. This condition is invariant by 
the simple transformations pi — > 1— pi, i = 1,2. One immediately verifies that 
the Schwarzian condition (94) is actually verified for the pullbacks (pi(x), P2ix)) 
occurring in (91), namely (27 x^, 1 — (1 — 3a;)'^/(l + 6x)^), or in (87), namely 
(271/(1+ 4a;)3, 27x7(1- 2x)^). 

Let us consider the modular curve (92) for u = x, seeing d as a algebraic function 
of X. The three root-solutions expand as: 

1 . 5 4 56 . 1691 5390 7 

voix) = 1 X X X X^ X + 

°^ ' 729 2187 19683 531441 1594323 

3 , f^,V3-- 




70 . /- 532\ 5 /1171 2297 . /-\ g 
243 ^^+ 729j ■ +(,1458 + 13122 ^^; ^''^ 

the third series expansion V2 (x) having its coefficients complex conjugate of the ones 
in the series expansion vi . These series expansions compose nicely: 

vi{v2{x)) = V2{vi{x)) = X, vo{vi{x)) = Vo{v2{x)) = vo{x), (96) 

Vi{vi(x)) = V2ix), V2{V2{X)) = Vi{x), Vi(vi(vi(x))) = W2 (^^2 (^2 (a^))) = X. 
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One verifies on these three series expansions (95) the condition corresponding to 
the pi{x) = X subcase in (94): 

1 8v{x)'^ -8v{x) +9 /dv{x)\^ _ 1 8x2-8x+9 

{V{X), 4 + • 1,(^)2 . _ 1)2 • \r^) ~ ■ x2.(x-l)2 • 

One also verifies that the three series expansions (89) and (90) also satisfy, as they 
should, the previous Schwarzian condition (97). One also verifies, as it should, that 
the composition of all these series satisfy (97). For instance 

~ ,~ , 1 4 20 . 274 g 86636 , 

^ ^ " 19683 177147 1594323 387420489 
1 10 1 S7 

voivolx)) = 1 x^ x'' x^ + (98) 

^ ^ " 14348907 43046721 387420489 ^ \ > 

satisfies (97). 

In other words, the Schwarzian condition (94) is another way to encode the 
modular curves (92) or (88), and in fact, an infinite number of modular curves. The 
emergence of Schwarzian derivatives should not be seen as a surprise. A relation like 
2Fi([l/3, 2/3], [1], v{x)) = Ri,2{x) ■ 2i^i([l/3,2/3], [1], x) is obviously stable by the 
compositionf of the puUback functions v(x). 

7.2. Second modular form example 

The integrality of series-solutions can be quite non-trivial like the solution of the 
Apery-like operator 

n = X ■ {I - llx -x^)-Dl +{l -22x -Zx^)- D:c - (x + 3), (99) 
or: x-VL ^ -x- (116*2 + lie + 3) -x^-{e+lf, 

which can be written as a HeunG function. Introducing a = 11/2 — 5 • 5^/^/2^ this 
(at first sight involved) HeunG function reads: 

HeunCi- — ■ - ^ + 1^ . ^^/\ 1, 1, 1, 1, « • x 

\ 2 2 2 2 ''''' 

. ^„e.„G(i- 11^ 5./^i-1.5./M. 1,1,1, 

= i:i:(:)7"r 

1 4- 3 • a; 4- 19 • -)- 147 • -f 1251 • -f 11253 -x^ -f 104959 • x*^ 
+ 1004307 • + 9793891 • x^ + 96918753 • x^ + 970336269 • x^^ 

9807518757 • -f 99912156111 • + 1024622952993 • 
+ 10567623342519 • 109527728400147 • x^^ + 1140076177397091 • x^^ 

+ 11911997404064793 • 124879633548031009 • x'^^ 

-h 1313106114867738897 -cc^^ + ■■■ (100) 

f The Schwarzian derivative S{f) = {/(x), x} is the well-suited derivative to take into account 
when composition of functions occurs. This can, for instance, be seen on the chain rule: {S{f o g))(z) 
= S(f){g(z))-g'(z)^ +S(g). 
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but actually corresponds to a modular form, which can be written in two different 
ways using two pullbacks: 

{x^ +12x^ + Ux^ - 12a; + 1)-^/^ 



X 2^1 



5 ^ 1728 -a;^- (1 - 11a; - x^) 

(a;4 + 12a;3 + 14a;2 - 12a; + 1)^ 



(1 + 228x + 494 a;2 - 228 a;^ + x^)"^^'^ (101) 



.12' 12^' ' " (1 + 228a; + 494a;2 - 228x3 + a;4)3 
Do note that these two pullbacks are respectively of the form: 
1728 -a;^ • (1 - 11a; ~ x^) 



(a;4 + 12a;3 + 14a;2 - 12 a; + 1)^ 



(102) 



1728 -x 1 - 11a; - a; 

o 



2 



and 



(a;2 + lOx +5)3 

1728-x - (1 - 11a; - x^)^ 
(1 + 228a; + 494 a;2 - 228 a;3 + a;4)3 

1728 -x^ 1 - 11a; 

o 



(103) 



(a;2 + 250.x +3125)3 

where one recognises [116] the two Hauptmoduls of the modular curve corresponding 
to r — >■ 5 • r: 

•^=<^) = WT^TTW- 

, ^ 1728 •a;5 /53 

= (x2 + 250 X +3125)3 ^ \- 

Also note that the fact that the relation between the two Hauptmoduls (104) 
corresponds to r — >■ 5 • r can be seen straightforwardly if one introduces the quite 
non-trivial Dcdekind-like parametrisation 

^(9) - 9 - ft (1 - '?")"^"'' ^here: £ a{n) ■ t^ = ^ ' ' ' 'J/' ~ 

n— 1 n— 1 

yielding: 

^.(i^li^ii^^^) ^ M.i'-^^^^^^y (105) 

More modular form examples of series with integer coefficients are given in 
Appendix H. The modular form examples displayed in Appendix A, corresponded to 
lattice Green functions [52]. Therefore, they have n-fold integral representations, and, 
after section (4.4) , can be seen to be diagonals of rational functions. In contrast the 
modular form examples displayed in Appendix H correspond to differential geometry 
examples discovered by Golyshev and Stienstra [117], where no n-fold integral 
representation is available at first sight. 
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7.3. Remark 



At first sight one might think, and this is almost suggested in the Uterature, that the 
integraUty of the series-solutions of the globally nilpotent operators, corresponds to 
some deep arithmetic property and, in the same time, that these integer coefficients 
have a deep "physical" meaning (instantons, ...). There is often some confusion in the 
literature between the concept of integrality of series (globally bounded series) and the 
concept of modularity [34] which suggests a connection with selected algebraic varieties 
(modular forms and elliptic functions [118], mirror maps, Calabi-Yau manifolds). 

Note that, along this "selected algebraic varieties" line, it is also worth recalling 
Krammer-Dcitwciler's counter-example [119, 120] to Dwork's conjecture of a globally 
nilpotent operator that cannot be reduced to an operator having hypergeometric 
solutions up to a puUback, which corresponds to the two following HeunG functions: 

HeunG{81, 1/2, 1/6, 1/3, 1/2, 1/2, Six), (106) 
a;i/2 . HeunG{81, 21, 2/3, 5/6, 3/2, 1/2, 81a;). 

These two HeunG arc solutions of a globally nilpotent linear differential operatorf 
but are not globally bounded: they are G-series that cannot be reduced to series with 
integer coefficients. However, this special example corresponds to a special arithmetic 
situation: we are dealing, here, with very special HeunG functions associated 
with periods of a family of abelian surfaces over a Shimura curve [120], and very 
special ODEs corresponding to uniformising linear differential equations of arithmetic 
Fuchsian lattices. 

The interest of demonstration of sections (4.2), (4.4) is to show that this 
integrality can be, in fact, the straight consequence of a quite simple form of the 
n-fold integral namely (49) with the analyticity condition (50), that can be expected 
to be seen in a quite general framework (for instance in enumerative combinatorics), 
far from the theory of elliptic curves, or the theory of Calabi-Yau manifolds and other 
Frobenius manifolds [121]. Let us now try to "disentangle" the concept of integrality 
of series (globally bounded series) and the concept of modularity. 



8. Integrality versus modularity 



8.1. Dijfeomorphisms of unity pullbacks 

Let us consider a first simple example of a hypergeometric function which is 
solution of a Calabi-Yau ODE, and which occurred, at least two times in the 
study of the Ising susceptibility n-fold integrals [32, 33] x'"^ a-nd Xd^\ namely 
4F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1], 256a;), where we perform a (diffeomorphism of 
unity) puUback: 



4i^3([i [1,1,1], — % ^ ) = 1 +16-a; (107) 

\2 2 2 2 1 + ci X + C2 x"^ + ■ ■ ■ / 

+ (1296 - 16 ci) • a;^ + (160000 + 16 - 16 C2 - 2592 ci) ■ x^ + ■■■ 



If the pull-back in (107) is such that the coefficients c„, at its denominator, are 
integers, one finds that the series expansion is actually a series with integer coefficients, 
for every such puUback (i.e. for every integer coefficients Cn). Furthermore, a 



f Generically HeunG functions with rational parameters are not globally nilpotent, which indicates 
that they do not have an integral representation as an integral of some algebraic integrand. 
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straightforward calculation of the corresponding nome q(x) and its compositional 
inverse (mirror map) x{q), also yields series with integer coefficients: 

q{x) = X +(64-ci)-a;2 +(c? + 7072 -C2 -128ci)-a;3 +•••, (108) 
x{q) = q + (ci - 64) • g2 + [c\ + 1120 + C2 - 128 Ci) • + . . . , (109) 

when its Yukawa coupling [32], seen as a function of the nome q, K{q) is also a series 
with integer coefficients and is independent of the pullback: 

K{q) = 1 + 32 • g + 4896 • g2 + 702464 • g3 (no) 

This independence of the Yukawa coupling with regards to pullbacks, is a known 
property, and has been proven in [88] , for any pullbacks of the difFcomorphism of unity 
form p{x) = a; + • • • 

Remark: Seeking for Calabi-Yau ODEs, Almkvist et al. have obtained [91] a 
quite large list of fourth order ODEs, which are MUM by definition and have, by 
construction, the integrality for the solution-seriesi analytic at x = 0. Looking 
at the Yukawa coupling of these ODEs is a way to define equivalence classes up to 
pullbacks of ODEs sharing the same Yukawa coupling. This "wraps in the same bag" 
all the linear ODEs that are the same up to pullbacks. Let us recall how difficult 
it is to see if a given Calabi-Yau ODE has, up to operator equivalence, and up to 
pullback, a hypergeometric function solution [32, 33], because finding the pullback 
is extremely difficult [32, 33]. We may have, for the Ising model, some „+iF„ 
hypergeometric function prejudice [32, 33]: it is, then, important to have an invariant 
that is independent of this pullback (we cannot find most of the time). 

Remark: The Yukawa coupling is not preserved by the operator equivalence. 
Two linear differential operators, that are homomorphic, do not necessarily have the 
same Yukawa coupling (see Appendix J). 



8.2. Yukawa couplings in terms of determinants 

Another way to understand this fundamental pullback invariance, amounts to rewriting 
the Yukawa coupling [122, 88], not from the definition usually given in the literature 
(second derivative with respect to the ratio of periods), but in terms of determinants 
of solutions (Wronskians, ...) that naturally present nice covariance properties with 
respect to pullback transformations (see Appendix J). 

We have the alternative definition for the Yukawa coupling given in Appendix J: 



where the determinantal variables W^'s are the determinants built from the four 
solutions of the MUM differential operator. This alternative definition, in terms 
of these W^'s, enables to understand the remarkable invariance of the Yukawa 
coupling by pullback transformations [33]. These determinantal variables Wm quite 
naturally, and canonically, yield to introduce another "Yukawa coupling" (which, in 
fact, corresponds to the Yukawa coupling of the adjoint operator (see J. 12)). This 
"adjoint Yukawa coupling" is also invariant by pullbacks. It has, for the previous 
example, the following series expansion with integer coefficients: 

K*{q) = 1 + 32 • q + 4896 • g2 + 702464 • g3 + (112) 

X Hence, these operators are G-operators. One can also, calculating their p-curvatures, check, 
directly, that the Calabi-Yau operators in Almkvist et al. [91] tables are actually globally 
nilpotent [28], thus yielding automatically the rationality of the exponents for all the singularities 
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which actually identifies with (110). The equality of the Yukawa coupling 
for this order-four operator, and for its (formal) adjoint operator, is a 
straightforward consequence of the fact that the order-four operator annihilating 
4-p3^[^j 5' 5' 5]' I-*^' -'^1' 256 is exactly self-adjoint, and, more generally, of the 
fact that the order- four operator, annihilating (107), is conjugated to its adjoint by a 
simple function (which is nothing but the denominator of the pullback) . 

8.3. Modularity 

This example, with its corresponding relations (107), (108), (110), (112) may suggest 
a quite wrong prejudice that the integrality of the solution of an order-four linear 
differential operator automatically yields to the integrality of the nome, mirror map 
and Yukawa coupling, that we will call, for short, modularity" . This is far from being 
the case, as can be seen, for instance, in the following interesting example, where the 
nome and Yukawa coupling K{q) do not correspond to globally bounded series, when 
the solution of the order-four operator as well as the Yukawa coupling seen as a 
function of x, Kix), are, actually, both series with integer coefficients. 

Let us consider the following 4^3 hypergeometric function which is clearly a 
Hadamard product of algebraic functions and, thus, the diagonal of a rational function: 

4^^3([^, i \, J], [1, 1, l],x) 

= (1 - .t)-1/3 * (1 - x)-l/2 ^ (1 _ ^)-l/4 ^ _ ^)-3/4 
= Diag((l - z,)-^/^ (1 - z,)-^/^ (1 - Z3)-l/4 (1 _ ,^)-3/4^^ 
It is therefore globally bounded; 
4F3([i, i, ^, [1, 1, 1], 2304a;) = 1 + 72x +45360a;=^ -f 46569600 

59594535000 H- 86482063571904 -f 136141986298526208 a;^ 
+ 226888189910421811200 + 394399917777684601926000 
4-708188604075430924446000000x9 +■■■ (113) 

Its Yukawa coupling, seen as a function of a;, is actually a series with integer 
coefficient in x: 

K{x) = 1 +480 a; + 872496 -f- 1728211968 a;^ 3566216754432 a;^ 

+ 7536580798814208 x^ + 16177041308360579328 x^ (114) 
+ 35105183794659521064960 a;^ + 76799014669577085362391024 a;^ 
169059790576811511759706311168 a:^ + ••• 
However, do note that the scries, in term of the nome, is not globally bounded: 

K{q) = 1 -I- 480 g -I- 653616 + 942915456 q^ + 1408019875200 q^ 

+ 2146833138536640 + ••• (115) 

12 

-1- 571436303929319146711343817202689132288 ^ + ••• 

In fact, the nome q{x), and the mirror map x{q), are also not globally bounded. 
Note that in this example, the non integrality appears at order twelve (for x{q), q{x) 
and K{q)). If the prime 11 in the denominator in (115) was the only one, one could 
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recast the series into a series with integer coefficients introducing another rescaling 
2304 a; — )• 11 x 2304 x. But, in fact, we do see the appearance of an infinite number 
of other primes at higher orders denominators in x{q), q{x) and K{q). 

8.4- Hadamard products of w„ 's 

Let us consider the two order-two operators 

which arc associated with two modular forms corresponding, on their associated nomes 
g, to the transformations q ^ and q ^ q^ respectively (multiplication of r, the 
ratio of their periods by 2 and 3), as can be seen on their respective solutions: 

2F^f[i, il ri], _64x) = (1 + 256x)-^/^- 2Fi([—, —1, [1], , 
2 lyi^, ^l, I I, j \ -r , ^ ^42' 12^' ^ ^' (1 + 256x)'^ 

X 1 M /r 1 5 , , , 1728 \ 

- ' -^^([i2'i2]'W'(rTw) 

= 1 -4x +100^2 -3600a;^ + 152100a;'' - 7033104 a;^ + 344622096 a;*^ 
- 17582760000 a;^ + 924193822500 a;^ - 49701090010000 a;^ + • • • 

,-1/12 /I 5 , _ 1728x 



((l+27a;)(l+ 243xf)" . ^F,([-, [I], 



(1 + 243x)3 (l +27x). 

/ \ -1/12 / 1 5 1728t-^ \ 

((1 +27a;)(l +3a;)3) . ([—,—], [1], ^ ^) (119) 

M -r ; y 2 1VL12' 12^' ^ J' (1 +3a:)3(l +27a;)/ ^ ' 

2-F'i([^, [1], -27x) = 1 -3a; +36 a;2 - 588 a;3 + 11025 a:^ - 223587 a:^ 

+ 4769856 a;^ - 105423552 a;^ +2391796836 a;^ - 55365667500 a;^ + ••• 



The relation between the two Hauptmodul puUbacks in (118) 
1728 a; 1728 a;^ 



(120) 



(1 + 256x)3' (1 + 16a;)3' 

corresponds to the (genus-zero) fundamental modular curve: 

1953125 u^v'^ - 187500 u^v"^ ■ {v + u) + 375 uv {iQu^ + 16 - 4027 uii) 

- 64 (w +u) • (u^ + 1487iiu + M^) + 110592 = 0, (121) 
The relation between the two Hauptmodul puUbacks in (119) 

1728 a; _ 1728 a;^ 

" ~ (1 +243a;)3(l +27x)' " " (1 +3a;)3(l +27x)' ^ ' 

corresponds to the (genus-zero) modular curve: 

262144000000000 u^v"^ ■ {u + v) + 4096000000 u^w^ • (27 + 27^^ - 45946™) 

+ 15552000 uv {u +v)- (u^ + 241433 uv + u^) (123) 
+ 729 [u^ + t;**) - 779997924 {u\j + uv^) + 1886592284694 u^v^ 
+ 281 1677184 w (u +t;) - 2176782336 = 0. 
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Similarly, one can consider the order-two operators a;„ associated with other 
modular forms corresponding to t — > n ■ t. The a;„'s can be simply deduced 
from Maier [123], for modular forms corresponding to genus-zero curves i.e. for 
n = 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 16, 18, 25. After a simple rescaling, one gets 
series with integer coefficients. For instance considering the linear differential operator 
£7 (annihilating the modular form /17) in Table 13 of [116] 

2 7x'^+65x+U7 Ax + 21 

' ~ ^ + 3(a;2 + 13:r + 49)x " ^ 9 (x^ + 13 x + 49) • x ' ^ ^ 



one has the modular form solution 



2.2/3 

1 11 , 380 , 3887 

1 ~ ^ X 



21 3087 1361367 200120949 



where 

(.t2 + 13x + 49) (a;2 + 5a; + l)^ /49\ 

h{x) = ^ '-^ jj{x) = j7[ — ] 

X \ X / 



- (x2 + 13. + 49Hx2 + 5. + l)3 ' - ^Kf)' (^26) 

The series (125) is globally bounded. Rescaling the x into 3^ 7^ - x = 441- x, Maier's 
linear differential operator (124) becomes 

^ r.2 , 1 + 195X+ 9261^2 21- (1 + 84x) 

^ " (1 +117x-f 3969a;2). a; " (1 + 117 a; + 3969 rr^) . 3. ' ^ '> 

which has the following series with integer coefficients: 

i,,(441 f,([l,|l, [11;-^) (128) 

= 1 -21.T -t-693a;^ - 23940 a;^ + 734643 a;'' - 13697019 x'^ 

- 494620749 a;^ + 83079255420 .t^ - 6814815765975 + ■■■ 

The two operators L02 and W3 have a "modularity" property: their series 
expansions analytic at a; = 0, (118) and (119), as well as the corresponding 
nomes, mirror maps are series with integer coefficients. The Hadamard product is 
a quite natural transformation to introduce because it preserves the global nilpotence 
of the operators, it preserves the integrality of series-solutions, and it is a natural 
transformation to introduce when seeking for diagonals of rational functions^. Let 
us perform the Hadamard product of these two operators. With some abuse of 
language [33], the Hadamard product of the two order-two operators (116) and (117) 

_ 4 , . (2064X- 1) , (19020 X- 7) , 

- +^1728x-l).x-^- +(1728x-l).x2-^- 

(4788 x- 1) 12 , , 

' -D, + — — ^, (129) 



(1728x-l)-x3 (1728x- 1) • x3' 

^ And, consequently, has been heavily used to build Calabi-Yau-likc ODEs (sec Almkvist et al. [88]). 
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is defined as the (minimal order) linear differential operator having, as a solution, the 
Hadamard product of the solution-series (118) and (119), which is, by construction, a 
series with integer coefficients. This series is, of course, nothing but the expansion of 
the hypergeometric function: 

4i^3([J, J, J, J], [1, 1, l],1728x) (130) 

= 2F^{[^, i], [1], -64x) * 2i^i([i, i], [1], -27 x). 

The Hadamard product of the order- two operator (116) with itself (Hadamard 
square) 

^ Dt +2p^^^^.Dl ^ (42496.- 7) 

(4096x-l)-x (4096a;- 1) • x2 ^ 

(9984^, 16 

(4096 a;- 1) • x3 (4096 a; - 1) • a;^ ^ ' 

is defined as the (minimal order) linear differential operator having the series-solution 
1 + 16 a; + 10000 a;2 -I- 12960000 a;3 +23134410000 x'' 49464551874816 a;^ +■■■, 

which is the Hadamard product of the solution-series (118) with itself. This series is, 
of course, nothing but the expansion of: 

\, J, [1, 1, 1], 4096a;) (132) 

= 2Fi([i, i], [1], -64a;) * 2Fi([i, i], [1], -64a;). 

This operator i/2,2 is a MUM operator. We can, therefore, define, without any 
ambiguity, the nomc (and mirror map) and Yukawa coupling of this order-four 
operator [33]. One finds out that the nomef, and the mirror map (and the Yukawa 
coupling as a function of the a; variable), are not globally bounded: they cannot be 
reduced, by one rescaling, to series with integer coefficients. 

Similarly, one can also introduce the Hadamard square of (117) 

_ 4 , . (891 a; -1)^ (1215a;- 1) 

^^"^ - +^729a;-l).a;'^^ +^(729^^1^^-^^ 

I (2295 a; -1) , 9 .^33^ 

+ (729x-l)-x3 + (729X- 1)- x3' ^ ' 
which has the hypergeometric solution: 

aF3{[^, ^' ^' 1' l]'729a;) (134) 

= 2Fi([i, i], [1], -27 x) * 2Fi([i, 1], [1], -27 x). 

Let us remark that the three linear differential operators (129), (131) and (133), 
are MUM and of order four. However, they are not of the Calabi-Yau type. 

f The nome of the Hadamard product of two operators has no simple relation with the nomc of these 
two linear differential operators. 
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8.5. Hadamard products versus Calabi-Yau ODEs 

This is not the case for other values of n and m. For instance one can introducel 
i?4,4 = UI4 -k L04, the Hadamard square of W4, which is an irreducible order-four 
linear differential operator, and has the hypergeometric solution already encountered 
for some 71-fold integrals of the decomposition of the full magnetic susceptibility of 
the Ising model [32, 33]: 

4ml, I' i ^' l]'256a;) (135) 

The associated operator having (135) as a solution, obeys the "Calabi-Yau 
condition" that its exterior square is of order five. 

Let us give in a table the orders (which go from 4 to 20) of the various 
Hm,n = Hn^m Hadamard products of the order-two operators associated with the 
(genus-zero) modular forms operators cOn and uj„i. 
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where the star * denotes Calabi-Yau ODEsjl. 



The following operators are of order four: H2.21 H2.3, -^2.6: -f^2,8i ^2,9i 

H3.3, H3 4, H3 Q, i/a.s, ^^3,9, ■•■ Their exterior squares, which are of order six, do not 
have rational solutions^. 

The following operators are of order six: i?2,5, -^2,7, ^^3,5, H^j, i?4,5, i?4,7, 
-f^s.S: ^5,61 -^5,81 ^^5.9i Hqj, Hj t, Ht sjHt q, ... Their exterior square, which are 
of order fifteen, do not have rational solutions (and cannot be homomorphic to higher 
order Calabi-Yau linear ODEs). 

Remarkably the following ten order- four operators -ff4.4, -ff4,6: -^4,8; -f^4,9, ^^6,6, 
-^6,81 -^^6.9, ^^8,8> -^8,9: -^9,9 (with a star in the previous table) are all MUM, and 

X To get the Hadamard product of two linear differential operators use, for instance. Maple's 
command gfun[hadamardproduct]. 

(j Recall that Calabi-Yau ODEs are defined by a list of constraints [88], the most important ones 
being, besides being MUM, that their exterior square are of order five. There are more exotic 
conditions like the cyclotomic condition on the monodromy at 00, see Proposition 3 in [91]. 
^ They cannot be homomorphic to Calabi-Yau ODEs. 
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are such that their exterior squares are of order five^: they are Calabi-Yau ODEs. 
Actually the nome, mirror map and Yukawa coupling series are series with integer 
coefficients for all these order four Calabi- Yau operators. The Yukawa coupling series 
of these Calabi-Yau operators are respectively, for _ff4.4 

K{q) = K*{q) = 1 + 32 • g + 4896 • + 702464 • + 102820640-/ 

+ 15296748032 • q^ + 2302235670528 • q^ + 349438855544832 • q'' 

+ 53378019187206944- + 8194222260681725696- (j^ + (136) 

which is Number?) in Almkvist et al. large tables of Calabi-Yau ODEs [91], and is 
the well-known one forft 4^3([l/2, 1/2, 1/2, 1/2], [1,1,1], 256 a;), for Hifi-. 

K{q) ^ K*{q) ^ 1 + 20 - g + 36 - + 15176 - + 486564 - 

+ 21684020- cf' \ 1209684456- q^ + 58513394904- 

+ 2921860726948- + 141376772107064- -f - - - (137) 

which is Numher\2)l in tables [91]. 

We give, in Appendix L, the expansion of the Yukawa coupling for a set of other 
iim,n that are Calabi-Yau: in particular their exterior square is order five (not six as 
one could expect for a generic irreducible order- four operator). It will be shown, in 
a forthcoming publication, that the fact that the order five exterior power property 
occurs means that these operators are necessarily conjugated (by an algebraic function) 
to their adjoints. Thus, the "adjoint Yukawa coupling" K*{q) is necessarily equal to 
the Yukawa coupling K{q) for these operators. 

Remark: The operator having the Hadamard product of the two HeunG 
functions HeunG{a, q, 1, 1, 1, 1; x) and HeunG{A, Q, 1, 1, 1, 1; x) as a solution 
reads: 

[x -l){x- a) {x - A) (.T ~Aa){Aa~ x^f ■ ■ (138) 
+ 2 {x^ ^ Aa) -Ua- x^ ■ Dl - U2 ■ x ■ -Ui-D^ +Uo, 

where the polynomials Un are given in Appendix K. 

The exterior square of this order- four operator (138) is of order five for any value 
of the parameters a, q. A, Q (instead of the order-six one expects for the exterior 
square of a generic irreducible order- four operator). 

The HcunG hmctions solutions of the form HeunG{a, q, 1, 1, 1, 1; x) are 
an interesting set of HeunG functions. They verify the following (six Mobius) 
identities [124]: 

(1 q X 
— , — , 1, 1, 1, 1; — 
a a a 

1 ^z' a a — q x 
HeunGi , ^, 1, 1, 1, 1; - 



\ — X \a— I'a— 1' '''' I — X 

^ HeunG(^—, ^— 1, 1, 1, 1; —) (139) 



1 — x/a Vl— a'a— 1' 



HeunGi 1 — a, 1 — g, 1, 1, 1, 1; 



(a-1) 



a — x \ a — X 

1 { a — 1 a — q x a— 1 

= -; HeunG[ , , 1, 1, 1, 1; 

1 — X \ a a 1 ~ x a 

§ They are conjugated to their (formal) adjoint by a function. 

ft Actually 1/4,4 is exactly the operator for 4-F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1], 256 x). 
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The ten linear differential operators denoted by a star * in the previous table are 
all of this form: they have the Hadamard product of two HeunG functions solutions 
of the form HeunG{a, q, 1, 1, 1, 1; a;) as a solution. Note, however, that this HeunG- 
vicwpoint of the most interesting Hm,nS does not really help. Even inside this 
restricted set of HeunG functions solutions of the form HeunG{a, q, 1, 1, 1, 1; x) 
it is hard to find exhaustively the values of the two parameters a and of the accessory 
paramctcrf q such the series HeunG{a, q, 1, 1, 1, 1; x) is globally boundcd:|:. or, just, 
such that the order- two operator having HeunG{a, 1, 1, 1, 1; a;) as a solution is 
globally nilpotent (see Appendix K). 

The order-four operators H3 4, are all MUM operators^, but, similarly to 

the situation encountered with -ff2,2, their nome, mirror map and Yukawa couplings 
are not globally bounded. 

Many Hm,n are not MUM, for instance the order-eight operator -ffi2,i2, or 
the order-six operator H^j, are not MUM. Concerning H^ r and as far as its six 
solutions are concerned, it is "like" the four solutions of an order-four MUM operator, 
together with the two solutions of another order-two MUM operator, but the order- 
six operator H^j is not a direct-sum of an order-four and order-two operator. We 
have two solutions analytical at a; = (no ln(a;)), two solutions with a ln(a;). A 
linear combination of these two solutions analytical at x ~ is, by construction a 
series with integer coefficients (the Hadamard product of the two series with integer 
coefRcients which are the initial ingredients in this calculation), when the other linear 
combinations are not globally bounded. 



8.6. "Atkin" transformations 

It is worth noting that the globally bounded character of some of these Hm,n, at a; = 0, 
is simply related to the globally bounded character at x = cx) of a conjugate operator, 
as a consequence of some simple homomorphism relation with their transformed by 
"Atkin" puUbacks: 

224 X 



x^'^- Hs^six) = H3,3(^x^ ^) • (140) 



but also for the order-six operators 



X 



X 



V2.i/„(,) . (141) 
. HrAx) = Hr.7{x ^ ■ x^^ (142) 



f The accessory parameter appears in many applications as a spectral parameter [124]. 

I Along this line see the paper by Zagier [96] on integral solutions of Apery-like equations. 

^ Note that the Hadamard product of two MUM ODEs is not necessarily a MUM ODE: the order-six 

operator H^ ^ is not MUM. 
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and the order-eight operators 



x'^ ■ -^12,12(2;) = -^12,12(2; 



1 



-^10,10(^2; 



or the order-ten operators 

X^/^ ■ i?10,10(.T) 
x'/^ ■ Hi3S3{x) 

or the order-twelve operator 

x^ ■ His,is{x) : 
or the order-fourtecn operator 

x^^^ ■ H25MX) = ^^25,25(2:^ 
to be compared with 



2432 X 
1 



26 



1 



2852 X 

1 



.3/2 



2833132 X 

1 



,7/6 



2232 X 

1 



nl/4. 



^2(2;) 



1 



nl/2 



aj4(x) 



tLi4( X 



212 X 

1 

^ 28^ 

1 

Wx 



X ■ ivsix) 
and: 



ujsix 



r.l/2 



,2/3 



,1/4 



,1/2 



■ X, 



2852 X/ 

• Loaix) 

X ■ UJq{x) - 
X ■ i^Jg{x) 



,5/2 



(143) 
(144) 

(145) 
(146) 

(147) 

(148) 



cjs \x 
ujcJx - 



1 

36^ 
1 



X 



1/3 



UJq\X 



23 322; 
1 

33" a; 



^ 3^) 



^5(2;) 

IU7(.t) 



UJr,\X 



1 



UJ7 X 



26 53 X 
1 



,1/2 



,2/3 



■ X, 

X, 

(149) 



3472 x. 

Note, however, that these properties (140) are no longer valid for the ofF-diagonal 
operators Hm,n, m ^ n. For instance, the other order-four Calabi-Yau operators 
1/4,6, -^4,87 Hi g are not globally bounded at x = 00: if one changes these operators 
by a X — > 1/x pullback, the new order-four operators after a well-suited conjugation 
by x^ (r rational number) are not globally bounded. For instance H^ q^x — >■ 1/x) 
requires to be conjugated by ^1/2 (Puiseux series): after conjugation by xi/2, the 
series analytical at x = is not globally bounded. 

Remark: Let us denote An the constant in the "Atkin" involution x An/x 
for modular form of order n (i.e. t — > nr). We denote by a;,i(x), the order-two 
operator associated with a modular form of order n. 

One has 

1 



x*"" • a;„(x) 



An Bl 



(150) 



where the i?„'s are integers, and, quite remarkably^, one also has for iJ„_„(x) the 

*l If one forgets, for a second, that sum and multiplication do not commute, it is tempting to see 
naively these conjugation results (151), on the Hn,nS, as a straight consequence of the relation, like 
(150), on the w^'s, since the H„^„'s are Hadamard squares of the uin's, the constant in the pullback 
involution in the right-hand side of (151) being the square of the constant in (150). This is not the 
case. 
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(151) 



n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


12 


13 


16 


18 


25 


An 


2i-2 


3^ 


2« 


5^ 


2^ 3^ 


7^ 


2^ 


3^ 


2^5 


2^3 


13 


2^ 


2 • 3 


5 


Bn 


1 


1 


1 


2" 


1 


3^ 


1 


2 


1 


1 


2^ 3^ 


1 


1 




rn 


1 

4 


1 

3 


1 

2 


i 

2 


1 


2 
3 


1 


1 


3 
2 


2 


7 

6 


2 


3 


5 
2 



iS. 7. 2-^1 ([1/A^, 1/-/V], [1], a;) hypergeometric functions 

Modular forms can always be written as the hypergeometric function 2^1 j up 
to an algebraic pre- factor, and up to a pullhack (see Maier [123]). Relations 
(130), (132), (134) and (135) underline the special role of modular forms that 
can be written as 2F1 with no algebraic pre-factor, and no pullbacks (see 
Table 15 in [123]), namely 2-Fi([l/4, 1/4], [1], -64a;), 2Fi([l/3, 1/3], [1], -27x), 
2i^i([l/2, 1/2], [1], —16a;). Along this line it is worth considering the hypergeometric 
functions 2Fi{[l/N, 1/iV], [1], —N^x) which always yield globally bounded series, 
together with their associated linear differential operators. The nome for these 
operators does not (gencrically) correspond to globally bounded series. One notes, 
however, that iV = 6 is "special", yielding series with integer coefficients for the 
hypergeometric functionf 



.^^r([J,^],[l], 



-432 a; 



(152) 



(l+ 432a;)-i/6.^^^([_L^^]jl]^ 



1728 X 



(1 + 432 x)2. 

= 1 - 12.x + 1764.t2 - 397488x3 107619876 a;'' - 32285962800 
-f 10342270083600 x^ - 3467404345579200 x^ + ■■■ 
but also for the corresponding nome 



q ^ X -120x2 +24660x3 - 6322720 x'' 1828573410 x'^ - 570359919024 x*^ 



-t- 187363061411720 x^ - 63912709875600960 x^ + ■■■ (153) 

as well as the mirror map: 

x{q) = q +l2Qq^ +4140(j3 + 166720 q'* - 6012210 + 1165528224(7*' 

- 178811454280 +29512658112000(7^ + ••• (154) 

In fact the hypergeometric function 2-^1 ([1/6, 1/6], [1], x) (as well as the 
hypergeometric function 2^'i([l/6, 5/6], [1], 1 — x)) is actually "special" (as was first 
seen by Ramanujan, see for instance Cooper [125]). In Appendix M the modular form 
character of (152) is made crystal clear in a quite heuristic way. 

Remark: Using the 2^'i([a, &], [1], a;) hypergeometric functions associated 
to modular forms]: (namely [a, fo] = [1/2, 1/2], [1/3, 1/3], [1/4, 1/4], [1/6, 1/6], 

t We use identity 2-Fi( [1/6, 1/6], = {1 - x)-^/^ ■ 2-Fi{[l/12, 5/12], [1], -4a;/{l - x)^), which 

singles out the known [113] puUback x — > — 4x/(l — x)"^ . 

t For [1/2, 1/2], [1/3, 1/3], [1/4, 1/4], [1/6, 1/6], see the Ramanujan's theories to alternative 
basis [125] for other "signatures". For [1/8,3/8] see the third example in Appendix A. 
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[1/3,2/3], [1/3,1/6], [1/6,5/6], [1/4,3/4], [1/8,3/8], [1/12,5/12], ...), one can build 
4F3 globally bounded examples by simple Hadamard products of these selected 2F1. 
We give in Appendix P a miscellaneous set of identities expressing HeunG functions, 
or modular forms, as the previously selected 2-F1 hypergeomctric functions with two 
pullbacks. 

8.8. Modularity and Hypergeomctric series with coefficients ratio of factorials 

As a consequence of the classification by Beukers and Heckman [111] of all algebraic 
„_F„_i's, the 8-PV hypergeomctric series 



1 2 3 4 1 2 1 

5' 5' 5' 5' 2' 3' 3 



2" 39 55 X 



1 7 11 13 17 19 23 29 

30' 30' 30' 30' 30' 30' 30' 30_ 
has integer coefficients, and is an algebraic function. The Galois group belonging to 
this function is the Weyl group W{E^) which has 696729600 elements [126]. It is an 
algebraic series of degree 483840. More precisely, it was noticed [127] by Rodrigucz- 
Villegas that the previous power series reads: 

E • (155) 



n=0 



(15n)! (lOn)! (6n)! 



which is precisely the series introduced by Chebyshev during his work [128] on the 
distribution of prime numbers to establish the estimate 

0.92 — ^ < 7r(x) < 1.11 — (156) 
log X log X 

on the prime counting function t:{x). 

Considering hypergeometric series such that their coefficients are ratio of 
factorials, a paper by Rodrigues and Villegas [127] gives the conditions of these 
factorials for the hypergeometric series to be algebraic (all the coefficients are thus 
integers). A simple example is, for instance the algebraic function: 

^ /.I 1 3. .1 2. 256 \ ^ (An\ „ 

71 — ^ ^ 

Along this line it is worth recalling Delaygue's Thesis [129] (sec also Bober [130]) 
which gives some results!] for series expansions^ such that their coefficients are ratio 
of factorials: 

4^^3([i i],[l,l,l];256..) = E^^-", (159) 



/rl 1 1 5, , , s q \ ^ (671)! (2n)! 

^^3([2'2'6'6]'r^'^'^]^ 2 ^ ^ EUw'- - (^^0) 



rl 1 1 5^ , ,1 ,8 o3 _ V- MU2n)! 

- (3n)!(n!) 

These ratio of factorials are integer numbers. The series expansion of (160) reads: 
1 + 240.T + 498960 a;2 + 1633632000x3 + 6558930378000 

+ 29581300719210240 x^ + 143836335737833939200 x^ + • • • 

J Necessary and sufficient conditions for the integrality of tlic mirror maps series. 
^ These series are not algebraic functions. 
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8.9. More Hadamard products: Batyrev and van Straten examples [90] 
8.9.1. A first auto-adjoint Calabi-Yau ODE 

An order-four operator has been found by Batyrev and van Straten [90] 

Bi = -3x-{7e'^ +76 +2)-{39 + l)-{39 + 2) (161) 

- 72a;2 • (3 61 + 5) • (3 61 + 4) • (3 61 + 2) • (3 61 + 1), 

which is conjugated to its adjoint: Bi ■ x — x ■ adjoint(Bi). 

Operator (161) is a Calabi-Yau operator: it is MUM, and it is such that its exterior 
square is of order five. Its has a sohition analytical at x = which is actually the 
Hadamard product of the previous selected hypergeometric 2F1: 

' 'l>W^H*(T^-^'(i|l'W;7Tfi^))^ (1^2) 



2 



^^1 



^3' 3^'^ " J VI +4a; " 'V^3' 3" ' " (1 +4x)3 

8.9.2. A second auto- adjoint Calabi-Yau ODE 

A second example [90] (seef page 34) of order-four operator corresponds to Calabi- 
Yau 3-folds in Pi x Pi x Pi x Pi: 

B2 = 6** -4a:- (56i2 +561 +2) • (261 + 1)2 

+ 64x2 • (261 + 3) • (26* + 1) • (26' + 2)^, (163) 

corresponding to the series-solution with coefficients: 



f2n\ ^ fnV {2k\ {2n-2k 



\ n I ^-^ \k J \ k J \ n — k 

^ ' k=o ^ ^ \ / \ 



2nV ^ /A , , 2n - 1 



(164) 



, ^Piiin^ -n, ~n\, [1, 1, ]; 1 

n / V 2 2 

Its Wronskian reads: 

w — i 0.3 , wi/2 _ \ ('I65'l 

(1 -64x)2 . (1 - 16a;)2 . a;6' (1 - 64 x) (1 - 16 x) ' ^ ' 

This operator is also a Calabi-Yau operator: it is MUM, and it is such that its exterior 
square is order five. This order five property is a consequence of B2 being conjugated 
to its adjoint: B2 ■ x = x ■ adjoint{B2). 

The series-solution of (163) can be written as an Hadamard product 

S = (1 -4x)-^/2 ^ i?eunG(4, 1/2,1/2,1/2,1,1/2; 16a;)2 (166) 

= 1 + 8 a; + 168 a:^ + 5120 x^ + 190120 x^ + 7939008 x^ + 357713664 x^ 4- • • • , 

the modular form character of HemiG{4:, 1/2, 1/2, 1/2, 1, 1/2; 16 x) being illustrated 
with identities (A. 2) in Appendix A. Its nome reads: 

q ^ x + 20 a;2 + 578 + 20504 x* + 826239 x^ + 36224028 x^ + 1684499774 x^ 

+ 81788693064x^+4104050140803x^+211343780948764x1° + ••• (167) 

Its mirror map reads: 

x{q) = q ~2Qq^ + 222 - 2704^"^ +21293(7^ - 307224 + 80402 

- 67101504 - 1187407098(7^ - 37993761412 +••• (168) 

t There is a small misprint in [90] page 34: (2 9 + 1) must be replaced by (2 6 + 1)^ in the ix term. 
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The Yukawa coupling of (163) reads: 
K{q) = K*{q) = 1 +4g +164g2 + 5800 + 196772 + 6564004 

+ 222025448 + 7574684408 q'^ + 259866960036 q^ + ■■■ (169) 

The equality of the Yukawa coupling with the "adjoint" Yukawa coupling, K{q) = 
K*{q), is a straight consequence of relation B2 ■ x = x ■ adjoint{B2). 

Do note that recalling Batyrev and van Straten [90], (see step2 page 496), and 
following Morrison [89], one can also write the Yukawa coupling as: 

Vl Vx(g) dq ) yl V dq ) ' ^''"^ 

where W4 is the Wronskian (165). From this alternative expression for the Yukawa 
coupling it is obvious that if the analytic series 2/0(2^)1 "^^^^ the nome (167) are 
series with integer coefficients, then, the mirror map (168) is also a series with integer 
coefficients, and, therefore, j/o seen as a function of the nome q, as well as x^ W^^^ , 
and, consequently, the Yukawa coupling is a series with integer coefficients (as a series 
in q or in x). 

The globally bounded character of the analytic series yo (x) together with the nome, 
thus yields the globally bounded character of the mirror map, Yukawa coupling, that we 
associate with the modularity. Similar results can be found in Delaygue's thesis [129], 
in a framework when the coefficients of hypergeometric series are ratio of factorials. 

In contrast the globally bounded character of the analytic series yo{x), together 
with the globally bounded character of the Yukawa coupling (seen for instance as a 
series in x) does not imply that the nome, or the mirror map, are globally bounded 
as can be seen on example (113) (see (114) and (115)). 



8.9.3. An operator non trivially homomorphic to B2 

Let us, now, introduce the order- four operator 

62= 256a;2- 6|2(2 6' + 3)(2 6'+l) 

- 4x- (261 + 1) (261-1) (56|2- 561 + 2) +{9-1^. (171) 

This operator is non-trivially:|: homomorphic to the Calabi-Yau operator (163): 

32- X- (26 + 1) = x-{29 +!)■ B2. (172) 

As a consequence of the previous intertwining relation, one immediately finds that the 
series-solution analytic at x = of this new MUM operator (171) is nothing but the 
action of the order-one operator x • (2 6* + 1) on the series (166), and reads: 

x-{2e+l)[S] = X -1-24x2 -f 840 +35840 a;'' + 1711080 -I- 87329088 a:^ 
+ 4650277632 + 254905896960 + ••• (173) 

It is obviously also a series with integer coefficients (the action of a; • (2 + 1) on the 
series with integer coefficient is straightforwardly a series with integer coefficients). 
More generally, the globally bounded series remain globally bounded series by 
operator equivalence (non trivial homomorphisms between operators: generically the 
intertwiner operators are not simple functions). 

I The intertwiners between B2 and B2 arc operators not simple functions. 
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The exterior square of the order-four operator (171) is an order-six operator which 
is, in fact the LCLM of an order-five operator £5 and an order-one operator: 

Ext\B2) = ® (d, - where: p{x) 



dx J' " ' (1 - 16x)(l -64a;) 

Operator B2 is non-trivially homomorphic to its adjoint: 

32- x^ ■ {20 +3) ■ {29 +5) = x^ ■ {2 9 +3) ■ {2 +5) ■ adjoint{B2). (174) 

The Yukawa coupling of this order-four operator (171), non-trivially homomorphic to 
(163), reads: 

K{q) = 1 -4g -140g2 -4040g^ -64436-^ + 1889332 -i- 

4- 88331368 — + 1652707624 69295027684 — + ■■■ (175) 

5 9 63 

The Yukawa couphng series (175) is not globally bounded. 

The "adjoint Yukawa couphng" of this order-four operator (171) reads: 

K*{q) = 1 + 12 g + 564^2 + 20440 + 865732 g'' + 37162444 

+ 8255346664 ^ -M 121762648248 ^ -f 72336859374772 ^ + ••• (176) 
5 15 21 

Again, the adjoint Yukawa couphng series (176) is not globally bounded. 

On this example one sees that the Yukawa coupling of two non-trivially 

homomorphic operators are not necessarily equal. The Yukawa couplings of two 

homomorphic operators are equal when the two operators are conjugated by a function 

(trivial homomorphism) . The modularity property is not preserved by (non-trivial) 

operator equivalence: it can depend on a condition that the exterior square of the 

order-four operators are of order five. The Calabi-Yau property is not preserved by 

operator equivalence. 

To sum-up: All these examples show that the integrality (globally bounded 
series) is far from identifying with modularity. All these examples have to be taken 
into account if one has in mind to build new conjectures combining these globally 
boundedness of various series with the concept of diagonal of rational functions: for 
instance, can wc imagine that being a diagonal of rational functions automatically 
yields that the nome or the Yukawa coupling are globally bounded series in q or x, 
etc ? 



9. Conclusion 

Seeking for the linear differential operators for the x''"'''s, we first discovered that 
they were Fuchsian operators [14, 16], and, in fact, "special" Fuchsian operators, 
namely Fuchsian operators with rational exponents for all their singularities, and 
with Wronskians that are 7V-th roots of rational functions. Then we discovered that 
they were G-operators (or equivalcntly globally nilpotent [28]), and more recently, 
we accumulated results [33] indicating that they arc "special" G-operators. There 
is, in fact, two quite different kinds of "special character" of these G-opcrators. On 
one side, we have the fact that one of their solutions is not only G-series, but is a 
globally bounded series. This special character has been addressed in this very paper, 
and we have seen that, in fact, this ^'integrality property [31]" is a consequence of 
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quite general mathematical assumptions often satisfied in physics (the integrand is 
not only algebraic but analytic in all the variables (50)). However, we have also 
seen another special property of these G-operators, namely the fact that they seem 
to be quite systematically homomorphic to their adjoints [33]. We will show, in a 
forthcoming publication, that this last property amounts, on the associated linear 
differential systems, to having special differential Galois groups^ and that their exterior 
or symmetric square, have rational solutions. This last property is a property of a 
more "physical" nature than the previous one, related to an underlying Hamiltonian 
structure [121], or as this is the case, for instance in the Ising model, related to the 
underlying isomonodromic structure in the problem, which yields the occurrence of 
some underlying Hamiltonian structure [121]. In general the integrality of G-operators 
does not imply the operator to be homomorphic to its adjoint, and conversely being 
homomorphic to its adjoint does not imply:]: integrality (and even does not implyf 
the operator to be Fuchsian). Interestingly, the x^"-''s, as well as many important 
problems of theoretical physics, correspond to G-operators that present these two 
complementary "special characters" (integrality and, up to homomorphisms, sclf- 
adjointness), and, quite often, this is seen in the framework of the emergence of 
"modularity" . 

Nomes, mirror maps, and Yukawa couplings are not D-flnite functions: they 
are solutions of quite involved non-linear (higher order Schwarzian) ODEs (see for 
instance Appendix D in [28]). Therefore, the question of the series integrality 
of the nomes, mirror maps, Yukawa couplings, and other pullback-invariants (see 
Appendix J) requires to address the very difficult question of series-integrality for 
(involved) non-linear ODEs, or, equivalently, the problem of non-linear recursions 
with integer sequence solutions. Note, however, as seen in Section 8.9.2, in particular 
in (170), that the integrality of the series yo{x) and of the nome q{x) are sufficient to 
ensure, provided the operator is conjugated to its adjoint (see (??)), the integrality 
of the other quantities such as the Yukawa coupling, mirror maps. However the 
integrality of the nome remains an involved problem. These questions will certainly 
remain open for some time. 

In contrast, and more modestly, we have shown that a very large sets of problems 
in mathematical physics (see sections (4.4), (5) and (6.1)) actually correspond to 
diagonals of rational functions. In particular, we have been able to show that the 
-fold integrals of the susceptibility of the two-dimensional Ising model are 
actually diagonals of rational functions for any value of the integer n, thus proving 
that the x^"^ 's are globally bounded for any value of the integer n. As can be seen in 
the "ingredients" of our simple demonstration (see (4.4)), no elliptic curves, and their 
modular forms [132], no Calabi-Yau [118], or Frobenius manifolds [121], or Shimura 
curves, or arithmetic lattice assumption [119, 120] is required to prove the result. We 
just need to have a n-fold integral such that its integrand is not only algebraic, but 
analytic in all the variables. 

Acknowledgments We would like to thank A. Engc and F. Morain for 

I See Appendix N which gives an example of a (hypergeometric) family of order-four operators 
satisfying the Calabi-Yau condition that their exterior square is of order five, and, even, a family of 
self-adjoint order-four operator, the corresponding hypergeometric solution-series being not globally 
bounded. See also Appendix O. 

f For instance the operator — x Dx — 1/2 (see page 74 of [131]) with an irregular singularity is 
self-adjoint. 
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Appendix A. Modular forms and series integrality 
First example: The generating function of the integers 



(A.1) 



-A /nV /2fc\ /2n-2fc 
t'o W ' UJ ' V n-k 

^ I • 2-f^i(^[2' """J' 2 ^' 

is nothing else but the expansion of the square of a HeunG function 

i7eOTiG(4, i i i 1 i; 16-a;) = 1 +2x + 12 a;^ + 104a;^ 

+ + 12348 + 150528 + 1914432 + ••• (A.2) 

solution of the order-two operator 

Hdiam = 0^ -2-x-{lO0^ +59 +1) + 16x2 • (2 61 + 1)2. (A.3) 

which corresponds to the diamond lattice [48]. This HeunG function (A.2) is actually 
a modular form] which can be written in two different ways: 

HeunG{A, i, i, i, 1, i; 16x) 

. (i-4.)-/v.^.([i.i,,|il;^) (A^4, 

, „ ^ 1 /9 ^ /r 1 1 1 r 1 108 X 

= (1 -16x)-l/2 . 2J^i [- -], [1]; 



/6' 3" (1 - IQxf, 

These two puUbacks are related by an "Atkin" involution x ^ 1/64/x. The 
associated modular curve, relating these two pullbacks (A.4) yielding the modular 
curve: 

4-y^z^ - 12y^ ■ (y + z) + 3y z ■ (Ay^ + A z^ - 127 y z) 

- A- {y + z) ■ {y^ + z^ +83yz) +432yz 0, (A.5) 

which is (y, z)-symmetric and is exactly the rational modular curve in eq. (27) already 
found for the order-three operator in [33] for the five-particle contribution x'^-* of 
the magnetic susceptibility of the Ising model. 

This result in [46, 48] can be rephrased as follows. One introduces the order-three 
operator which has the following 3F2 solution 

(4-^^- ^^4[3'2'3]'[^'^]' (4-^)' (^-^^ 

f Gcncrically HeunG functions arc far from being modular forms. They are even far from being 
solutions of globally nilpotent operators (they gcncrically have no integral representations [133, 134]). 
There is a relation between these operators being finite-gap [135] and their globally nilpotence. 
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associated with the Green function of the diamond lattice. Along a modular form line 
lets us note that this hypergeometric function actually has two puUbacks: 

- 4 /,! 1 2_ , 27 x'^ 



4. (a;2 _ 1) ^V^3' 2' 3^' ^ ' J' 4- {x^ - 1)3 
These two puUbacks related by the "Atkin" involution x — > 2/x: 

"(^) = (4-^x^)3 ' ^(^) = "0 = 4. (1^^-1)3 ' (^-^^ 

corresponding, again, to the modular curve (A. 5). 
Second example. The HeunG function 
HeunG{~3, 0, 1/2, 1, 1, 1/2; 12 • x) (A.9) 

= (I +Ax)-^/'^- HeunG(A, 1,-, - 

^ ^ ' V ' 2' 2' 2' ' 2' 1 



16 a; 



4x/ 
+ 1! 

1823094 a;* + 18406752 a;^ + 189532980x1° 



1 +6x^ +24x''^ + 252x'' +2G16x^ + 19320 x^ + 183456 x^ 



is solution of 

Heuufcc = 0"^ -2x-6'- (46* + 1) - 24 • x^ • (2 6i + 1) • (6* + 1), (A.IO) 

The square of (A.9) is actually the solution of an order-three operator (see equation 
(19) in [48]) emerging for lattice Green functions of the face-centred cubic (fee) lattice 
which is thus the symmetric square of (A.IO). This hypergeometric function with a 
polynomial pull-back can also be written: 

HeunG{-i, 0, 1/2, 1, 1, 1/2; 12 • x) 

= 2^^i([^, J],[l]; 108-x2-(l+4x)) (A.ll) 

^^ io ^-l/2 r:,(A li ni 108-x- (1-H4x)2n 

= (l-12x) '•2-^1-,-,!; ^ — , 

^ ' V^6' 3^' ^ ^' (1 - 12x)3 /' 

where the involution x ^ — 1/4 • (1 +4 x)/(l — 12 x) takes place. The modular curve 
relating these two puUbacks reads exactly the rational curve (A. 5) already obtained 
in [33]. 

Third example. The HeunG function i7eunG'(l/9, 1/12, 1/4,3/4, 1, 1/2; 4x) is 
solution of the order-two operator corresponding to the simple cubic lattice Green 
function 

H,c = 0"^ - X • (40 -f- 20 61 4- 3) + 9 • x^ • (4 61 -h 3) • (4 6* 1). 

The square of this HeunG function is a series with integer coefficients which identifies 
with the Hadamard product of (1 — 4x)~i/^ with a modular form : 

HeunGil/9, 1/12, 1/4, 3/4, 1, 1/2; 4x)2 (A.12) 
= (1 - 4x)"i/2 * HeunG{l/9, 1/3, 1, 1, 1, 1; x) 

= 1 +6x -1-90x2 1860x3 -f 44730 x'' 1172556 x^ + 32496156 x^ 

4- 936369720 x^ -h 27770358330 x* -f 842090474940 x^ + ■■■ 

The HeunG function HeunG{l/9, 1/12, 1/4, 3/4, 1, 1/2; 4x) is globally bounded: the 
series of i?eMnG'(l/9, 1/12, 1/4, 3/4, 1, 1/2; 8 x) is a series with integer coefficients. 



Diagonals of rational functions 50 

One can also write this HeunG function in terms of a 2^1 ([1/6, 1/3], [l],a;) 
hypergeometric function up to a simple algebraic pullback (with a square root), or 
in terms of a 2-F'i([l/8, 3/8], [l],a;) hypergeometric function: 

i?ewiG'(l/9,l/12,l/4,3/4,l,l/2;4a;) = Ca^" • ([1/8, 3/8], [1]; P2) , with: 



= 9.(l|l2.)^ -(^-^^-+^-^^-^^-)^^^)' - (TTtSf'^- 

P2 = (1 -42 a; + 352x2 - 288x3) + (1 _ 4 3;) . (1 _ 20 x) • (1 - 36 .t)1/^ 

Do note that taking the Galois conjugate (changing (1 — 36 x)^/^ into —(1 — 36 x)^/^) 
gives the series expansion of 3-^2- ffeunG(l/9, 1/12, 1/4, 3/4, 1, 1/2; 4x). This 
shows that there exists an identity for 2^i([l/8, 3/8], [l],x) with two different 
pullbacks, namely the previous P2 and its Galois conjugate, these two pullbacks being 
related by a (symmetric genus zero) modular curve: 

5308416 • z'' + U2368 ■ ■ {y + z) + 512y^ ■ {27 y^ + 27 z^ - 2737Axy) 
+ 192 y z • (y + z) • (y^ + z^ + 10718 y z) + / + z** + 3622662 y^ 
- 19332 • yz • (?/2 + z^) + 79872 •yz-(y +z) - 65536- yz = 0. (A.13) 

Revisiting the examples. In a recent paper [136] corresponding to spanning 
tree generating functions and Mahler measures, a result from Rogers (equation (36) 
in [136]) is given where the two following ^F^ hypergeometric functions take place: 

,,F.([|,2,I,l,ll,|2,2,2,21.^,,^ 

I 1,11,12.2.2.2], (A.14) 

The corresponding order-five linear differential operators (annihilating these two 5F4 
hypergeometric functions) are actually homomorphic (the intertwiners being order- 
four operators). The relation between these two pullbacks y = 256x3/9/(x + 3)'^ 
and z = 256x/9/(l +3x)*, remarkably yiwes, again, the previous (y, z)-symmetric 
modular curve (A.13). 

The order-five linear differential operator, corresponding to the first 5i^4 
hypergeometric function, factorizes in an order-one operator, an order-three operator 
and an order-one operator, the order-three operator being, in fact, exactly the 
symmetric square of an order-two operator: 

• • (x-9)t+3)^ ■ ""'^ 

where the order-one operators read respectively 

dx V(9x2 + 14x +9) • (x + 3)4/' dx V x3 

and where the order-two operator W2 reads: 

_ „2 o (6- x^ +7x +3) 3 3x +2 

~ +^9x2+14x+9)-x'^^ +4' (9x2 +14x +9). X- (^.15) 

We have a similar result for the order-five linear differential operator 
corresponding to the second ^F^^ hypergeometric function. 
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Another solution of this order-five hncar differential operator reads: 
(a; +3)4 f x-9 ^ I 3, ,^ ^, 256x3 n ^ 



(x + 3)- X 2' 4^' ^ ' "9- (.T + 3)4, 

The expansion of the 3F2 hypergeometrie funetion in (A. 16) is globally bounded 
(change a; — >■ 9 x to get a series with integer coefficients) . 

Recalling the two previous pullbacks we have, in fact, the following identity: 

/rl 1 3, , , 256 x^ 
3- (1 + 3x). 3F2 [-, -, -], [1, 1], 



4' 2' 4^' ^ ' "9- (x + 3) 

= + [Ml, ^^^). (A.17) 

However this 3F2 hypergeometrie function is nothing but the square of a 2F1 
hypergeometrie function 

3i^2([J, ^, J], [1, 1], x) = 2Fi([i [1], x)'. (A.18) 

Thus, the previous identity (A.17) is nothing but the identity on a 2F1 hypergeometrie 
function with two different pullbacks: 

A 3, 256x3 



(l + 3.x)i/2. ^fJ[-,\[11 -^P^) 



9 • (x + sy 

1 + 3) -^^Os'S^'M' 9.(l + 3.x)0 - (^-19) 

The expansion of (A.19) is globally bounded. One gets a series with positive integer 
coefficients using the simple rescaling x —J- 36 • x. Note that the two pullbacks can be 
exchanged by the simple "Atkin" involution .x o 1/x, being related by the modular 
curve occurring for the simple cubic lattice, namely (A. 13). 

We have a similar result for the other 5^4 hypergeometrie functions popping out 
in [136]. 

For instance, for the diamond lattice one gets an expression (see eq. (50) in [136]) 
where the two following ^F^ hypergeometrie functions take place:]:: 

,5 3 4 , , , 27x4 \ ^ 

3, 2^3, 1,1], [2, 2, 2, 2], j^—^y (A.20) 

These two pullbacks can be exchanged by the simple "Atkin" involution x O 2/x. 
These two pullbacks have been seen to be related by the (genus-zero) (y, z)-symmetric 
modular curve (A. 5): 

-12y^ ■ {y + z) + 3y z {Ay^ + A z'^ - 127 y z) 

- A- {y + z) ■ {y^ + z^ + 83yz) +432yz = 0. (A.21) 

Similarly to (A.17) we have an identity between two 3F2 hypergeometrie functions 
(namely 3F2([2/3, 1/2, 1/3], [1, l],z)) with the two pullbacks (A.20), and these 3F2 

J Note a small misprint in eq. (50) of [136]: one should read —272^/4/(1 — 2^)'^ instead of 
-272-1/4/(1 - 2^)3. 
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hypergeomctric functions being the square of 2F1 hypcrgeometric functions, one finds 
that the "deus ex machina" is the identity similar to (A. 19): 

The series expansion of (A.22) is globaUy bounded. Rescahng the x variable as 
a; — >■ 4x, the series expansion becomes a series with positive integer coefficients (up 
to the first constant term). 

For the face-centred cubic lattice one gets an expression (see eq. (52) in [136]) 
where the two following hypergeometric functions take placef: 

/r5 3 4 , , , a;^ • (x + 3)\ , , ^ 

'^^[3' 2' 3' ^' ^' ^' 4 )■ ^^-^^^ 

This example is nothing but the previous diamond lattice example (A. 20) with the 
change of variable x — > -~3x'^/{x'^ — 4) in (A.23). Therefore, the two pullbacks in 
(A.23) are, again, related by the modular curve (A. 5). The two pullbacks in (A.23) 
can actually be seen directly in the following identity (equivalent to (A.22)): 

/ I 1 X- jx +3)^ _ 2N1/2 p (A li Ml ■ +3) \ 

Finally, the equation (17) of [136] on Mahler measures, the two following 4F3 
hypergeometric functions take place: 

.f,([y,l,ll,|2,2,21, 

.f,([y,l,ll,|2,2,21,^). (A.24) 

These two previous pullbacks can be exchanged by an "Atkin" involution x ^ —8/x 
and are related by the (genus-zero) (y, z)-symmetric modular curve: 

8y^z^ -I2y^z^ ■ {y + z) + 3y z ■ {2y^ + 2 + 13y z) 

- {y + z) ■ {y^ + +29y z) +27 yz = 0. (A.25) 

The underlying identity on 2F1 hypergeometric functions with the two pullbacks 
(A.24) read: 

-2-{x-2)- 2fJ[11],[11 



3' 3^' ' " (x + 4)= 

The series expansion of (A.26) is globally bounded. Rescaling the x variable as 
X ^ —8 a;, the series expansion becomes a series with positive integer coefficients. 

t There is one more misprint in [136]: the puUback —x{x -f 3)/{a; — l)** must be changed into 
x{x + 'i)/{x - if. 
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Appendix B. Another logarithmically bounded series 

Let us display other logarithmieally bounded series than (2). The hypergeometric 
function 2Fi{[N/3, 1/6], [7/6], 9x) is not globally bounded but is such that the order- 
one operator 69+1 acting on it, is a series with integer coefficients for every integer 
value of N: 

{6 9+1)[,F,{[^,1]aI],9^)] = i^^o([f], [],9x) 

= l + 3x+18a;2 + 126x^ + 945a;'^ + 7371x^ + 58968x^ + ••• for = 1, 

= 1 + 6. 1 + 45x2 + 360 + 2970 x'' + 24948a;5 + 212058 + ••• for = 2, 

= l + 9x + 81x2 + 729x^ + 6561x'^ + 59049x^ + 531441x^ + ••• for TV = 3, • • • 
Similarly 

U = 3^.([i 4i],[|^],64x), (B.l) 

which is such that the action of the order-one operator 7 6' + 1 changes it into a 
globally bounded function 

{79 + l){U) = 2i^i([J, ^], [^], 64x) = 1 + 7x + 190x2 +7068x3 
+ 303924 x"^ + 14208447 x^ + 701448594 x^ + 35983401900 x^ + ••• 

Appendix C. <f>^'^(w) as diagonals 

The family of simple integrals was introduced in §4 of [82], as a way to 

simplify the study of the singularities of the Ising integrals x''"''- By definition, they 
are equal to 

^ J 12 r'^'^ d(f) 1 

where 

x(w;,0) = (C.2) 

1 - 2m; cos(0) + ^(1 -2zi;cos(0))2 -4^2 

By an easy change of variables, it follows that 

= +4-*d'H. (C-3) 



where 



*^"'H = - / F^[w,t)- ^f=, (C.4) 



1 /"^ , , dt 



TT 



-1 vT^ 



2 ' 



the algebraic function Fn{w, t) being defined by 



where 



1 ~h[w,ty^ ^ ■h{w,Tn-i{t)) 
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and where Tjn{t) is the m-th Chebyshev polynomial of the first kind, that is, the 
unique polynomial of degree m such that cos{mt) ~ r,„(cost). 

In order to express ^^^\w) as the diagonal of an algebraic function in two 
variables, it is sufficient to use the following general result: 

// F^w, t) is a bivariate power series, then the univariate power series 

"1 Fjw, t) 

is the diagonal of the generalised power series^ 

The only non-trivial point in the proof of this fact is the classical integral evaluation 
^ 1 for \u\<l. (C.9) 



J-1 



dt (C.7) 



TT y_i (1 - ut) ■ VT~^ VT^ 



2 



Expanded proof: Letting F{w,t) ~ J2e>o /^(O ' "^^j the series 'if{w) is equal 

to 

while the series G{w, t) is equal to 



It follows that the diagonal of G is equal to 



Diag(G)H ^ Y.[tO]Ii^.^^. (C.12) 



f>0 

To prove that 4* = Diag(G'), it thus suffices to show that for any power series /, 

By linearity, it thus suffices to prove that for any no n- negative integer s, 
\ t^ 1 

-• / -^==-dt = [^1^==f (C.14) 
This follows from the classical integral evaluation 



i. [ '^^ = ^ for \u\<l. (C.15) 

example, when n = 2, the previous construction shows that 

= 1 16^^) (C.16) 



1 
2 

is equal to the diagonal of the algebraic function 



1 -2w + ^(l-2w)2 -AwH^ 1 
2 • ■ - 2w)2 - AwH^ ^ 2' 

§ In the sense of [58]. 



(C.17) 
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Appendix D. Creative telescoping: computing ODEs for diagonals 

The notion of diagonal of rational function is ubiquitous in combinatorics [137]. 
Its importance comes from the fact that many operations on power series with 
a combinatorial relevance (Hadamard products, constant terms, or positive parts, 
of Laurent power series, etc) can be encoded as diagonals. A classical result by 
Lipshitz [57] predicts that diagonals of rational functions are D-finite. The question 
is then: how to obtain algorithmically a differential equation satisfied by the diagonal 
Diag(/) of a given rational function /(xi, . . . , The question can be reformulated 
in terms of computing a multiple integral with parameters, over an algebraic surface 
("vanishing cycle" or "evanescent cycle" in Deligne's terminology [64]), and thus can 
be attacked from a geometric viewpoint. 

A first answer to this algorithmic question is provided by Lipshitz 's result [57]: if 
denotes the rational function F = /(xi, a;2/a;i, . . . , Xn/a;n-i)/(xi • • • x,i_i), and if 
the following equality, called the creative telescoping equation^ 

"9x„)'~^'* dxi ^ ^ dxn-i' ^^'^"^ 

admits a solution (i, gi, . . . , .gn-i), where P (called telescoper) is a linear differential 
operator with coefficients in Q[a:;„], and where . . . , Rn-i arc rational functions 
in Q{xi, . . . ,x„) (called certificates), then P annihilates the diagonal Diag(/) of /. 

Several algorithms exist for solving equation (D.l). A common weakness of 
currently known algorithms for solving (D.l) is that they are not able to compute the 
telescoper L without computing the certificates (i?i, . . . , i?„_i). This is unfortunate, 
since in practice only the telescoper is really needed, while the size of the certificates 
is much more important than that of the telescoper. 

Lipshitz 's initial argument requires the construction of a non-zero operator 
annihilating F which involves all the partial derivatives d/dxi. This reduces the 
resolution of (D.l) to that of a linear system over Q. The big practical issue with this 
approach is the size of the linear system, which is about several millions even for the 
simple rational function / = 1/(1 — xi — X2 — a^a). A much more efficient algorithm 
for solving equation (D.l) is Chyzak's extension [65, 138] of the Zeilberger's celebrated 
creative telescoping method [139], although the computational complexity of Chyzak's 
algorithm is not yet well understood. The most efficient implementation of Chyzak's 
algorithm is due to Koutschan [140]. 

More generally, for ?i-fold parameterised integrals of D-finite functions, Chyzak's 
creative telescoping algorithm delivers a system of PDEs. For two variables 
(anisotropic Ising model), one will get a system of PDEs corresponding to two 
"telescopers" , that can be written in the following form: 

^^i^'^'dx) ^ 51 ■ ■^^(•^'^'a^) " ^ Inix, y) ■ D"^. 

n—Q n—Q 

where the Pn's, and the Qn's, are polynomials of the two variables x and y. 

Note that, in practice, the down-to-earth physicist's guessing techniques we have 
used in our various papers [14, 16, 26, 82], which amount to getting]: the linear ODE 

X Strictly speaking, the correctness of the linear ODEs obtained by "guessing" is not mathematically 
guaranteed. However, one may be convinced on the correctness of the ODE, since, in practice, one 
has longer series than what is used in the guessing. Also, some properties (as global nilpotcnce, 
expected known structures) are retrieved. 
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from the series expansion of the diagonal (or, in general, of the parameterised integral) 
is mueh more efRcientf than the creative telescoping approach. This is moral, since 
time consuming computations are the price to pay in order to guarantee the correctness 
of the ODE. 

Appendix E. Christol's theorem in more heuristic terms 

Let us give a sketch in heuristic terms of how the main theorem in [49] is proved. 

The first step is purely algebraic-geometric. The algebraic function F involved in 
the integral representation 

f{x) = / F{x]Xi, . . . ,Xn) ■ dxi---dxn, (E.l) 
Jc 

lives on a complex (n + l)-nianifold V, or, more precisely, on a family of smooth 
complex n-manifolds Vx- One applies to it the so-called embedded resolution of 
singularities [142]. This process uses a succession of blowing up which is theoretically 
explicit but seems to be inaccessible for computation. 

Roughly^speaking, we so obtain a new family of manifold V with Vx = Vx for 
X 0, and Vq, a divisor with normal crossing, namely, a union of "smooth algebraic" 
n-manifolds Dj that meets "transversally" . In particular, if non void, an intersection 
Di^ n • • • n Di^ of m distinct U^'s is of (complex) dimension n — m + 1. It is obvious 
that, at most, n + I divisors I?i's can intersect at a given point of Vg. 

Moreover, if there are really n -I- 1 divisors Di intersecting at the point P of Vq, 
then, choosing P as origin and equations Xi ^ oi Di {0 < i < n) as new variables, 
the equation of V becomes, at least locally, Xo...Xn = x. Applying the (algebraic) 
change of variables (Ai, . . . , A„) = ip{xi, . . . , x„) to (E.l) (x coming in the picture 
as a parameter), one gets 

/(x) = / F(x;Ai,...,A„). ^ ... (E.2) 

for an algebraic function F. If we are lucky, the cycle f{C) is (honiotopic to) Cp, 
the vanishing cycle around P, and (E.2) is an avatar of (53). Then 

/ = Diag(F) with F(Ao, A„) = F(Ao . . . A„ ; Ai, . . . , A„), (E.3) 

(up to a multiplicative constant) and / is the diagonal of an algebraic function (in 
n + 1 variables), hence, the diagonal of a rational function (in 2n -I- 2 variables). 

Actually the computation of <p(C) is inaccessible. To find hypothesis under which 
it is possible to conclude, we turn to the Picard-Fuchs equation Ly because it does not 
depend on the cycle C. Moreover it is stable under birational maps like the blowing 
up used for the desingularisation. All reasoning we will do from now do concern all 
solutions of Ly and cannot be done by considering only, for instance, the minimal 
order linear differential equation of /. 

Reverting the process, we conclude from formula (E.2) that integration on Cp 
gives a solution of the Picard-Fuchs equation in the ring of diagonals of rational 
functions (proposition 11 in [49]). Let us recall that the Picard-Fuchs linear differential 

f Serious programming improvements of the creative telescoping method have been developed 
recently [140, 141], and it is now possible to get the linear ODE for the isotropic Ising model x'^'i 
from creative telescoping calculations. 



Diagonals of rational functions 



57 



equation is given by the derivation ^ acting (through derivation under the integral 
sign) on the space of n differentials modulo exact ones and a solution of this 

"differential module" is a C(x)-lincar application from this space to some function 
space (here the diagonals of rational function) that do commute with ^ . 

It is difficult to decide, a priori, whether, or not, such a solution is zero on 
the particular differential F{x; Xi . . . Xn)^Yr ' " begin with. But we can 

assert that it is non zero for differentials G{x; Xi, . . . , Xn)^^ ■ ■ ■ such that 
G(0, ...,0) ^ 0. So, we consider the P-residue, for P — Di Ci ■ ■ ■ Ci Dn, which, 
roughly speaking, associates to a given differential the coefficient of • • • it 
contains. Then the Poincare residue map associates to a differential the family of its 
residues in all the point P of Vq which are the intersection of n + 1 divisors Di (this 
set could be void)§. 

The last step is to connect the Poincare residue and the monodromy filtration on 
the space of solutions of this differential module. This is more or less contained in [143]. 
Actually this paper shows how to compute subspaces of differential with logarithmic 
poles of given order by means of spaces built from the monodromy filtration. A by- 
product of this construction (cf theorem 12 in [49]) says that a differential of i/"(14) 
the Poincare residue of which is (i.e. it has a zero residue for each P) is in the kernel 
of any solution of maximal (monodromy) weight for Ly ■ 

As a consequence, the kernel of a solution of maximal weight for Ly contains the 
intersection of the kernel of the solution corresponding to integrate on the vanishing 
cycles Cp. But the monodromy filtration is characterised by its "dual" filtration on 
H"{Vx) given by corresponding kernels. So we conclude that the solution obtained by 
integrating on C is in the span of solutions obtained by integrating on the vanishing 
cycles Cp. Hence it takes its value in the set of diagonals of rational functions (in 
2n + 2 variables) . 

When Ly is MUM a simpler argument is the following : by hypothesis the solution 
associated to C is of maximal weight and the differential, we started with, is not in 
its kernel because / 7^ 0. So there is, at least, one point P such that integration 
on Cp of that differential is not zero and gives a diagonal of rational function g. In 
particular g is analytic (near zero) and one can conclude by unicity, up to a constant, 
of the analytic solution for Ly- 

Appendix F. Other hypergeometric "blind spots" for Christol's conjecture 

Let us give a list of 3^2 that are not algebraic hypergeometric functions]:, that are not 
obviously Hadamard product of algebraic functions, but actually correspond to series 
with integer coefficients: 



and then are in H"{Vx)- The naming logarithmic pole comes from this remark. It is much less 
obvious to prove that the Poincare residue is actually well-defined on H"{Vx) but it does. 
I The 2F1 case is well-known. 




(F.l) 



§ It is rather easy to convince himself that differential ■ • ■ ■ ■ ■ ^ " 



are exact ones for h > 1 
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where the four integers (Ni, N2, N^; Mi) read respectively: 

[1,2,7;2], [1,2,8;2], [1,4,5;1], [1,4,7;1], [1,4,7;2], [1,4,8;2], [1,5,8;1], 

[1,7,8;1], [2,4,5; 1], [2,4,7;1], [2, 5, 7; 2], [2,5,8;1], [2, 5, 8; 2], [2,7,8;1], 
[4, 5, 7; 2], [4, 5, 8; 2]. 

The series expansion of the first eandidate reads: 

3J^2([^, ^, ^,], 1], a^x) = 1 +21 a; +5544x2 + 2194500^3 

+ 1032711750 x'' + 535163031270 + 294927297193620 

+ 169625328357359160 + 100668944872954458000 + • • • 
Other examples read for instance: 

3F,[[^, I ^,], [i, 1], 7^t), 3^^2([^, ^, ^,], 1], ll^x). (F.2) 

Do note that, even if these various hypergeometric functions look very much 
ahke, the linear differential operators that annihilate them are not equivalent (no 
homomorphisms between these operatorsf or their symmetric powers). 

One can also try to find, systematically, hypergeometric functions that are 
not algebraic hypergeometric, that are not obviously Hadamard product of algebraic 
functions, but actually correspond to series with integer coefficients. Note that some 
of these 4F3 are deduced from the previous 3F2, as a consequence of a Hadamard 
product by an algebraic function: 

= 3i^2([^, I ^,], [| 1], S^x) ^ (1 -23x)-i/^ = 1 + 42 X +55440x2 

+ 131670000x3 + 402757582500 x'' + 1419252358928040 x^ 

+ 5475030345102361680 x^ + 22492318540185824616000 x^ +••• 

Appendix G. Proof of integrality of series (77) 

Let us sketch the proof of the integrality of series (77), namely, the integrality of 
coefficients (79). For each power of the integer number q = a, term like 4 + 9n 
is periodically divisible (period p) by q. In order to have the ratio (80) be an integer, 
one needs the numerator to be divisible by this factor q before the denominator. The 
case p = 3 is an easy one. The other prime p do not divide 9. One needs to find 
the first case of divisibility, namely the first integer n such that 4 + 9 n = kq (this 
corresponds to the smallest k). If dq ~ 1, mod. 9 then k = Ad, mod. 9. In other 
words, the smallest k is the rest of 4d, mod. 9. Consequently, we have replaced the 
calculations, for every integer g, by a, finite set 0/ calculations for d = 1, 2, 4, 5, 7, 8. 
Let us use this approach for the ratio (80). 

Remark: The terms n + 1 are always the last to be divisible by q. Hence, one 
can forget its factors. However, one needs as many factors at the numerator than at 
the denominator. For the other terms, the following table of the rest of d • a gives the 

f Associated with those various hypergeometric functions with the same singularities at x = 0, 1, 00. 
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(G.l) 

One finds out that this is always a factor of the numerator, before the occurrence of 
a factor at the denominator. 



Appendix H. Integrality of differential geometry modular form series 

Appendix H.l. Golyshev and Stienstra examples [117] 

From a differential Geometry viewpoint, Golyshev and Stienstra gave a set of selected 
order-three linear differential operators in [117]. The Wronskians of all these Golyshev 
and Stienstra examples, displayed in [117], are square roots of simple rational 
functions. Consequently, the differential Galois groups of the order-three operators 
displayed in [117] will be 0(3, C) instead of S0{3, C). 

Furthermore, all these Golyshev and Stienstra order-three operators are 
symmetric squares of order-two linear differential operators. 

For instance for G5, it is the symmetric square of 

■W = d2 + l- 66x- 32x^ _ 3-ix + l) 

^ ^ (1-44.T -16a;2).a; (1 - 44 x - 16 x2) • x ' ^ 

This Hcunian second order operator "Hs is actually solvable in 2F1 hypergeometric 
function with a (modular) puUback 

Fi(ii;^^]Ai];Pu) (H.2) 



K 



1/4 



'12' 12' 

To see this, one can, for example, calculate the nome q in terms of the variable x from 
the series solutions yo{x) and yi{x) in x of (H.l), i.e. from q = exp{yi{x) / yo{x)) , then 
obtain a series expansion in x for using the well-known expansion of the modular 

invariant j{q) = 1728/Pu = l/q + 784 + 196884g -I- 21493760^2 -\ in terms of q, 

and then recognise the algebraic equation satisfied by using for example Maple's 
command gfun[seriestoalgeq]. One obtains: 

(144 -f 216 X + 1)3 . p2 _ 1728 x (3456 x^ + 7776 x^ - 12600 x^ 

+ 1890x2 -80 a; + 1) • P„ + 2985984 a;^ = 0, (H.3) 



then: 



864 X ■ (1 - 80 X + 1890 x^ - 12600 x^ + 7776 x^ + 3456 x^) 
" ^ (1 + 216 x-f 144x2)3 ^ ■ 



^^^-■;^'^-ni-18x)(l-36x) , 
(1 + 216X+ 144x2)3 ^ ' ' 



5- (36x- 13) + 60- (1 - 44x - 16x2) 
1 + 216 X + 144x2 



2U/2 
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They correspond to the following series expansions: 

Pu = 1728 a; - 1257984 • a;^ + 575828352 • a;^ - 214274336256 • a:^ 

+ 70880897026368 • a;'^ - 21731780729723904 • + ■■■ , 
K = -125 + 28500 a; - 6123000^2 + 1318794000 -283968657000 a;'' 

+ 61 147607046000 x^ - 13166982207738000 a;^ + ••• 

Note that the series for k, and Pu, have integer coefficients. 

Introducing the rational parametrisation of curve — (1 —44 a; — 16a;^) = 0, 
in order to get rid of the square root (1 — 44a; — 16a;^)^/^. namely 

M ^ , (m' - 125) 

(125 + 22m +m')' ^ (125 + 22/1+^2)' ^ ■> 

the corresponding two puUbacks and k's reading 

1728/1 5 (125 + 22 /i +/i2) 



Pu = 
P, 



(5+10m+m')'' (5+10/i+/i2) ' 

1728 /i^ 125 (125 + 22 /i + 



(3125 + 220 m +Ai^)3' (3125 + 250 M +Ai^) 

where it is straightforward to see the "Atkin" symmetry fi i — > 125/ fi: 

.( V K...(lt^ = K...( V 



The relation between the two pullbacks Pu and Py corresponds to a (rational) 
modular curve. One immediately recognises the modular curve corresponding to the 
elimination of fj, between (the two Hauptmoduls) P„, and Py, which is well-known 
to correspond to q — >■ g^, or r ^ 5 • r (namely the fundamental modular curve 
*(j(t), j(5-r)) = 0). ^ 

All these results for equation (H.l) can be summarised in the following equation 

Pbio-b ■ a;)^^ • Fullback 7^5,3; — >■ 05(05 • a;) • p^ia^ ■ x) = ^5, (H-''') 

where we denote by uj^ , the order- two operator corresponding to the modular solution 
-D5 • 2F1 {[y2^ I-*^]' ' ''^'icre I?5 and 75 arc given by Maier (see tables 4 and 12 

in [123]). The expression of C5(a;) is given in (G.4): c^{x) = x/{x'^ + 22a; + 125) 
and P5(a;) is given below. 

This means that the order-two linear differential operator, associated with 
Golyshev and Stienstra example G5, when pull-backed by appropriate functions 05(3;), 
and P5(a;), is simply the operator corresponding to the modular solution j^. 

In a similar way, one can perform the same calculations for other examples given 
in [144]. The results are displayed in the following table with the same notations as 
before: 

p„(a„ • a;)"^ • Fullback T-LmX ^ Cn{an ■ x) -pnian ■ x) = a;„, (H.8) 
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Cji (x) 


Pri(2;) 


n — A ■ R 


2 


x/{x + 64)2 


(a; + 64)^/'' 


212 


3 




(a; + 27)i/a 


36 


4 


+ 16)2 


(a; + 16)1/2 


28 


5 


a;/(x2 + 22x+ 125) 


(x2 +22a;+ 125)1/"^ 


2253 


6 


x/{x + 9)/{x + 8) 


((a; + 9)(a; + 8))i/2 


23 32 


7 


xl{x^ + 13X + 49) 


(x2 + 13a; + 49)1/3 


72 32 


8 


a;/(a; + 8)/(a- + 4) 


((x + 8)(:e + 4))1/2 


25 


9 


x/(a;2 + 92: + 27) 


(a-2 +9:r + 27)1/2 


33 



Operators "Hn for n = 6, 7, 8, 9 are given in the next subsection (for n = 2, 3, 4, 
see [144]). 

Let us note that aU these examples are associated to genus zero curves. Another 
example Gn given in [144] which is associated to a genus one curve will be considered 
in detail in Appendix I. 



Appendix H.2. More details 

The order-three operator Gg is the symmetric square of 



X- (l-34a; +a:2) 4a; -(1-342; 

Let us introduce the rational parametrisation of — [1 — "iA x + x"^ ) ~ 0: 
/i /i2 - 72 



(a. + 9)-(a* + 8)' ^ (^ + 9).(/, + 8)- 

With this new parametrisation the operator (H.9) becomes 

Lf, = 4 • (^j + 9)2 . (/I + 8)2 . 6i2 -10-(M + 9)-(/i + 8)-Ai + M^ (H.IO) 

which is covariant by the "Atkin" involution /i -;->■ 72//x, a; being invariant by this 
involution: 



x(^) = a;(M), = -y(M), (H.ll) 

It is worth recalling the rational parametrisation of the modular curve r — > 6 • r 
namely [123]: 

. _ (m +6)3. (/x3 +18^/2+ 84 M + 24) 

" /.-(a^ +8)3.(^+9)2 ' ^"^-'^^ 



, {fi + 12)3 . _^ 252 Aj2 + 3888 a* + 15552) . / 72 



^3 . (// + 252 Aj2 

m'^-Ia^ +8)2. ifi + 9) 
The solutions of (H.IO) read 

'{^l+8)^^l +9)Si/i2 /I 12. je 



^^,f[±,±],[^];_Z^), (H.13) 

VW 12^' ^3^' 1728/ ^ / 



but can also be written as: 

(/i +8)2 (^ +9)2 xi/4 / I 5, 1728 



(Ai + 6)(/i3 +18^2^ 34^ ^ 24)/ 'V42'l2''''' je / 

36 . (x + 8)2 (a; + 9)2 xV* 



.(Ai + 12) . (/x3 + 252 //2 + 3888 /i + 15552)/ ^ " V^12' 12^' ^ ^' j'^ 
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The third order operator 6*7 is the symmetric square of the second order hnear 
differential operator 

^ (l-27x)(a; + l)a; (1 - 27 a;) (x + 1) a; ^ ' 

Introducing the parametrisation of the rational curve 

- (I +a;) • (1 - 27x) = 0, (H.16) 

namely 

/X /i^ - 49 



where one verifies the existence of an "Atkin" involution: 



x[j) = y{j) = -y{til (H.I8) 



With this change of variables the second order differential operator (H.f5) reads: 

= (/i^ + 13 /i + 49)2 • 0^ - 2 • (/i^ + 13/i +49) • /i (H.19) 

It is worth recalling the rational parametrisation of the modular curve r — > 7 • t 
namely [123]: 

+i3m +49)-(a^2 + b^l + 1) 

J7 = , (M.^U) 

, (/i^ + 13 At + 49) • ^ 245 /i + 2401) /49 
y7 \ ^1 

The solution of (11.19) reads: 

'0? +13Ai +49)Si/i2 /I 12. j7 



V /i J ' ''^VW 12^' ^3^' 17287' 



'7*^ • (/i^ + 13^ +49)^1/12 / I 12 j4 



\ y V12' 12" '3" 1728 

The order-three operator Gs is the symmetric square of 

Us ^ Dl + 1-36. + 32.2 -2 , -. (H.21) 

^ X • (1 - 24a; + 16a;2) a; • (1 - 24 a; + 16 a;2) ^ ' 

Furthermore one has an "Atkin" symmetry x O 1/16/ x. Introducing the 
parametrisation of _ (1 _ 24 a; + 16 a;^) = 0, namely 

11 I? -32 , ^ 



(a. +4).(a. +8)' (a* +4).(/i +8)' 

the linear differential operator (11.21) becomes 

= (a^ + 4)2 . (a. + 8)2 . 02 _2.(a. +4).(m +8)-M +4-^', 
this operator being covariant by the "Atkin" involution which leaves x invariant: 

x(^) = .(/.), y(^) = -2/(a.). (H.23) 



Diagonals of rational functions 



63 



It is worth recalling the rational paramctrisation of the modular curve 
representing t — > 8 • r, namely [123]: 

(// + 16 ^3 + 80 fi^ +128^1 + 16)3 
• (/i + 4)2 + : 



J8 = .. ... , .^2^. , ON ' (H-24) 



., _ (// + 256 Ai'^ +5120Ai-^ + 32768 A* +65536)-^ _ . /32 



I a ) ■ ^^^l[l2' 12]' [3]' 1728 J' (^-'^^ 



A solution of (H.23) reads 

'(/z + 8)_^M + 4)4^1/12 ^ 1 1 ^ ^2^^ is 
M 

but can also be written 

((M+4).(,+8))^^^.F,([i,i],[l];-^i:l^), 
or equivalently, using Gauss-Kummer identity: 

Finally, the order-three operator Gg is the symmetric square of 



^ X- (1 - 18x - 27x2) 4a; • (1 - 18x - 27x2^ ^ ' 



With the paramctrisation of the rational curve y2 — (1 — 18 a; — 27a;2) = 0, 
namely 

~ ~2 m I nr. I 2/ = ~2 l~n , r,-7 ' (H.27) 

/i^ + 9 /i + 27 /i^ + 9 Ai + 27 

the linear differential operator (H.26) becomes 

= +9fi +27f -el -6-{fi^ +9fi +27)- -27 

which is covariant by the "Atkin" involution leaving x invariant: 

x(^) ^ a;(M), y(^) = -Viti). (H.28) 

It is worth recalling the rational paramctrisation of the modular curve r — > 9 • r 
namely [123]: 

Qi +if■{^i^ + 9^?+2^^x + 3) 

J9 = / 2 I n n^r^ ' (H.29) 

Ai • (/.t"^ + 9 Ai + 27) 



(/i + 9)3 • (Ai3 + 243 Af2 + 2187 + 6561) / 27\ 

= ..g ...2 , n.. , o^x = -^^^TJ- 



+911 +27) ^ ■"'^11 

A solution of (H.28) reads 



(a.!±iii±^I)!)-..,(,>,^u|,:^), 

but can also be written as: 

((m' + V +27))"^ . ,F.{ii.ll [11; - "■"■';;>■+"' ). (H.31) 
/27- (/.t2 + 9At +27)\i/2 /I 1 „2^9^_^27n 



.F,([-,-],[l];-27^ 



^^2 y ^ ^V^3' 3^' ' " ■ ^3 



The relation between these two puUbacks corresponds to the modular curve 
x^y^ - 270 • x^y^ + 972 ■ xy ■ {x + y) - 729 ■ {x^ + xy + y^) = 0. (H.32) 
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Appendix I. Higher genus modular forms 

For non-zero genus modular curves, we have generalisations of these structures 
associated with an "Atkin" involution of the form z A/z, which correspond 
to the introduction of the so-called Atkin^ modular polynomials, or star-modular 
polynomials [146]. 

Appendix I.l. A genus-one curve 

Let us first consider the order-three Golyshev and Stienstra linear differential operator 
Gil given in [144] which is the symmetric square of 

_ 2 (625 - 12750 X - 30800 ~ 3150 - 4512 x^) 
" ~ ^ + X- (125 - 1900x -40x2 -188x3)- (8a; -f 5) ' " 
6 • (125 700 X 105 a;2 + 188 a;3) 
~ (125 - 1900x - 40x2 _ 188x3) • (8x-l-5)' *■ ' -* 

Performing the same calculations as for (H.l), one deduces that the solution 
of (I.l) can be written in terms of pull-backed hypergeometric function, namely 
2^1 (ll^'^]' [1]' '^ith Pu — n2S/J, with J satisfying the following algebraic 

equation: 

244140625x^2- 7^- 5- (132645814272 x" + 372815032320 x^° -f- 1405869696000 x^ 
-(- 5229172080000 .x^ - 225383400000 x^ - 5599578600000x6 4- 339591656250 x^ 
+ 1103850000000.x'* - 349421875000 x^ 42861328125 x^ - 2363281250 x 
-f- 48828125) ■ x- J + (78336 x'' -I- 181440 x3 + 561600 x2 + 144000 x -I- 625)^ = 0. 

Then, when we perform the puUback x — >■ 5/(5 z — 3), we obtain the modular 
polynomial of order-eleven givenf by Morain [146], or by Elkies (see eq. (49) in [147]), 
associated to a genus-one modular curve: 

J2 -Qi(z)- J +(z'' + 228^3 -f- 486^2 - 540 z 4- 225)3 = 0. (1.2) 

where: 

Qi{z) = z^^ - 55 z^" -f 1188 - 12716 +69630 z'^ - 177408^6 

+ 133056 z^ 132066 z'' - 187407 z^ 40095 z2 24300 z - 6750. 
The j-invariant of the genus-one (J, z)-curve (1.2) reads^: 

122023936 _ 496^ 

^™ " 161051 " "TT^- ^ ' 

Appendix 1. 1.1. Atkin modular polynomial and modular forms for r — > llr 
Let us focus on the previous Atkin modular polynomial (1.2): 

^l^{z,j) = f -Qi{z)-j + (z^-h 228z3-h486z2 - 540z-f 225)3. (L4) 

Note that there is no associated 'Atkin" involution of the (rational) form z — > A/ z, 
since the curve $ii(z, j) = is a genus-one curve. 

f Atkin's work wc are interested in, cannot be found in papers but in emails [145]. 

f See <I>*j^(_F, J) in the last equation of subsection 2.3.2 of [146]. The variable F in [146] is z here. 

% Use the command algcurvesp Jnvariant] in Maple. 
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The elimination of z, between the two solutions and j_ of ^li{z, j) = 0, 
yields the modular curve 

<I>n(j+,J-) = 0, (1.5) 

which is a quite large (146 monomials) j_ )-symmetric polynomial of degree twelve 
in j+ (resp. j-). The discriminant of $ii(j+, j-) in j_ reads (we denote here j- 
by j): 

- 11" ■ / • (j - (12)^)^ • • (1.6) 

where Q reads: 

Q = U + (15)^) • (j - (20)3) . {j + (96)3) . (ggo)3) . _ (255)3) 

X {f - 425692800.7 +9103145472000)2- {f + 117964800 j - 134217728000)^, 

and where P is a polynomial that factors into the product of seven polynomials of 
degree two, fourteen polynomials of degree four, one polynomial of degree six and four 
polynomials of degree eight: 

r = pi^' • • • • pi'^ ■ ■ ■ pi''^ ■ Fi^^ • p^^ ■ ■ ■ . (1.7) 

In the j+ = j- = j limit it becomes a polynomial of degree 22, which simply 
factors as follows: 

$ii(j, j) = -(j +(32)3). q2 

X if - 1122662608 f + 270413882112 j - 653249011576832). 

where, besides the Complex Multiplication value [148] j = (255)^ = 16581375, one 
recognises a large set of Heegner numbers [17], corresponding to the following integer 
values of the j-invariant: (20)^, -(15)3, -(32)3, _(96)3, -(960)3. Note that the 
j-invariant of the genus-one (j^, j_) modular curve also reads the same j-invariant 
as the one for the genus-one curve (1.2), namely —122023936/161051 (see (1.3)). 
Recalling the well-known expansion of the j-function as a function of q: 

j{q) = ~ +"^4"^ + 196884 g + 21493760 +864299970^3 (jg) 

+ 20245856256 + 333202640600 + 4252023300096 q^ + ■■■ 
one verifies immediately that it does correspond to t — > 11 r: 

^n{jiq), jiq'')) - 0. (1.9) 

Recalling the modular discriminant of Weierstrass, A{q) (i.e. the 24-th power of the 
Dedekind eta function up to (27r)^2 factor), it reads: 



A(g) = q 11(1 



q-r- 



Let us introduce the Eisenstein series^ £^2(9) and the Euler-likc function p{q): 

q dHGjq)) A(<7ll) 
E2{q) = XT 1 , where: G{q) = 



24 dq ' ■ A{q) ' 

p{q) = (1-9")' -11(1 (I-IO) 

71—1 n— 1 

X Quasi-modular form [149]: G2((ar +6)/(cr + d)) = [cr + d)^ ■ G2(r) - {cr + d) c/A/Tr/i. 
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and, finally z{q) (denoted F in Morain [146]): 

3 12 £'2(0) 1 9 A 

z{q) = 7+^ 7^ = - +5 +17g +46g2 +116g^ + 252g-* 

5 5 p{q) q 

+ 533g^ + 1034g*^ + 1961 + ••• (I.ll) 
One verifies that the two relations on the Atkin modular polynomial (1.4) are satisfied: 
$ti(z(g), j(<z)) - 0, and: $ti(z((7), jCz")) = 0. (1.12) 

Introducing x{q) — l/z{q) one has the following series expansion with integer 
coefficients: 

x{q) = q -5g2 +8g3 -g** -17g^ +62^^ -176g^ + 339g** - 386g^ + ••• 
Let us now show an identity characteristic of modular forms: 

where A(x) is the algebraic function 

A{x)^ 112 2 7321 - 87612x4- 73206x2 + 21060x3 - 23175x4 



112 ^(2:)4 112 1 + 228 x + 486x2 - 540x3 + 225x4 

Appendix 1.1.2. A change of variables 

Let us rewrite (1.4) in this variable x = 1/z, and in the Hauptmodul H = 1728/j, 
introducing a new star-modular polynomial Pii(x, H) = ■ 2^2. $*^(l/x, 1728/i7): 

Pi\ {x, H) = (1 + 228 X + 486 x2 - 540 x^ + 225 x*)^ • (L14) 
- 1728- Qi(x) -x • + 17282x^2^ ^j^h: 

Qi(x) = 1 - 55x + 1188x2 - 12716x3 +69630x'* - 177408 x^ + 133056 x^ 
+ 132066 .x^ - 187407 x^ + 40095 x^ + 24300 x^" - 6750x". 
The two Hauptmodul solutions of polynomial P^^{x, H) expand respectively as: 
1 H, 



739 x^ +349254x3 - 135092042 x* + 46600204623 x'^ 



j{q) 1728 
1 u 

Jiq^ ^ TtIs " +55X^2 _^ 2^837 2,13 _,_484;L;L ^14 ^ j^j^Qgggg ^15 ^J j^g-, 

+ 23244727x^6 + ... 

The corresponding expansions of 2i^i([l/12, 5/12], [1]; Hi) and 2-Fi([l/12, 5/12], [1]; H2) 
read respectively two series with integer coefficients: 

2Fi([l/12, 5/12], [1]; iJi) = 1 + 60 X - 4560x2 +614400x3 - 95660400 x'' 

+ 16231863060 x^ - 2905028387700 x*^ + • • • + 20000242239261022140 x^ 
- 3953288123422938241560 x^" + 791518845663517087144740 x" + ••• 
2Fl([l/12, 5/12], [1]; iJa) = 1 + 60 x" +3300x^2 _^ 110220 x" +2904660 x" 
+ 66599940x^5 + 1394683620x^6 + 27425371380 x^^ + ••• (L16) 
Defining A{x) as the ratio of these two series expansions: 

2^l([^, ^], = A(x).2Fi([l, A], [1];H2), (L17) 
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one can easily see that this ratio A{x) is solution of the genus-one algebraic equation: 

112 + ^(,j.)4 ^(2.) ■ y-°> 

It is worth noting that this genus-one algebraic curve (1.18), in the two variables 
X and y = has the same j -invariant, namely —122023936/161051 (see (1.3)), 

as the genus-one (x, j/)-curve (1.2) corresponding to the Atkin-modular polynomial. 

The ratio A{x) then reads: 

Aixf = ^,ere: (1.19) 

5{x) 

a{x) = 45^2 - 246a; +61, P{x) = (1 + x) • (1 - 17.t + 19.t2 - 7x^), 

j{x) = 7321 - 87612 a; + 73206 a;^ +21060 2;^ - 23175 a;^ 
6{x) = 1 +228 a; +486^^ - 540a;3 + 225x''. (1.20) 
Its series expansion is a series with integer coefficients: 
A{x) = 1 + 60 a; - 4560 x^ + 614400 a;^ - 95660400 x^ + 16231863060 x^ 

- 2905028387700 + ••• + 20000242239261022140 a;^ (1.21) 

- 3953288123422938241560^1" + 791518845663517087144680 a;" + ••• 

Introducing the polynomial C{x) = 1 + 8 .t — 9 a;^ + 10 .t^ — 6 a;'', one has the 
following relations on A{x): 

dlnA{x) 60 -({x) 



dx /3{xy/^ ■ 5{x)' 

Appendix 1.1.3. Order-four operator 



(1.22) 



Let us introduce the order-four linear differential operator L4 annihilating 
2-Fi^[l/12, 5/12], [1]; Hi). Quite remarkably this linear differential operator L4 also 

annihilates 2Fi{\^/12, 5/12], [1]; H2). Therefore, this order-four operator is not a 

MUM operator (it has two series-solutions (1.16), analytic at a: = 0). 

The symmetric square of L4 is of order nine, instead of the order ten one can 
expect generically. Furthermore, the exterior square of L4 is of order six, but, 
remarkably, this exterior square factorizes into a direct sum of an order-one operator. 
Ml, an order-two operator, m!^\ and an order-three operator, which is the symmetric 
square of another order-two operator, M^^: 

ext^{Li) = Ml ® Af^i^ ® Sym^{MP), (1.23) 

where 

^ 1 dhi(ii(x)) . , , , S(x) 

D.--, — wrth: ,(-) = (124) 

the Wronskian of the two order-two operators M^^ and M^'' reading respectively: 

- 4 i^f ,r V W^KMf')^ = 4 ,3 - (1-25) 

x^ ■ pyxy ■ o{x) x^ ■ p(x) ■ o(xy 

This factorisation is a consequence of relation (1.17) and indicates that L4 is very 
special. 
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The Hauptmoduls Hi or H2 arc not rational functions of x, they are algebraic. 
Therefore 2^1 (^[1/12, 5/12], [1]; Hi) has no reason to be solution of a second-order 
operator. As far as factorisation of hnear differential operators in operators with 
polynomial coefficients, L4 is irreduciblef. Let us show that such a reduction to a 
second order linear differential operator actually exists. 

Appendix 1.1.4- ■ Order-two operator: operator ujh 
Introducing the algebraic function A{x) 
Aix) 



5(x)i/8 . A(x)i/2 - 120 • a{x) ■ /3(a;)V2 

• (5(a;)3/8 /-I d (lM(a;)) \ 



C(X)V2 

which has the series expansion 

117 64635 2 21853425 . 32050683795 4 
= 1 +^-' + 128 ^ 
12299248285371 9718868161850799 . 
256 ^ + K)24 ^ (^-2^) 

one finds that 

^ ^ 1 5 . ^ \ ^ 3 75 2 1335 3 

-^(^)-^^([t2'i2]'W'^0 = 1 +2-" 

111795 4 2559789 . 124177119 . 

is actually solution of an order-two linear differential operator: 
_ 2 l-24;r + 4.x2 + 30x3-21a;'* 

= ^^'^ + a;(a; + l)(l -17x +19x2 - 7 x^) ' 

2 



3 ix~l){7x^-x-2) 



(1.29) 



Ax - (x + l)(l -17a; +19a;2 -7x3)' 
Note that 

Mi'^ = ^ • ^11 • A(.), with: X{x) = ^ ■ (1-30) 

The series (1.28) is globally bounded, the rescaling x ^ Ax changing this series 
into a series with integer coefficients, 

1 +6 a; + 150 +5340 a;^ + 223590 x'' + 10239156 + 496708476 x^ 

+ 25083657720 x^ + 1304819854470 x^ +■■■, (1-31) 

solution of the order-two operator (puUback of (1.29) by x — 4x ): 

^2 1 - 96 X + 64 x2 + 1920x3- 5376 x'' ^ 

~ ^ X • (l + 4x) • (1 - 68 x + 304x2 -448x3) ' 



6(l-4x)(l + 2x-56x2) 



X • (l + 4x) • (1 - 68X + 304x2 - 448x3)' 

t The command DEtools[DFactor] in Maple. 



(1.32) 
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Appendix 1.1.5. From order-two operator to order-four operator 



This result can be revisited as follows. The function (^[1/12, 5/12], [1]; i/i), 
known to be solution of the order-four operator L4 is also solution of the order-two 
operator: 

where: = S{x)^^^ ■ wn • ^^^^^^g (1.33) 



X ■ S{x) ■ f3{x) X ■ S{x)^ ■ fi(x) 

^2 _ 1 d\TL{}j.{x)) ^ ^ ^ po{x) 



4 dx X ■ J(a;)2 • /3(x) ' 

with fi{x) is given in (1.24) and 

(a;) = 1-1- 147 a; - 4313 - 7083 x^ + 14073 x^ + 4125 x^ - 19395 a:^ 

+ 12555 a;^ - 3150 a;®, (1.34) 
Pq{x) = 1 -19 a; -M389a;^ -6497 a;^ 27603 a:'' - 37155 a;^- 18369 

+ 45477 a;^ - 17280 x^ - 270 (1.35) 
and where the order-one operator 51i reads (using (1.22)): 

dln^(a;) / pi(x) \ d'^\nA(x) 

60 • /3(a;)V2 . C(a;) ^o^^ ^ rfln(^(.T)) x <5(a;)5 



x) • 5(x) V ^ 4 dx J' ■ ' ' a;4 • C(a;) 



The function 2^'i(^[l/12, 5/12], [1]; i/2), known to be solution of the order-four 
operator L4 is also solution of the order-two operator: 



1 



wii • A{x) ■ A{x) = \ ■ wii • A{x)^/^ 



A{x) ■ Aix) - ^ ^ ^ ^ A{xy/ 



1 /dln^(a;)\2 r^i 



4 V / 2 

These last two order-two operators (1.33) and (1.37) can be written as iln ±ili/2 
where: 

1 fd\nA{x)\2 



fill — '^ii 



4 V da; 



X • (5(x) • /3(x) X • (5(x)2 • /3(x) ' 

with 

qoix) = 1 -49x 909x2 - 7877x3 + 31323x4 -44025 x^ - 10089 x*^ 

+ 39237 x^ - 13680 x® - 1350 x^. (1.38) 

These last two order-two operators (1.33) and (1.37) are not linear differential operators 
with rational coefficients, but with algebraic coefficients: there are /3(x)^/2 terms. 



Diagonals of rational functions 70 

One can verify directly that the order-four operator L4 can actually be seen as 
the direct sum of these last two order-two operators with algebraic coefficients (1.33) 
and (1.37): 

- (^(-)^'' • • A^) ® ( IRv^ ■ ■ ^(-)''') 

= (f^n -^) ® (f^ii +^). (1.39) 

Note that this result is a particular case of a more general result. Let us consider the 
direct-sum of the two order-two operators with algebraic coefficients depending on one 
parameter u 

U{u) = {nu +u - ^) (B {nn - u- (1-40) 

The order-four operator is actually a linear differential operator with rational 
coefficients for any u. It is of the form 

i4(tt) = Ma-l^ii +u^-Mi-Ni, (1.41) 

where M2 and Mi are respectively order-two and order-one operators with rational 
coefficients and the order-one operator A^i reads: 

1 dln(p(x)) . , , , (5(a;)5 

Ni = D.^-. w:th: p(.) = (1.42) 

We have the following general result: the exterior square of a direct sum of the 
form (1.39) is a direct sum of three order-one operators and the symmetric square of an 
order- two linear differential operator. Let us denote , W~ and W the Wronskian 
of respectively 

One has the following direct sum decomposition for an arbitrary order-two linear 
differential operator ujh: 



ext'^[Ld^^ ® w^ij = (L43) 



(3) 

where the order-three operator M2 can, for instance, be written as a conjugation of 
a symmetric square (see also (L30)): 

.(3) _ ..,,^2fdA{^'^',-fdAix)Y^ 



M, 



2 



The order-one operator — ^'^^Y^ in (L43) actually corresponds to (L33) together 
with (L24). However we see from (1.43) that the irreducible order-two operator 
we found in the direct sum decomposition (L23), can in fact be decomposed in a direct 
sum of order-one operators with algebraic coefficients. 
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Appendix 1.1.6. Back to Golyshev and Stienstra operator 

Do note that the order-two operator "Hn, previously encountered with a genus-one 
situation (see (I.l)), is conjugated to a pullback of (1.29): 

= ■ ^ rrV) • \/5 + 3^. (1.45) 

V5 + 3x V 5 + da;/ 

The puhback ofHuby >5x reads: 

n 4512 + 630x3 + 1232x2 + 102 a; -1 



Hii^x ^ 5a;) = Dl 



(188a;3 + 8a-2 + 76a;-l)- (Sx + l) 
6(188a;3 + 21x2 + 28a;+l) 



(I- 



(188x-3 + 82;2 + 76x-l)- (8x + l)- x' 
which has as a solution the series with integer coefficients: 

1 +6 a; + 204 a;2 + 8790 + 445170 x"* + 24577236 a;^ + 1436107596 a;^ 

+ 87310665684 x'' + 5466252149820 a;* + • • • (1.47) 

Appendix 1.1.7. Hadamard products 

The Hadamard square of wn is a linear differential operator of order ten. The 
Hadamard square of "Hn (or its pullback by x 5 a;) is also a linear differential 
operator of order ten. This order-ten operator is not MUM. 
Its head polynomial is 

a;^ • (1 -64 a;)- (96256 a;^ - 512 a;^ + 608 x + 1) • (35344 a;^ + 28512 a;^ + 5792 a; - 1) 

X (35344a;3- 14288a'2-8a;-l)- P29, 

where P29 is a polynomial of degree 29. One verifies easily that Hadamard's 
theorem [62] on the location of the singularities is verified. If one denotes xi, X2, 
X3, the three roots of polynomial 188a;'^ + 8a;2 + 76a;— 1 (see (1.46)), one sees that the 
three roots of polynomial 96256 a;^ — 512 a;^ + 608 x + 1 are nothing but — a;i/8, —X2/S, 
—X3/8, the three roots of polynomial 35344 x^ + 28512 a;^ + 5792 x—1 are nothing but 
a;2, X2, x\, the three roots of polynomial 35344 a;"^ — 14288^2 — 8x — 1 are nothing 
but xia;2, a;2 xa, xi X3, and of course, besides x = and x — 00, 1/64 = (—1/8)^. 



Appendix 1.2. Other higher genus modular forms 

Similar calculations can be performed for the modular forms associated with genus- 
one modular curves, corresponding, for instance to r — > A^r, for N = 17, 19. For 
= 23 the modular curve is a genus-two curve [150, 151]. These detailed analysis 
and calculations will be given in a forthcoming publication. 

The A'^ = 23 genus-two case requires to introduce the order-two operator \: 
2 1 - 12x +4x2 +5x3 _ 33^4 +353,5 _ 283-6 
^^23(2;) - D^, + ^. (1 _8x +3x2 _7^3). (1 _3.2 +3,3) ■ ^- 

1 - X - x2 + 12 x^ - 15 x-* + 14 x^ 



X- (1 -8X +3x2 „7 3;3). (1 -j;2 +3,3)' 

This operator has the analytic solution with integer coefficients: 
1 +x +3x2 +13x3 +67 X* +375 x^ +2223 x*^ +13713 x^ + 87123 x^ 

I where one notes that the polynomial 1 — + has the root — ^/ P, where P = 1.324717958 ... 
is the smallest Pisot number [152]. 
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The pullback of 0)23 (x) by x — > l/x gives an order- two operator with two analytie 
sohitions (no logarithm): 

5 2 11 . 85 4 3499 , 2041 . 18317 7 



14 196 1372 57624 57624 672280 

332455 o 21994361 „ 

x^ H x^ + ■■■ , (1-48) 

19765032 1383552240 ^ ' 



and 



„2 



5 , 3 4 251 . 137 507 ^ 24007 o 

x' H a;-^ X* x'' x^ x** 

14 98 1372 1372 9604 470596 

144083 



(1.49) 



+ ... 

6588344 

These two series are not globally bounded (but a hnear combination of these two 
series may be globally bounded ...). The fact that the two previous solutions have no 
logarithmic terms excludes any relation like (149). We encounter the same situation 
with the order-two operators J)ii(x), cJi7(x), wig(x), W29(x), W3i(.x), W4i(x), (D47(x), 
a;59(x), and d'7i(x), the corresponding two series at x = cxj having no logarithmic 
terms yielding the same obstruction for a relation like (149). These various order-two 
operators correspond to higher order genus modular curves [153], namely genus-one 
for uJii{x), uJir(x), ujiq(x), genus-two for W29(x), (I'3i(x), genus-three for a)4i(x), 
genus-four for a;47(x), genus- five for ujr^g^x), and genus-six for (I'7i(x). 

Note that all these higher-genus a)„(x)'s are simply homomorphic to their adjoint. 
They are such that 

i?(x)i/^^ • adjoint{Cj„{x)) = w„(x) • i?(x)i/^^ (1.50) 

where R{x) is a rational function, and where M ~ 2, except for n = 19 where 
M = 6. 

Appendix 1.3. Order-two operators associated with modular Atkin equations 

Let us give the explicit expressions of some cj„'sfor n = 17, 19, 23, 31, 39, 41, 47, 59, 71. 
The higher genus a;„'s read 

with: 



All = 


84: x^- 120x3- 


16x2 -f 


96x-4, 




Bn = 


3(x- 1) {7x^ - 


x-2). 






Cu = 


4x • (7x3 - 19 X 


2 + 17x 


-l)(x + l). 




An = 


448 x^ - 576 x'^ - 


}- 280x3 


-1-224x2 - 264x- 


f 16, 


Bn ~ 


168 x^ - 180x3 - 


f 67x2h 


-41x- 20, 




Cn = 


16x(x- l)(8x4 


-4x3-1- 


3x2-|-10x-l), 




Ai9 = 


264 x^ - 354 x'^ - 


- 258 x3 


+ 300x2 + 93 X- 


9, 


Big = 


110 x''- 130x3- 


~ 69x2 ^ 


- 56 X + 6, 




Cig = 


9x • (2x + 1) (x 


-fl)(4x3- 12x2-f 10x- 


-1), 
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A23 ^ 28 - 35 + 33 x'' - 5 - 4 + 12 X - 1, 
B23 = 14 x^ - 15 x'^ + 12 x^ - .T^ - X + 1, 
C23 = X- (7x3-3x2 + 8x-l)(x3-x2 + l), 

A29 = 504x^ -64x^- 896x^-480x^ + 400x3 + 512x2 + 72x- 
B29 = 315 x^ - 35 x^ - 441 x* ~ 213 x^ + 109 x^ + 103 x + 4, 
C29 = 16x- (x + l)(7x^ - 8x^-8x'* -2x3 + 12x2 + 4x-l), 

^31 = 108 x*^ + 343 x^ + 477 x^ + 235 x^ + 28 x^ - 6 x - 1, 
B31 = 60x^ + 161 x* + 180x3 + 69 x2 + 5x-l, 
C31 = X- (x3 +3x2 + 4x+ l)(27x3 + 17x2 - 1) 

Ail = 704 x^ + 960 x*^ - 360 x^- 1472 x*^- 672 x^ + 480 x"* 

+ 720 x3 - 128 x2 - 120 X + 16, 
Bii = 616 X® + 756 x^ - 289 x^ - 969 x'"^- 353x^ + 265x3 

+ 261x2 - 39x- 12, 
C41 = 16x- (x-l)(8x® + 20x^ + 15x^-8x^-20x'*-10x3 

+ 8x2 + 4x- 

A47 = 66xi°-154x^ + 190x^-135x^ + 52x^ + 56x^-57x'^ 

+ 60x3 - 22x2 + 9x-l, 
S47 = 66 x^ - 140 x*^ + 158 x^ - 103 x^ + 35 x^ + 33 x" 

- 28x3 + 24x2 - 6x + l, 

C47 = X- (Ilx5-6x^ + 15x3-5x2 + 5x-l)(x5-2x'' + x3 

+ x2 — X + 1), 

A5g = 308 xi2 _^ 5242.114. 1632 x^° + 2640 x^ + 3520 x^ + 3816 x^ 
+ 3232 + 2128 x'^ + 1008 x^ + 280 x3 - 24 x - 4, 

B59 = 385 x" + 720 xi° + 1748 x^ + 2600 x*^ + 3128 x^ + 3024 x^ 
+ 2217 x^ + 1216 X* + 448 x3 + 86 x2 - 4 X - 4, 

C59 = 4x- (x3 + 2x+l)(llx^ + 24x^ + 46x^ + 61x^ + 60x^ 
+ 44x'' + 21x3 + 4x2-2x- 1), 

Aji = 88 x^^ - 30 x" - 280 x^^ - 195 x" + 420 x^° + 671 x^ - 5 x' 

- 666 x^ - 444 x^ + 91 x^ + 231 x^ + 95 x3 + 4 x2 - 6 x - 1, 
B71 = 132 x^3 _ 42 2;i2 _ 372-^11- 248 x^° + 475 x^ + 717 x^ 

+ 25 x^ - 563 x^ - 362 x^ + 29 x'' + 112 x3 + 38 x2 + X - 1, 
C71 = X- (llx^ - 4x^ -18x^-5x'' + 11x3 + 7x2 -l)(x^-2x' 
-3x'* + x3 + 5x2+4x+l). 
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Consider an order-four MUM linear differential operator. Let us introduce the 
determinantal variables Wm = dct (Mm) which are the determinants! of the following 
mx m matrices Mm, m = 1, • • • , 4, with entries expressed in terms of derivatives of 
the four solutions yo{x), yi{x), 2/2(2;) and y3{x) of the MUM linear differential operator 
(see Section (5) for the definitions). One takes Wi{x) ~ yo{x) and: 



Mo 



yo 
y'o 



yi 
y'l 



where: 
Since q 



M. = 



dx ' 





yi 


yi 


y'. 


y'l 


y'2 










y'l 


y'i . 








yt 







Ma 



yo 


yi 


2/2 


y3 


y'o 


y'l 


y'2 


y's 


y'i 


y'i 


y'i 


y'i 


y'^' 


y'l' 


y'i' 


y'i' 



dx^ 



yi- 



exp{yi/yo), and hence, 



dq 



wl 

W2 



d 

dx 



yl 
W2 



d_ 

dx ' 



and thus 



dq 



dx"^ 



yl dyo 
dx 



d_ 

dx 



yt dW2 



wi 



dx 



d 

dx ' 



(J.l) 



(J.2) 



(J.3) 



we deduce, after some simple algebra, an alternative definition for the Yukawa coupling: 



K{q) 



q ■ 



dq 



(J.4) 



to be compared with the other previous alternative expression previously given (170) 
for the Yukawa coupling 



K(q) 



x{qf ■ 



1/2 



q ^ dx{q) Y 

x{q) dq ) 



1/2 



yl 



dq ) ■ 



In fact from (J.2) we deduce 



dq I 



Wi 



yl_ 



(J.5) 



(J.6) 



(J.7) 



and, so, (J.5) is compatible with (J.4) if the following identity is verified: 

wi = Wa - yl ^ Wi - Wf. 

This identity is in fact specific of order-four operators conjugated to their adjoint (see 
below (J. 17)). Therefore we prefer to use definition (J.4) for the Yukawa coupling, 
instead of the more restricted definition (J.5). 

Let us assume that the pullback p{x) has a series expansion of the form 

p{x) = \ - x'' - A{x), (J. 8) 

where the exponent r is an integer, where A is a constant, and where A{x) is a 
function analytic at .t = with the series expansion: 

A{x) = 1 + ai - x + a2 - x^ + • • • 

X For an order-four operator the Wronskian is Wi. 
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The determinantal variables W^'s transform very nicely under puUbacks p{x) of the 
form (J. 8): 

{Wi{x), W2{x), W3(x), Wi{x)) (J.9) 
(W^{p{x)), ^ ■ W^ipix)), • Wsipix)), ^ • W,{p{x))), p' = 

One can show that the nome (51) of an order- A'^ operator transforms under a 
pullback p{x): 

q{x) — > Q{x) with: X-QixY = q{p{x)). (J. 10) 

From the covariance property (J.9), and from the previous transformation q — > 
A • q'^ for the nome, one easily gets the transformation of the Yukawa coupling seen as 
a function of the nome K{q) — > K{X ■ q''): 

- . Kiqipix))) ^ KiX.QixY). 

For A = 1 and r = 1 (i.e. when the pullback is a deformation of the 
identity transformation), one recovers the known invariance of the Yukawa coupling 
by puUbacks (see Proposition 3 in [154]). 

One finds another pullback invariant ratio, namely: 

K* = ' I , (J.12) 

which is, in fact, nothing but the Yukawa coupling for the adjoint of the original 
operator. 

Another invariance property is worth noting. Let us consider two linear 
differential operators ili and ^2 of order that are equivalent, in the sense of 
the equivalence of linear differential operators. This means that there exists linear 
differential operators intertwiners Ii, I2, Ji, J2, of order at most — 1 such that 

VLi-h = 12-^2, and: Ji ■ Qi = SI2 ■ ^2- (J. 13) 

Let us assume that one of these intertwiners is a linear differential operator of order 
zero (a function), then the previous homomorphism between operators amounts to 
saying that the two operators are conjugated by a function: 

Q2 = pix) ■ ni ■ p{x)-\ (J. 14) 

which correspond to changing the four solutions as follows: — > p{x) ■ yi. In 
such a case (quite frequent as will be seen in forthcoming publication) the previous 
determinant variables transform, again, very nicely under the "gauge" function p{x): 

[Wi, W2, W3, Wi) ^ (pix) ■ Wi, p{xf ■ W2, p{xf ■ W^, p{xf ■ Wi). (J.15) 

It is straightforward to see that the Yukawa coupling and the "dual Yukawa" , 
K* , are invariant by such a transformational. Two conjugated operators (J. 14) 
automatically have the same Yukawa coupling. 

Do note that the Yukawa couplings for two operators, which arc non trivially 
homomorphic to each other (intertwiners of order one, two, ...), are actually different. 

^ K and K* (and their combinations) arc the only monomials W^^ to be invariant 

by (J.9) and (J.15). 



Diagonals of rational functions 



76 



The (pullback invariant) Yukawa coupling is not preserved by operator equivalence (see 
subsection (8.9.3)). 

Remark: The definition of these determinantal variables WiS heavily relies on 
the MUM structure of the operatorf. The four solutions are not on the same footing: 
the log filtration imposes a natural order between the four solutions the definition 
of WiS relies on. It is worth noting that if one permutes the four solution j/i, one 
would get 24 other sets of (VFi, W2, W3, W4) which are actually also nicely covariant 
by pullbacks, thus yielding a finite set of other "Yukawa couplings" or adjoint Yukawa 
coupling K* also invariant by pullbacks. In fact these "Yukawa couplings" , and other 
K*, can even be defined when the linear differential operator is not MUM. and they 
are still invariant by pullbacks. 

Appendix J.l. Pullback-invariants for higher order ODEs 

These simple calculations can straightforwardly be generalised to higher order linear 
differential equations. We give here the invariants for higher order linear differential 
operators. 

Let us give, for the n-th order linear differential operator the list of the Kn 
invariants by pullback transformations: 

_ ■ W3 _ W! ■ Wj _ Wl'^ ■ W5 _ W?^ ■ We 

A n-th order linear differential operator has Kn as an invariant by pullback 
transformation, as well as all the with m < n. K3 is the Yukawa coupling, 
and one remarks, for the order-four operators, that the other pullback invariant K* 
(see (J. 12)), which is actually also the Yukawa coupling of the adjoint operator, is 
nothing but K^/K^. 

For order-four operators conjugated to their adjoint (see (J. 14)) (i.e. operators 
homomorphic to their adjoint, the intertwincr being an order zero differential operator, 
a function), one has the equality 

A'4 = A'l, i.e. K3 = K\ or W| = ■ W^, (J. 17) 

to be compared with the equality in Almkvist ct al. (see Proposition 2 in [88]) 

2/02/3- 2/3 2/0 = 2/1 2/i- 2/2 2/i, (J-18) 

which is satisfied when the Calabi-Yau condition that the exterior square is of order 
five is satisfied. 

If a linear differential operator ^4 verifies condition (J. 17), its conjugate by a 
function, p{x) ■ ■ p{x)~^, also verifies condition (J. 17) (their Yukawa couplings are 
equal) . 

The condition (J. 17) is not satisfied for linear differential operators homomorphic 
to their adjoint with non-trivial intertwiner (of order greater than zero). For instance 
the order-four operator (171) does not satisfy condition (J. 17). 



f Or recalling the powerful result of Steenbrink [143] that the filtration by the logs is the same as 
the Hodge filtration, the definition relies on this filtration. 
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Remark: Concerning order-three operators and hypergeometric functions. It is 
worth noting the four examples 

3i^2([i i [1, 1], 64x), 3i^2([i i [1, 1], 108x), 

3i^2([^, i ll [1, 1], 256z), 3i^2([^, ^, [1, 1], 1728x), 

which are such that their series expansions, as weU as their associated nome, mirror 
map (compositional inverse of the nome), are series with integer coefficients, the 
previous invariant being the constant 1. 

This however comes as no surprise since the four corresponding operators are 
all symmetric squaresf of order-two operators. Their solutions of the form yoix) = 
u{x)'^, yi{x) = u{x) • v{x), and y2ix) = v{xY /2, give automatically = 1. 

Appendix K. Quasi-Calabi-Yau ODE associated to the Hadamard 
product of two HeunG functions 

The operator having the Hadamard product of the two HeunG functions 
HeunG{a, q, 1, 1, 1, 1; x) and HeunG{A, Q, 1, 1, 1, 1; x) as a solution reads: 

{x -l)ix~ a) {x - A){x - Aa){Aa - x'^f ■ ■ (K.l) 
+ 2 (x^ - Aa) -Us- x^ ■ Dl ~ U2 ■ x ■ -Ui-D^ +U^, 

where: 

[/3 = 5 - 4 (a + 1) + 1) • x^ + {iaA^ + ia^A- Aa + iA + ia) ■ x^ 
+ AAa{a + l){A + l) ■ x^ - Aa (5 a'^A + 5 aA^ + 9Aa + 5a + 5A) ■ x^ 
+ AA^a^ {a+l){A + l)- X -3a^A^ 

U2 = -25a;*^ + (14^a+(14 + Q)a+(14 + g)yl+14 + g + Q)- 

-2 (iAa{a + A) + Aa{Q + q-2Q) + 'iA + qA + ia + Qa) ■ x^ 

-Aa{4:2Aa + A2A + qA + Qa + A2a + Q + q + A2)-x^ 

+ 2 Aa{llAa{a + A) + 2 Aa{q + Q -2) + 11 a +\lA + 2Qa + 2 qA) ■ x'^ 

+ A^a'^{-iQAa + ^Oa + iOA-Qa~qA-Q-q + iO)- x^ 

-2A^a^ [12 Aa{a + A) + Aa{Q + q + 2b) + Qa + qA + l2a+\2 A) ■ x^ 

+ J^a^(\\Aa + 14a +14 A + Qa + + Q + g + 14) • x -7 A^a^, 

Ui = -15x^ +2{2Aa + 2A + 2a + Qa + qA + Q + q + 2)- x'^ 
-2 {Aa{Q + q- 23) + qA + Qa) ■ x^ 
-QAa{2Aa + 2A + 2a + Qa + qA + Q + q + 2)- x^ 
- A Aa {Aa{a + A) - 2 Aa{Q + q - ^) - 2Qa - 2 qA + a + A) ■ x^ 
+ 2A^a^ {ISAa +Qa + qA + ld,a+lSA + Q + q+l%) - x^ 
-QA^a^ {2Aa{a + A) + Aa{Q + q + b) + 2 A + 2a + Qa + qA) ■ x^ 
+ 2 A^a^ {2Aa + 2A + 2a + Qa + qA + Q + q + 2)- X - A'^a'^, 

t Recall that 3-P2{[l/2, a, 1 - a], [1, 1], x) = 2Fi{[a/2, 1/2 - a/2], [1], xf. 
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Uo = a;^ -gQa;6 - (-g2^ +3Aa -aQ2) . a;5 

-aA{2qA +2q +2Qa +2Q - qQ) ■ 
+ Aa(3Aa(2g + 2(7 + l) + 6ga + 6 - 2 aQ^ - 2 A) 
- A^a^ (6 + 6 q + 6 Qa + 6 Q - gQ) • 

+ A^a^ (A a (2 Q + 2 g - 1) + 2 Qa + 2 + aQ^ + A) • x - qQ . 

This order-four operator satisfying the Calabi-Yau condition (its exterior square is of 
order five). The solution of this order-four operator, analytic at x = reads: 

Qq (2 ga + 2 g + ^2 _ a) (2 + 2 g + Q2 _ A) , 

= ^ + AT'"- + 16^51^ + 

the nome reads 

and the first terms of its Yukawa coupling K{x) read: 

Even inside this restricted set of HcunG functions solutions of the form 
HeunG{a, g, 1, 1, 1, 1; a;) it is hard to find exhaustively the values of the 
two parameters a, and of the accessory parameter g, such that the se- 
ries HeunG{a, g, 1, 1, 1, 1; x) is globally bounded. The HeunG function 
HeunG{a, q, 1, 1, 1, 1; 6 a;) becomes a series 1 + A^i a; + N-zx'^ + ■ ■ ■ with integer 
coefficients A^i, N2, for 

^ 7V2.(2g-l) ' ' q ^ ' 

These are necessary conditions. One can find some other necessary conditions on 
the integer coefficients of the series. Besides g = 1/2, a = 2q, a = l/3g(13g — 
6)/(2g — 1) (i.e. N2/NI = 5/Aoi N2/N^ = 4/3) that require some specific analysis, 
one finds that one has necessarily the following relation among the first four coefficients 
Ni, N2, N3, N4, namely: 

27 • + 10 • {N^ -9N2)- Ni- N3 + 8 • {AN^ -3N2)- 

-9- Ni - {Nl ~6N2) = 0. (K.5) 

One verifies easily that (K.5) is actually verified for (52) and (100). 

Note that the ratio Np/Nf are rational expressions of a and g, invariant by 
(a, g) — >■ (1/a, g/a), of the form P{a, q)/q^, where P(a, g) is a polynomial. 

However, finding the values of the accessory parameter g, such that 
HeunG{a, g, 1, 1, 1, 1; x) is globally bounded remains difficult. Is it possible to 
finds such values of the accessory parameter g for any given integers Ni and 
N2 ? For instance for Ni = 2, N2 = 10 the series with integer coefficients 
HeunG{—8, q, 1, 1, 1, 1; 8 a;) is actually a series with integer coefficients for g = —2. 

Furthermore it is difficult to find the values of the two parameters a and g 
such that the order-two operator having HeunG{a, g, 1, 1, 1, 1; a;) as a solution is 
globally nilpotent. For instance, if one restricts to a = 9/8, it is hard to show that 
the only rational number value of the accessory parameter g yielding global nilpotcnce 
is g = 3/4. For g a rational number one gets an infinite set of divisibility conditions. 
The prime 5 must divide the numerator of (g + 3) (g + 1) (^ + 2) (g^ + 4g + l), the prime 
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7 must divide the numerator of (q^ + 6q + 6) {q + Q) {q + 1) (q + 4) (q^ + 4 g + 1), etc. 
With q ~ N/D where N and D integers < 1000 and with the conditions emerging 
from the p-curvature calculations for the first primes < 42. One finds that q ~ 3/4 
is the only value of the accessory parameter corresponding to global nilpotencc. 

Restricting to HeunG functions solutions of the form HeunG{a, q, 1, 1, 1, 1; 60 ax) 
for integer values of a and q it is hard to find the integer values of a and q such that 
the corresponding series are globally bounded. 

Appendix L. Yukawa coupling of Calabi-Yau ODEs 

Let us give the expansion of the Yukawa coupling for a set of other Hm.n that are 
Calabi-Yau: in particular they are MUM and their exterior squares are of order five 
(the "Calabi-Yau condition" ) . 

For H4,s the Yukawa coupling reads: 

K{q) ^ K*{q) = 1 -f 16 • g -f 352 • + 33280 • + 2058528 • g'' 

4- 123766016- q^ + 7347718144- q^ + 439489011712- q^ 

26579639900960- q^ -f 1616513123552128- q^ +•••, (L.l) 

which is Number 36 in Almkvist et al. large tables of Calabi-Yau ODEs [91]. 
For i?4,9 the Yukawa coupling reads: 

A'(q) = K*{q) = 1 + 12 - q - 324 - q^ - 29544- q^ - 1314756 - q" 

- 12971988 - q^ + 2033927928 - q^ + 146587697352 - q^ 

+ 4172739566652- q^ - 77469253445544- q^ + -•- , (L.2) 

which is Number 133 in tables [91]. 
For i?6,6 it reads: 

A(q) = A'*(q) = 1 + 37 - q -4523- q^ +412327- q^ - 33924139 - q^ 
2662557912 - q^ - 203154013049 - q^ + 15217617773948- q^ 

- 1125153432893483 - q^ +82390368380951296- q^ + -- - , (L.3) 

which is Number 144 in tables [91]. 
For i/g.s it reads: 

A(q) A*(q) = 1 + 24 - q -2012- q^ + 139056- q^ - 8227932 - q^ 
+ 468328024 - q^ - 25856580632 - q^ + 1402012096656 - q^ 

- 74994891745116 - q® + 3972880128014736 - q^ +---, (L.4) 

which is Number 176 in tables [91]. 
For iJg.g it reads: 

A(q) = K*{q) = 1 + 18 - q - 1674- q^ +88209- q^ -4801770 - q" 
+ 239279643- q^ - 11680323039 - q^ + 558685593414- q^ 

- 26379917556714 - q^ + 1233104626297710 - q'' +•••, (L.5) 

which is Number 178 in tables [91]. 
For Hs,8- 

K{q) = K*{q) = 1 + 16 - q - 864 - q^ +47104- q^ - 1890528 - q"" 
+ 80502016- q^ - 3118639104 - q^ + 123287486464- q^ 

- 4691784791264 - q^ + 179585946086272 - q^ + - - • , (L.6) 
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which is Number 107 in tables [91]. 

For i?8,9 the Yukawa coupling reads: 

K{q) = K*{q) = 1 + 12 • g - 756 • + 27192 • - 1144644 • q^ 

+ 39948012 • q^ - 1377082728 • q^ + 47882164776 • q'^ 

- 1608623259588 • +53732432848152-99 +■■■, (L.7) 
which is Number 163 in tables [91], and for iJg.g it reads: 

K{q) = K*{q) = 1 + 9 • g - 567 • + 20205 • q^ - 615735 • 
+ 19431009- q^ - 608213043 • q^ + 18406651167- 

- 542566460727 • g® + 15865350996861 • g^ + • • • , (L.8) 
which is Number 165 in tables [91]. 

Appendix M. Modular form character of 2^"i([l/6, 1/6], [1], x) 

The modular form character of (152) is clear on the remarkable (and intriguing) 
identity (58) (or (59)) in Maier [123]: 

,F,{[l,l],[l],V,) (M.l) 

/ (a; + 60) (a; + 80) (a: + 96) \-i/6 / I 1 
■ I (x- + 48) (x + 120)2 J ■ ' H^6' 6^ 

where the two puUbacks read respectively (o denotes the composition of functions) 

X (x + 60)2 (x + 72)2 (x + 96) 
' 432 (x + 48)(x + 80)(a;+ 120)2 ^ 
X 1728 a; ^ x 1728 a;2 642 



a; - 1 (a; + 16)3 a; - 1 (a; + 256)3 G ' 

^ a:2(x + 48)2(x + 72)(a; + 120) 

^ 432 (.T + 60) (a; + 80)2 (x + 96)2 ■ > 

X 1728 a; 642 ^ 1728^2 

o — ■ o = o — ■ o H, 



where 



a;-l (a; +16)3 H x-l (a; + 256)3 

^ 8x- (a; + 96) ^ ^ -64a;. (a; + 48) 

(.T + 72)(a; + 60)' (a; + 72) (a; + 120) ' ^ ' ' 

The algebraic relation between these two pullbacks {u, v) = {Vi, V2), corresponds 
to the (genus zero) modular curve: 

6912 uu- {u'^v* + l) ~3uv- (80{u^+ v^) + 529999 uu) • (v^u^ + l) 

+ 16 (u + u) • (uv + l) ■ {u^ + v'^ - 241 u v) ■ (v'^ -uv +1) 

~6uv- (u + v)- {uv + l)- (6A0{u'^ + v^) - 652959 uw^ 

+ UV- (^6912 {u^ + v^) - 1589997 (u^ v + v^ u) + 21958300 v^ u^^ = 0. 
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Other identities are worth noting on these pullbacks, for instance: 

-Ax ^ 1728a; ^ 1728a;2 4096 

o Vi = 7 ——T o t = o — — , and: 



(1 - xf [x + 16)3 (a; + 256)3 F 

-Ax 1728a;2 1728a' 4096 

o Vf = -, TT ° F = -rr o , where: 



(l-x)2 - (x + 256)3 (a; +16)3 F 

^ _ a; • (a; + 80) (x + 48) {x + 96) 
(a; + 120) (x + 72) (x + 60) 
It can also be illustrated through the identity: 

p (,l il Ml ^ 



(1 - x)i/6 (1 + 8x)i/3 /I 5^ 27x2(1 -x)2(l + 8x) 



(l + 7x+x2)i/4 ^ H42' 12^' ^ J' 4(1 + 7x 
(1 - x)i/6 (1 + 8x)i/3 /I 5 -108x(l -x)(l + 8x)2 



(1 - 4x)i/2 'Vl2' 12" ' " (1 - 4x) 

where the last two pullbacks are related by the fundamental modular curve [17, 32, 123] 
(corresponding to r — > 2r), as can be deduced from the identities on these two 
pullbacks: 

27 x2 • (1 - x)2 (1 + 8 x) 1728x2 16 (1 + 8 x) 



o 



4(l + 7x+x2)3 (x + 256)3 (x-l)x ' 

108x- (1 -x)(l + 8x)2 _ 1728x 16(l + 8x) 



o 



4x)6 (x + 16)3 (x-l)x 



(M.6) 



(M.7) 



The relation between the hypergeometric function 2^'i([l/6, 1/6], [1], x) and 
2Fi([l/12, 5/12], [1], x), can also be understood from the Kummers's quadratic relation 
on 2Fi([l/6,5/6],[l],x): 



.FAiiliiH^) - .F,([l,A],[l],4x.(l 



Appendix N. Calabi-Yau condition versus integrality: a Saalschutzian 
hypergeometric family of operators 

The Calabi-Yau condition that the exterior square of an order-four linear differential 
operator is of order five is a fundamental condition defining Calabi-Yau ODEs [88, 90, 
91]. Let us consider the following Saalschutzian hypergeometric function 

iFsdb - c + d, b, c, d], [e, b + d, 1 + b + d - e], x), (N.l) 

and, let us introduce the order- four linear differential operator M4(6, c, d, e) which 
annihilates this hypergeometric function. It is a straightforward exercise to verify that 
this linear differential operator satisfies the Calabi-Yau condition: its exterior square 
of is actually of order five for any values of the parameters b, c, d, e. This order- four 
operator is almost self-adjointf : 

^b+d^i . j^j^^f^^ ^ adjoint{AU{b, c, d, e)) • x^+'*"^ (N.2) 

t If one denotes by N the denominator of the rational number b + d—l , the symmetric A^-th power 
of M4{b, c, d, e) is up to a normalisation, self-adjoint. 
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For generic rational parameters, the series expansion of the Saalschutzian 
hypergeometric function (N.l) is not globally bounded. 
Let us restrict to the condition d = 1—6: 

4F3([1 - c, 6, c, 1 - 6], [e, 1, 2 - e], x). (N.3) 

The order-four operator is still such that its exterior square is of order five, but it 
is now a self-adjoint operator (see (N.2)). Again, one sees that, for generic rational 
parameters, the series expansion of the Saalschutzian hypergeometric function (N.3) 
is not globally bounded. Actually the four solutions are not of a MUM form. Besides 
the analytic a.t x — 0, non globally bounded, solution (N.3) 

bed - (b — c + d) .\ 
' + e.(5+rf).(l+5+rf-e) - + (^-^^ 

the three other solutions have no logarithm and are, for generic rational parameters, 
Puiseux series: 

i_e (. (1 + - e) (1 + c - e) (1 + b - e) (1 + b - c + - e) 
^ 'V {e-2){l +b + d-e){b + d + 2- 2e) 

x^-''-''-{\ +•••), x'-''-''-{l + •••)• 

When the exponents 1 — e, 1 — 6 — d, e ~ b — d are integer values one can recover 
solutions with logarithms for the operators. 

Appendix O. Integrality of a one-parameter Saalschutzian 
hypergeometric family of operators 

Let us consider the order- four linear differential operator {0 = x ■ Dx): 

Af4(/i) = 16 • 6*2 . (6» - 1)2 (0.1) 

- x • {2 6 +l f ■ (2 6* - 1 + ^) • {26 -1 - fl). 

We have a one-parameter family of operators depending on fi^. 

The Wronskian of M4{ii) is independent of /i, and reads l/{x — l)'^/x^. This 
order-four operator is, for any rational value of /i, non-trivially homomorphic to its 
adjoint (with order- two intertwiners) : 

adjoint{Mi{pL)) ■ {x - 1) • L2{fi) = L2{fi) ■ {x - I) ■ Ahifi), (0.2) 

A'Uifi) • (a; - 1) • M2{n) = AMfi) ■ {x - 1) ■ adjoint{Mi{fi)), (0.3) 
where L2{fJ-) and M2{p-) are two self-adjoint order-two linear differential operators: 

L2it^) = Dl /'f^ , (0.4) 

4 /i^ {x — 1) • a;^ 

M2{p) = IQ-x'^-Dl +64-a;3.D^ +36x2+- 

{x — 1) ■ /i^ 

The order- two operator M2{p-) has simple hypergeometric solutions: 



2 /i 2 /i fl 

(1 • . ,fM^±1, ^±1], [^], x). 

\ 2 fl 2 fl fl J 



Diagonals of rational functions 



83 



M4(/i) = Ahiti) ■ {1 - x) ■ L2{fi) + , T n \.^ ■ (0-5) 



note that, gcncrically, the hypergeometric function 2-^1 ([ ^^^ , ^-f^], ^) does 

not corresponds to a globally hounded series, as can be seen, for instance, for ~ 
-5/11. 

The two intertwining relations (0.2) and (0.3) can, in fact, be seen as a straight 
consequence of the fact that the order-four operator M4 [fx] can remarkably be written 
in terms of these two self-dual operators: 

This order-four operator (0.1), or (0.5), does not satisfy the Calahi-Yau 
condition. Its exterior square AAq[^) is of order six, with a, not only rational function 
solution, but a constant solution. It is actually the direct sum of Dx and of an order- 
five hnear differential operator M.^{pL): 

Mein) = Dx ® M^iii), where: (0.6) 

Mz{^Ji) = 7^5(0) + • Mf + / • M'-i^ + / • M'i^ • i?, + / • • Dx. 

where the A^r™'''s are linear differential operators of order n. For /i^ = 1 the order- 
five operator A^5(/x) becomes the product:]: of two order- two operators and of Dx'. 

M^[±l) = N2-M2-Dx. (0.7) 

For the other odd values of {fi = ±3,±5, •••) the order- five operator M^ifi) 
factorizes into the product of an order-two, order-one and an order-two operator. 

As will be seen in a forthcoming publication, the fact that the exterior square 
has a rational solution is also a consequence of the decomposition (0.5). We have the 
following general result. Any order-four linear differential operator of the form 

M4 - M2 • co(x) • L2 + (0.8) 

Co (a;) 

where L2 and A/2 are two (general) self-adjoint operators 

L2 = a2{x) ■ Dl + • Dx + a^{x), (0.9) 

ax 

M2 = b2{x) ■ Dl + • Dx + bo{x), (O.IO) 

da; 

is such that its exterior square has A/a2{x) (A is any constant) as a solution. In the 
case (0.5), the solution is the constant solution: 1/02(2;) = 1. Instead of A/4, we 
can introduce 

M4 = co{x) ■ M4 = co{x) ■ M2 ■ cq{x) -12 + A, (0.11) 

which annihilates the same solutions as M4. It is worth noting that a decomposition 
like (0.8), provides, in fact, interesting results on the spectrum of A/4. For simplicity 
let us consider the example (0.1) with its decomposition (0.5), or on Al4{fi), the 
decomposition 

Ahil^) = {l~x)-M4 = Miifi) + (0.12) 

where the order-four operator A/4(/i) factors: 

mifi) = (1 - x) ■ Ahifi) • (1 - x) • L2ifi). (0.13) 

X It is also a direct sum M2 ffi {N2 ■ D^) where M2 and N2 are two order- two operators. 
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This order-four operator Mi{^) is, also, (non-trivially) homomorphic to its adjoint: 

adjoint{Mi{n)) ■ L2{ti) = i2(/^) • Mi{^i). (0.14) 

It is clear that /i = ± 1 needs to be analysed separately. For instance, for /i = 1, the 
order- four operator M/i{^) factors in a direct sum © M3 where the order-three 
operator has the three solutions: 
5 5, 



9x2. 4F3([1, 1, -,-], [2,3,3], x) + lex- ln(x) _ (q^^^^ 

/TI" „ /rl 1, r 1 \ dX 



Let us consider the solutions of L2{fi) which can actually be expressed in terms 
of hypergeometric functions: 



= (1 -.) ■ .(2^-^^)/^/^ ■ 2fJ[ '^' +^ -'^'^ -^ ], [^], 

V 4/i 4/i 2 ^ / 

where p reads p = 2 • (2 /i^ — 1)^/^. One immediately deduces from the 
decomposition (0.12), that ^'i and ^1 are two eigenf unctions] of the order-four 
operator M4(/^), with the same eigenvalue (/i^ — l)^//!*: 

Af4(Ai)-*. = 1 = 1,2. (0.16) 

Introducing a rational parametrisation for /i and p, namely 

u2 -2u -h2 ^ u2 _4,^^2 ,^ , 

' ^= .2-2 ' (^-^^^ 

one finds that, even for rational values of the parameter u, that forces /i and p 
to he rational numbers as well, the hypergeometric functions in the two previous 
eigenfunctions ^'^ do not correspond (generically) to globally bounded series. More 
generally, for rational values of the parameter u, the four solutions of Mji{fi) are 
Puiseux series of the form ■ A{x), where r is a rational number and A{x) are series 
analytic at x = 0. None of the four A{x) corresponds to a globally bounded series. 

Appendix 0.1. Saalschutzian hypergeometric solution 

In fact the order- four linear differential operator Mi{p), or M4{p), has the 4F3 
(Saalschutzian) hypergeometric solution 

X- ,Fs{[^^,^-±^,l,ll[Wlx), (0.18) 

which expands as: 

9 25 

^ "64^^"^^^^ + ^^'^' +^ (m-1)(m+1)(m-3)(m + 3)- 

1 995 

{p-l){^i + ^){^l^3){p + 3){p-5){fi + 5) ■ x^ (o.i9) 



9437184 



441 

(/i-l)(Ai + l)(M-3)(M + 3)(Ai-5)(/. + 5)(/i-7)(M + 7)- + 



268435456 

f See the spectral theory of ordinary linear differential equations, in particular for self-adjoint 
operators. 
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On this expansion it is clear that the hypergeometric function truncates into a 
polynomial for fi any odd integer (positive or negative). Furthermore, this series 
(0.19) is globally bounded for any rational number ^l. 

Let us take a simple rational value for /Lt, namely /i = 1/3. The series expansion 
(0.19) reads 

In 125 o 42875 - 94325 . 41544503 . 

X + -x^ H x-^ H x^ H x^ ^ x^ + • • • (O.20) 

8 2592 1679616 5971968 3869835264 ^ ' 

The series expansion of the nome (as well as the mirror map) is not globally bounded: 

17 . 22195 . 1913687 , 2195016283 , 

= a; H X ^ x* H x H x^ + • • • 

^ 48 124416 17915904 30958682112 

In contrast, the series expansion of the Yukawa coupling K{x) is a globally bounded 

series 

1 43 11 31517 2 522821 3 

^^'^ X " 144 " 288 ^ 3359232 ^ ~ 3869835264 ^ ^ "' 
Actually the rescaling a; — )• 2^ • 3'^ • a; changes 2** • 3^^ • K{x) into a series with integer 
coefficients: 

2^ • 3^ • X(2^ • 3^ • a;) = - - 129 - 7128 a; - 756408 - 4705389 a;^ 

X 

+ 58331013489 x^ + 38259799407522 a;'^ + 19576957591348938 x^ (0.21) 
+ 9193736880930978297 a;^ +4149261387452007788523 a;* + ••• 

Appendix 0.2. A non-trivially equivalent operator for /i = ± 1 
Let us introduce the order-four operator 

Ni = 16-6I- (6* -2) • (6* - 1)2 (0.22) 

- A - X - e ■ {e - I) - {2e - if. 

which is nothing but M2 ■ D'^ (see (0.7)). 

This operator is non-trivially homomorphie to Ma{± 1), that is (0.1) for = ±1: 

Ni-A2 = (^) • A2 ■ (i^) • M4(± 1), (0.23) 
where: A2 = {1 - x) ■ 9 ■ {d - 1). 

Operator A'4 does verify the Calabi- Yau condition: its exterior square is of order five. 
The operator iV4 can also be written as a direct sum: 

1 



or, in terms of = x ■ D^, the direct sum: 

e ® {9 - I) ® {iQ - 9- [e - I) ~ Ax - {26 ~ 1)2). (0.25) 

Thus, besides the constant solution and y{x) = x, its solutions can simply be written 
in terms of hypergeometric functions, for instance 

x-{l-x)- 2F^([%ll [2];x) (0.26) 



2' 2' 

4.-(i-.)-|:(.^r([iiL[i];. 
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This series is globally bounded. Changing a: — > 16 x turns this series into a series 
with integer coefficients: 

16 X +32x^ + 192. + 1600. t"* + 15680 + 169344 + 1951488 

+ 23557248 X* + 294465600 x^ + ••• (0.27) 

Appendix 0.3. A family of non-trivially equivalent operators for odd integer fj. 
It is tempting to try to generaUse (0.24), or (0.25), with a simple /i^-ansatz: 

Dl e M2(p) where: 7^2 (a*) = DI - , _ (0.28) 
or, in terms ol 9 ~ x ■ Dx, the direct sum: 

e ® [6 - I) ® (l6- 61 • (6* - 1) - 4x- (46*2 - 46* + ^?)^. (0.29) 

In order to compare this ansatz operator with the initial operator (0.1), we shghtly 
rewrite it as 

16 • 6*2 . (6» - 1)2 _ a; • (2 6* + 1)^ • {AO"^ - 40 + 1 - ^?), 
and we, also, rewrite this ansatz operator as follows: 

Codd{fi) = l&-e- {9 - If ■ {0 - 2) (O.30) 
- 4 • X • 6* • (6* - 1) • (40^ - 46* + n^). 
One finds that such an ansatz is actually non-triviallyf homomorphic to operator 
(0.1) for any odd integer values (positive or negative) of fi, and that its exterior square 
is actually of order-five for any value of ^ (Calabi-Yau condition). 

For fJ, — 2i/(l + t^), with t a rational number (hence < 1, the solutions 
of (O.30) read (besides the constant solution and y{x) = x) simple hypergeometric 
functions, for instance: 

Soix) = X ■ 2Fi{[j^, [2], x), (0.31) 

together with the solution§ Si{x) ~ So{x) •ln(x) + Si{x) where iSi(x) is analytic at 
X ~ 0, and solution of an order-four operator Af^, product of two order-two operators, 
Af4 = A/'2 (/.*)• 2(^)5 where A/'2(/i) is an order- two operator homomorphic to A^2(m) 
(see (0.28)): 

AA2(a.) • 4- (2^^-l) = -■(Dx-^^^^^)-M2{t^), (0.32) 
x^ X V ax y 

, , , 1 — X -I- ^2 X 

where: p{x) = ^ jj^. 

(1 — x) • X'^''' 

The series expansion of (0.31) is globally bounded for any rational valued of the 
parameter t: 

+2/0+2 I 1 W+2 



1 t^ 2 1 r (2t2 + 1) {t^ + 2) 3 

^2 (l+i2)2 • ^ +U {l+W 

1 t2(3f2 + 2)(2f2 + 3)(2f2 + l)(f2 + 2) 
"^144 (l + i2)6 



„4 



X 



(0.33) 



f The intertwiners are linear differential operators of order three. 

§ This solution can also be written MeijerG{[[], [{2 + t'^)/{l + t'^), {2t^ + 1)/{1 + t^)]], [[0,1], []],x), 
i.e. as a MeijerG function. 

^ In contrast the series S\ {x) is not globally bounded for the rational values of t. 
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For instance for fi ~ 4/5 the hypergeometric solution (0.31) reads So{x) — 
X ■ 2^i([l/5, 4/5], [2], x), corresponding to the globally bounded solution series: 

2 , 18 . 231 4 17556 , 1474704 . 

X H x^ H x-^ H x^ H x^ H x^ + • • • 

25 625 15625 1953125 244140625 

The rescaling x — >■ 5'^ x turns this series into a series with integer coefficients. 

Remark that the two solutions Sq{x) and Si{x), for ^ = 4/5, can be replaced 
by the two solutions well-suited for x large (z = ^/x): 

The two hypergeometric functions 2i^i([l/5, 6/5], [2/5], z) and 2-Fi([4/5, /5], [8/5], z) 
do not correspond to globally bounded scries. We have a similar result for the other 
rational values of t. 

Do note, however, that for the other rational values of ^ we have a drastically 
different situation: the operator (O.30) is no longer globally nilpotent^, as can be seen 
on the solution of D'^ - fi'^/4:/x/{l - x) in (0.28) 

(i-.).,^.( r i^ui (0.34) 



which has the series expansion 
^ +1^^'- ('^'+8)--' +^A^^.(/.2 + 8)(m2 + 24).x4 

+ • (m' + 8) (m' + 24) (/? + 48) • + • • • (0.35) 

For rational values of fi that are not of the form /i = 2t/(l + t^), the solution-series 
(0.35) of the operator (O.30) is not globally bounded as can be checked with the two 
rational values of ^ such that \fj.\ < 1 and > 1 respectively. One can verify that 
the solution-series (0.35) of (O.30) for ^ = 1/3 

1 2 73 3 15841 4 6859153 5 4945449313 g 
72 15552 ^ ^ 6718464 ^ ~^ 4837294080 ^ ^ 5224277606400 ^ " ' 
and for /i = 3: 

9 o 51 . 561 4 31977 . 948651 . 40791993 7 

X + -x^ -\ x"^ ^ X* H x^ H x^ H x^ + ■ ■■ 

8 64 1024 81920 3276800 183500800 

are not globally bounded. 



Appendix 0.4- Seeking for equivalent operators for other values of fi 

The operator (O.30), which is non-trivially homomorphic to (0.1), (or (O.30)) for 
odd integer values of the parameter ^, is not valid for even integer values of /j,. For 
example, for /i = 0, operator (0.1) actually verifies the Calabi-Yau condition, its 
exterior square being nothing but M^in) for = (see (0.6)). Operator A^5(0) 
reads 

16 • • {9 - 1)2 - X ■ {20 + 1)2 . {29 - 1)2. (0.36) 

It is not of the form (O.30), which cannot encapsulate all the operators homomorphic 
to (O.30) satisfying the Calabi-Yau condition. As previously remarked, the series 

I For instance one sees explicitly on (0.34) that its exponents are not rational numbers, therefore 
operator (O.30) cannot be a globally nilpotent operator. 
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expansion of the solution x ■ 4F3^[i, i, |, |], [1,2,2], x) is globally bounded, and 
can be turned into a series with integer coefRcients with the rescaling x — )■ 256 a;. 

Note that (0.36) can be used as a "seed" to get a family of /z^-dependent operator 
verifying the Calabi-Yau condition. The linear differential operator 

C{fi) ^ 16 • 6*2 . (61 - 1)2 (0.37) 
- X - {29 + 1 - /i) • (20 + 1 + /i) • (261 -1 -fi)- {20 - 1 + ^i), 
is such that its exterior square is actually of order five. Very fortunately, operator 
(0.37) is non-trivially homomorphic to (0.1) for any integer value (positive or 
negative, even or odd) of parameter fi. Again the intertwiners are of order-three. 

For instance, for /i = 2 one has the intertwining relation between operator (0.1) 
and (0.37) 

C{2) ■ ((86i2 . (261 - 3) + 461 + 1) - 2x ■ {20 - 3) (261 + 1)^^ (0.38) 

= (^{80^ ■ {20 - 3) + 461 + 1) - 2x ■ {20 - 3) {20 +1){20 +3)) • ]Vh{2). 

For larger values of the integer fi the intertwiners become more and more involved. 
They are still of degree three in 9 but of higher degree in x. 

Operator (0.37) has simple hypergeometric solutions for any value of /i: 

For any rational value of /x, the series expansion of the solution of (0.37) analytic 
at a; = 0, namely (0.39) 

X +^ (A^-l)(M + l)(/i-3)- (/i + 3)a;2 (O.40) 

+ (m - 1) (m + 1) (m - 3)' {^l + 3f{^l-5){^l + 5) ■ x' +■■■ 

is globally bounded^. 

For instance, for ^ = 1/3, the series (O.40) 

10 2 2800 3 1078000 4 53953900 5 26709338656 g 
^"^81^ 59049 ^ 43046721 ^ 3486784401 ^ 2541865828329 
can be turned into a series with integer coefficients after the rescaling a; — > 3^ a;. 
However, the series expansion of the nome for (0.37) for fi = 1/3 is not globally 
bounded: 

19 2 250403 3 218211473 4 281241377443 5 

^^M^ ^ 1417176 ^ ^ 2066242608 ^ ^ 4016775629952 ^ 

1456188325082179 g 3167628271177596809 7 

29282294342350080 ^ ^ 85387170302292833280^ ^ 
In contrast, the Yukawa coupling K{x) for fi = 1/3 is actually globally bounded 
^ 1 49 2179 ^ 1508129 2 47590097 

~ X ^ 162 ^ 59049 ^ ~ 172186884 ^ ^ 167365651248 ^ ^ " ' 
Actually the rescaling a; —>■ 2 -3^ a; on 2- 3^ A' (a;) turns the previous series expansion 
into a series with integer coefficients: 

2- 3^ • K{2- 3^x) = - -441 - 78444 a; - 27146322 a;^ + 1284932619 a;^ 

X 

+ 27674475754905x4 + 59119113109746798 a;^ + 100896041483693939736 a;^ 
+ 158984355721045048019613 a;^ + 241323001023828827752150059 a;^ + ••• 

t Operator (0.37) is thus globally nilpotent for any rational value of fi. 
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to be compared with (0.21) for operator (0.1). We have similar resuh for /i = 4/5 
(the series (O.40) can be turned into a series with integer coefficients after the rescahng 
X ^ 2* • S^x, and the series expansion of the nome is not globahy bounded). 
The Yukawa coupling K(x) for n = 4/5 can be changed into a series with integer 
coefficients: 

2^ ■ 3- 5^ K(2^ ■ 3- 5^ x) ^ - - 1057980 - 92574954000 x (0.41) 

X 

- 51733629839745000x2 + 74509092036686778685920 + ••• 

Note that, for any odd integer, (positive or negative), (O.40) is a terminating 
series, reducing to a polynomial. 

As a byproduct, we see, for any odd integer value of ^, that operator (0.1) is 
(non-trivially) homomorphic to two different operators (O.30) and (0.37), such that 
their exterior square is of order five. As it should these two operators (O.30) and 
(0.37) are, for any odd integer value of /i, homomorphicf (with intcrtwiners of order 
three). For fi ~ 3 one gets the homomorphism: 

CoddiS) ■ (40- {9- If ~ Ax - {9-1)- {9-2)- {0 + 2)^ (0.42) 

= (4{9~2) - {9-lf - X - {0-1)- {40^ -40 + 9)) • C(3). 

For other values of /i one gets slightly more involved intcrtwiners that are no longer 
of degree one in x. 

The homomorphisms, for odd integer values of /i, between a non globally nilpotent 
operator (O.30) and the globally nilpotent operator (0.37), or between operators with 
non globally bounded and globally bounded solutions, may seem misleading. It is 
important to note that, for odd integer values of /i, the operator (0.37) is non longer 
irreducible^ , that the intcrtwiners between (O.30) and (0.37) (sec (0.42)) arc not 
globally nilpotent, and, furthermore, that the globally bounded infinite scries (O.40), 
reduces to a polynomial. The intertwining relation between (O.30) and (0.37) then 
matches this terminating hypcrgcomctric scries (O.40) with the y{x) ~ x solution of 
(O.30). 

Do note that, given an order-four operator such that its exterior sguare has a 
rational solution ("extended Calabi-Yau condition"), finding an order- four operator 
(non-trivially) homomorphic to the first operator such that its exterior sguare is of 
order five (Calabi-Yau condition) is an extremely difficult task^, even if one assumes 
decompositions like (0.8), (0.11). We will address this difficult question of the 
reduction of the "extended Calabi-Yau condition" to the "Calabi-Yau condition" in a 
forthcoming publication. 

f Note that for even integer values of fi, operator (0.37) is irreducible when (O.30) is reducible. 

Therefore operators (O.30) and (0.37) cannot be homomorphic for even integer values of fi. 

t For any odd integer value of fi the operator (0.37) is the product of four order-one operators such 

that their wronskian is a rational function. Operator (0.37) is thus globally nilpotent. 

^ In particular because, as we have seen, this reduction to an operator satisfying the Calabi-Yau 

condition is not unique. In contrast, starting from an operator satisfying the Calabi-Yau condition, 

like (0.37), it is straightforward to get operators that are non-trivially homomorphic to this operator, 

and are such that their exterior square has a rational solution (just perform the LCLM of this operator 

with any order-three linear differential operator). This will be explained in a forthcoming publication. 
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Appendix P. Modular forms and selected 2F1 hypergeometric functions 
with two pullbacks 

We display here a (non exhaustive) hst of miscellaneous identities (between modular 
forms and their representations as 2^1 hypergeometric functions with two pullbacks) 
that we often encountered in our studies of the Ising model, lattice Green functions, 
or Calabi-Yau ODEs: 

2F,([^, 108-x2.(l + 4x)) (P.l) 

^ ' ^ ^V6' 3^' ^ ^' (1 - 12a;)3 

= 1 +6.x2 +24. T''' + 252 +2016 + 19320 a;*^ + 183456 
+ 1823094 .T^ + 18406752 a;^ + 189532980x1° + ••• 



(l + --+a:Y/^-2F,([i,i],[l]; ^ 



1 1^ 2 7a; (1 + x) 
(1 +2a-)6 



The series expansion of (P.2) are globally bounded and can be changed into a series 
with integer coefficients after the rescaling x ^ Ax. This identity is nothing but 
identity (P.l) after a change of variable. Another example corresponds to HeunG 
functions of the form HeunG{a, q, 1, 1, 1, 1; x), such that the two parameters a and 
q are associated with fixed points of the symmetries (139) of these HeunG functions : 

i7ewriG(i, 1,1, 1,1,4 X-) = iJewnG(2, 1, 1, 1, 1, 1, 8 x) 



^- HeunCi-l, 0,1,1,1,1, 

1 -4x \^ 



4x 



4x 



1 -4x " '^^4' 4^' ' " (1 -4x)4 

rl 1 



2F,{[-, -], [l],64x.(l -4x).(l -8x)2) (P.3) 
(1 - 256x+ 5120x^ - 32768 x^ + 65536 x'')"^/* 



1 5, ^^^^ x.(l -4x)-(l -8x) 



2 



^ ^^^1^12' 12^' f^^' "^''^^ (1 - 256x + 5120x2 - 32768x3 +65536x^)3 
(1 - 16 X + 80 x^ - 128 x^ + 256 x'^)-'^/'^ 

X F Al fll 1728 x^.(l-4x)^.(l-8x)2 x 

' 'V42' 12^' ^ J' (1 -16x +80x2 _i28x3 + 256x4)3^* 

= 1 +4x +20. x^ + 112 x^ +676 x^ +4304 x^ + 28496 x^ + 194240 x^ 
+ 1353508 X® + 9593104 x^ + ••• 
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The previous HcunG function can also be written 

2^^if[i i], [l],16x-(l-4x)) = -±-.,FJ[l-,h[ll ' 



^2 



2' 2" ' " ' V 1 -4x ' 'V^2' 2" ' " (1 -4x) 

(1 - Ax) - {1 -8x)- 2F^{[^, [2], 16x.(l -4x)) (P.4) 
2x ^ /,3 3, _ 16x2 



f.(l|.il.|21, 



(l-4x)3 V^2' 2^' ^ " (1 -4a;)2 
Using the identity 

9x2 a;' - (1 +2x)2 

* — 144- 



fl6x • fl - 4x)^ o f — ^ 

V ^ V Vl + 4x + 40x2; (1 + 4x + 40x2)2 

/ 163^2 ^ ^ -3x- (1 +2x) ^ 
= ((T-I^M (l-4x)2 ) (P-^) 



one easily deduces from (P.4) a linear functional identity between the same 
hypergeometric function with three different rational pullbacks: 

18x2. (1 + 4X + 40x2)3. (1 -4x)6. [2], pi(x)) 

- (1 +2xf . (1 -4x)^. [2], p,ix)) (P.6) 

+ (1 +8x)2 . (l+4x + 40x2)3. (1 +2x)*'. 2i^i([| [2], P3(x)) = 0, 
where the three pullbacks read respectively: 

16x2 N / 9x2 
(1 +2x)4 V(l -4x)2/ ° VI + 4x + 40x2- 

x2 . (1 +2x)2 



. N „„„ ( 16x2 N 

Pi X = 1296.- — = -TT 

^ 1 +2x 4 V 1 -4x 2y 



'P2W1 = 144 



(1 +4x +40x2)2 

'l6x . (1 - 4x 

X. (1 +2x) . (l+4x + 40x2 
(1 -4x) 



9x2 



1 + 4 X + 40 x2 



P3(x) = -48 . ^ Z ' 



( , A / -3x. (1 +2x)\ 

(l-4x)2 )■ 



Note that each of the three terms in (P.6) corresponds to series with integer coefficients, 
their series reading respectively: 

18 x2 - 216 x^ + 2160 x^ - 25344 x^ + 223236 x^ - 1810512 x^ + . . . , 

- 1 + 14x - 190x2 + 2524x3 - 26732 x" + 270184 x'^ - 2650164 x^ + . . . , 

1 - 14 X + 172 x2 - 2308 x^ + 24572 x* - 244840 x^ + 2426928 x^ + . . . 

Another HeunG{a, q, 1, 1, 1, 1; x) example, such that the two parameters a 
and q correspond to fixed points of the symmetries (139) of these particular HeunG 
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functions, read: 
HeunGl^ , , 1, 1, 1, 1; 9- 

= 2Fi([^, [1]; 81- .T- (1 -27x +243x2; 
1 ^ /,1 1, 1 - 27a; + 243a;^ 

2- 



X \ 6 6 ^ ' 729 a;3 

= 1 +9x +81x2 +567x3 + 729x4 -72171 - 1764909x^ - 28284471 
- 343842327 X® - 2859802119 x^ - 2072088459 x^° + 523309421259 x" 
+ 13407709577211x12+226522478442087x^3 +•••, (P-8) 

where the coefficients arc actuahy integers, their sign having a period six. This example 
is nothing but revisiting (H.31) with a x — >■ — 81x change of variable. 

Among the fixed points of the symmetries (139) the case a = for any value of 
q, corresponds to a special limit. Let us write q as q = (1 — r2)/4, one has 

Hma^oHeunGia, ^-j^, 1, 1, 1, 1; ^^ax) = 2^^i([^y^, ^^], [1], 2^ x), 

which is, of course, a series with integer coefficients for any even integer values of r 
(it is a simple polynomial expression for odd integer values of r). 

More identitiesf read: 

iJeMnG(-3, 0, 1/2, 1, 1, 1/2; 12 • x) ^ 2i^i([^, J], [1]; 108 • x^ • (1 + 4x)) 

, 1/9 /rl li r 1 108 • X • (1 +4x)2 

= (l-12x)-i/2 . 2J?i [-, -],[1]; - ^ ' 



.6' 3"'^" (l-12x)3 
1 +12x2+48x3 +540x4 +4320x5 + 42240 x^ + 403200 x^ 
+ 4038300 X** + ••• , 

Fe«nG(4,i,i,i,l,i;64x) = 
^ (l-16xrV..,^,M 



6' 3^' ^ ^' (1 - 16x)3 

s 1 /9 /rl 1-, r , 432 X \ 

= 1 + 8 X + 192 x2 + 6656 x^ + 275968 x^ + 12644352 x^ + 616562688 x^ 
+ 31366053888 x^ + 1645521666048 x® + 88371818921984 x^ + 



3F2 



i i 2 108x_ 

3'2'3^'^ ' ^' (l-16x) 

1 - 16x „ /A 12,,^ , 108x2 



X /A 12,, 

-• 3^2(13, ^.gUl, 1], 



l-4x 2' 3" (1 -4x)3. 

= 1 -12x -36x2 _ 1922:3 - 1380x4 - 11952x5 - 116928 x^ 
- 1242624 x^ - 14006628 X® - 164954640 x^ + ■ ■ ■ , 

f Not to be confused with Goursat-type identities on hypergcometric functions (see for instance 
section 4 in [155]). Here, the hypergcometric functions have the same parameters but different 
pullbacks. 
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2^1 ([2' 3]' W'^) ^ (1 -16x +16x2)-i/4 (p.ll) 

„ /r 1 5 , r 1 X • (1 — X) \ 

^^^^([12' 12]' W' -'''■w-rtrik^)- 

This globally bounded series becomes a series with integer eoefScients after the 
rescaling x — > 16a;. 

2Fi([J, J], [l],x) (P.12) 



(1 +3x)-i/^ • 2Fif[-^, [1], 27 



X • (1 — x) 



'12' 12^' ' " (1 +3a;)3 
This globally bounded series becomes a series with integer coefficients after the 
rescaling x 64 x. 

2^^if[i |, [l],x) (P.13) 



3 3' 



1/4 /I 5 , , , x^ -(l-x) 



(9 -8x) 

This globally bounded series becomes a series with integer coefficients after the 
rescahng a; — > 27a;. 

A few hypergeometric functions occur in the analysis of the Yang-Baxter 
integrable hard-hexagon [156, 157, 158]: 



,5 



/r 1 5 , , , X l-f 11a: - a;^ \ , ^ 

2-F^i(hT, 1 ; -1728 5^ ^ ^) P.14 

V42'12J'^J' (l- 228x + 494a;2 + 228a;3-f a-4)^/ ^ ' 

1 - 228 a- + 494 x^ + 228 a;^ + a;^ \ i/4 

a:4 - 12a:3 + 14a;2 -I- 12x+ 1 / 

1 5 , , , a;5(l-^lla:-x2) 

X 2Fi{ —, — , 1 : -1728 !^ ^ 

V12'12J'^J' (x4- 12a;3 + 14x2 + 12a; + l)' 

and the globally bounded algebraic hypergeometric functions 2-f'i([l/6, 2/3], [1/2], a;) 
2Fi([l/4, 3/4], [2/3], a;): 

2^i([^, ^], [^]; 27a;) = 1 + 6 .a; + 105 a;^ + 2184 a;^ + 48906 a;^ -f 1141140 x'^ 
6 3 2 

+ 27335490 a;^ + 666865800 x^ + 16488256905 x^ + ■■■ (P.15) 
2Fi([i, ^], [^]; 64x) = 1 + 18x + 756.x2 +37422^3 -f 1990170x'' 

4- 110198556 x^ 6261870888 x^ + 362293958520 x^ +■■■ (P.16) 
Note that the two second-order operators 

, 1 (llx-3) 1 1 



6 X • (x-1) 9 X - (x-1) 

, 2 (3x-l) 3 1 , , 

3 X • (x — 1) 16 X • (x — 1) 

annihilating respectively the two algebraic functions 2-F'i([l/6, 2/3], [1/2], x) 
2^1 ([1/4, 3/4], [2/3], x), are such that their symmetric sixth power are (non-trivially) 
homomorphic, the intertwiners being order-six operators. 
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